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ABSTRACT 

This publication vas produced as a teaching tool for 
college c'-3rtistry* The book is a text for a computsr«based unit on 
the chemistry of acid-base titrations, and is designed for use with 
FOETRAN or BASIC. computer systems, and with a programmable electronic 
calculator, in a variety of educational settings. The text attempts 
to present computer programs that are relatively free of reliance on 
specialized large computer systems programs. The case^study approach 
presented is highly research and laboratory oriented. Similar subject 
matter is conventionally taught in most introductory college 
chemistry courses, but this text material attempts greater depth of 
instruction through utilization of the computational resources of a 
computer. (Author/BT) 
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ABSTRACT 

The 223-page book, Case Studies in Systems Chemistry : Carboxyllc 
Mil Equilibria » was written, published, and distributed for evalu- 
ation as a teaching tool for college chemistry. The book Is a text 
for a computer-based unit on the chemistry of ac^d-base titrations, 
and Is designed for use with FORTRAN or BASIC computer systems, and 
with a programmable electronic calculator. In a variety of educational 
settings. The text Is novel In that It attempts to present computer 
programs that are relatively free of reliance on specialized large 
computer systems, programs that can be used by individual students in 
an open-ended manner. This case-study approach is highly laboratorv 
oriented, with a strong research flavor. Similar subject matter is 
conventionally taught in most Introductory college chemistry courses, 
but this text material permits significantly greater depth of 
instruction through utilization of the computational resources of a 
computer. 
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INTRODUCTION 



Collegiate Instruction In chemical equilibrium and chemical 
v-lnetlcj, the two broad areas of systems chemistry, has been limited 
In sophistication by the formidable barrier Imposed by the repetitive, 
time-consuming, and boring arithmetic required for comparing data from 
experiments with predictions of theory. Comparison of data and theory- 
a central aspect of the appreciation of science as well as a critical 
aspect of the learning of science— has thus traditionally been given 
scant attention In Introductory college courses dealing with the 
subject matter of systems chemistry. The wide availability of com- 
puters and computer terminals on college campuses has vastly expanded 
possibilities for exploration of the relationships between laboratory 
data and theory in systems chemistry, although these possibilities 
have not yet been exploited. This research project was designed to 
..ake available, to both instructors and students, resources needed for 
utilization of their available computer facilities for college-level 
instruction in systems chemistry. 

Work has progressed in the development of case studies in several 
areas of systems chemistry, with open-ended descriptions of how com- 
puters can be used to facilitate confrontation of student data and the 
predictions of student-formulated models. A detailed case study in 
the field of the solution chemistry of acid-base equilibria, appropri- 
ate ror instruction in an honors section of general chemistry, has 
been developed. as a prototype for other similar text materials pro- 
jected for coupled equilibria, kinetics, and simultaneous kinetics and 
equilibria. The coup let e case study, entitled Carboxylic Acid 
Equilibria, is submitted as a portion of this final report. 

Of strong interest to the investigator has been the presentation 
of this material in forms applicable to a wide range of computer 
systems. The case study includes programs written especially for an 
IBM 1130 system with a plotter, using FORTRAN language, and also 
related programs written especially for a remote typewriter terminal. 
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using BASIC language. A program for a Wang 700 series programmable 
calculator Is also Included, but our efforts to use such small 
computers has been disappointing and we do not feel that continued 
work with programming our Wang 700 or our Olivetti P602 Is warranted. 

The most Ingjortant aspect of this research Is to alert chemical 
educators to the new possibilities for teaching what has become a 
classical subjecc in general chemistry, analytical chemistry, and 
physical chemistry. The investigator presented a talk on systems 
chemistry (see Appendix A) to the New England Association of Chemistry 
Teachers at their August 1972 meeting, and copies of this talk were 
distributed to all participants at the 1972 Mount Holyoke Conference 
on Chemical Education sponsored by the Division of Chemical Education, 
American Chemical Society. It was announced at these meetings that 
materials being prepared under this grant would become available 
during the following year, and could be sent to chemlGtry teachers 
upon request. 

rffiTHODS 

An effective way of demonstrating to teachers of chemistry that 
a new approach to teaching can be used is to present examples of how 
it is being used. Likewise, one effective way of showing students how 
to use a technique is to show it being used, and then ask the student 
to use the technique in a different but closely related manner. For 
both reasons, the case study approach has been utilized in this 
research to provide examples of the use of computers to facilitate the 
coiiqjarison of models and data. 

Each computer program is illustrative of a wide variety of 
similar programs, so that the program need not be followed as a 
recipe, but can be varied and extended by the student. The, programs 
of the case study are meant to open new vistas, not to exhaust possi- 
bilities. A wide range of programs utilizing small desk-top program- 
mable electronic calculators were investigated, with detailed 
programming attempted with an Olivetti P602 and a Wang 700. One Wang 



program is presented in the con?)leted case study, but our efforts with 
four experienced student programmers and our attempts at stimulating 
student usage has been disappointing and generally negative. Workable 
programs for these calculators can be written, but there is an inherent 
lack of flexibility that inhibits student experimentation with these 
programs. We have found great enthusiasm for the IBM II30 programs, 
particularly with graphical plotter output, and good response to the 
BASIC programs, contrasting markedly with the responses to the smaller 
and less versatile equipment. 

The original programs were written by the students, in many cases, 
and many modifications of the programs have been suggested by students. 
The materials of the case study have been tested by classroom- laboratory 
use in a college general chemistry course for students with strong 
science preparation, and in an independent study project by a secondary 
school senior who was con?)lecing a second year of chemistry. Most of 
the students involved in formulation and testing of this material are 
girls, and almost all are science majors. 

The completed case study, Carboxylic Acid Equilib ria, has been 
?ent to 30 interested high school and college teachers who have 
expressed interest in this research, and some suggestions have been 
received. Corrections and changes arising from thelr^omments will be 
incorporated in a set of revision sheets that are toj^ incorporated 
in the case study in June 197^. Preliminary responses indicate that 
the computer prog ams are indeed usable in computer installations 
other than the Smith College facilities. 



RESULTS 

The case study begins with a detailed statement of its perform- 
ance objectives. Three ge^ieral objectives are stated, followed by a 
list of specific objectives written in behavioral terms. The 
Investigator has used this case study for teaching a unit in general 
chemistry, and has evaluated student performance in terms of the 
ability to meet the stated general objectives. Such student evaluation 



Is highly satisfactory with an honors-pass-fall (or Its equivalent) 
grading system. The performance tends to be either full achievement 
or failure, and the Investigator sends students nith a performance of 

failure - back to the text, lab, and computer for further work. 
Hedlocre work gets equated with failure, and the middle range of 
grades gets eliminated. The acceptance of such an evaluation scheme 
probably carries with it a rejection of a bell-shaped distribution 
of grades. 

The case study is an Integrated discussion of theory and methods 
of treating laboratory data. Computer programming details are intro- 
duced as needed. The approach is intended to appear highly pragmatic: 
methods are introduced whenever there is a need. The goal set before 
the student is the arranging of meaningful confrontations between 
data and theory. 

Eng)hasis is placed on model building and on simulation of these 
models. Students are reminded that a computer gives a novxce chemist 
power to test models that was unavailable to researchers in former 
decades. Students are encouraged to test assertions that appear in 
older research articles, as well as assertions made by textbooks, 
teachers, and this case study. Some teachers may find this open 
approach threatening, as indeed may some students. But the approach 
offers an opportunity for students to become liberated from reliance 
on textbook dogma, and to test important ideas for themselves. 

Distinction is made repeatedly between macroscopic and microscopic 
descriptions of chemical systems, emphasizing that macroscopic 
descriptions are phenomenological and therefore potentially accessible 
in terras of experimental observables, whereas microscopic descriptions 
are necessarily made in terms of theories and thus are accessible only 
in terms of models. Perhaps the discussions add motivation for the 
students to achieve facility with simulation techniques, for simulation 
of models provides an easy method of dealing with microscopic 
descriptions. 

The chemical literature contains vast numbers of research reports 
dealing with acid-base equilibria with only a few containing as much 



8 



detailed quantitative analysis as in this case study. There is thus 
2 reservoir of material in the Journals for use in teaching about (and 
in learning about) the equilibria involving carboxylic acids. Early 
in the case study (on page 61), there appears a listing of biblio- 
graphical sources providing entry for a student into the primary 
research literature. Additional comments are made throughout the text 
to lead the student to the library. 

The initial computer programs are presented in FORTRAN language, 
with much of the output to be obtained with a plotter. This is by 
far the most satisfactory method that has been employed by the 
investigator. However, - equate results can also be obtained using a 
time-sharing terminal and employing BASIC language, and several 
appropriate BASIC programs are presented, beginning with page 90. 
In every case, effort is made to show a pro£ram that "works", but at 
the same time to encourage the student to adapt that program to the 
particular specialized needs of the moment. 

In teaching this material, the investigator has repeatedly found 
that students miss many important aspects of both theory and experi- 
ment when they do a single experiment with a monoprotic acid. Thus 
repetition (by titrating a diprotic acid and interpreting that new 
data) of the whole procedure has a great deal of pedogogical value. 
Happily, there are many significant features that are new to the 
diprotic acid case, thus adding some additional interest. The new 
features center largely about the use of both macroscopic equilibriurc 
constants and microscopic equilibrium constants to describe the same 
chemical system, and the possibility of having alternative but indis- 
tinguishable descriptions of the same system. These features of acid- 
base equilibria are explored in Chapter 2. Finally, in Chapter 3, 
polyprotic acid equilibria is introduced in a research context, leaving 
many possibilities open for the inquiring student and the imaginative 
teacher. 



9 

-8- 



CONCLUSIONS 

It is clear that computer use by students in an introductory 
general chemistry course can be used to integrate laboratory and cla. 
room, as evidenced by the accompanying textual materials and by their 
successful use in the investigator's classes. It is the Judgment of 
the investigator that satisfactory evaluation of student performance 
can be made by using the three general objectives stated at the 
beginning of the case study. We feel that computer systems comparable 
to the IBM 1130 with plotter provide flexibility and performance that 
is completely adequate for the instruction envisioned with this case 
study, but that further efforts to utilize programmable calculators 
such as the Wang 700 or the Olivetti P602 are not warranted. Undoubt- 
edly there should be continuing efforts to increase the portability 
-.nd transferability of computer programs, and other writers should be 
encouraged to develop better ways of i^riting programs that can be ru: 
on a variety of different computers with a minimum of editing. 

The con?)leted case study is evidence that the availability of 
computing facilities for students can change the substantive content 
of introductory chemistry courses. The wide distribution of such case 
studies can be an aid to stimulating similar activity on many campuses, 
and to exchanging curriculum information between teachers. The case 
study approach also provides a medium for self study in which a 
student, with modest laboratory equipment and access to a computer, 
can learn a substantial amount of science without requiring extensive 
interaction with an instructor. 
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Appendix A 

Now Kuglnm! A j« :-.oo i nt. { on oi" ChemiHtry ■roncl-.or:-. : : Aun";*' <• , 107. • 

SYSTEMS CHEMISTRY: putting it all together 

George Fleck 
Department or Uh.-mistry 
Smith College 

Northampton, Massachu-.etts OIO60 

Humpty Vumpty 6 at on a. voa.ll 
Hampty Vumpty had a Qnzat iall) 
Kit the. hofUz^ 
And all thz king' 6 me.n 

Couldn't put Humpty Vumpty back togzzhzn again. 

Poor Humpty Dumpty got taken apart rather violently. It's 
easy to smash an egg. With some practice, an egg can be disassem- 
bled with more finesse. The art of separating yolk and white is 
not so hard, unless you insist on using just one hand. And with 
some more sophisticated techniques, the egg can be separated with 
finer resolution and a variety of pore compounds can br isolated 
from an egg. To take apart an egg — even to take apart an egg and 
separate it into many of its constituent compounds — is a straight- 
forward task. But clearly more than a kirg's command is needed to 
put it all back together again. 

So it is with descriptions of an egg. Analysis — the description 
of an egg in terms of its constituents — is straightforward. We 
can find out what is in the egg, and we can characterize each of 
the components in terms of its molecular structure. But synthcs is 
— explaining how the molecular constituents function together as 
a system to create a Plymouth Rock chick — is an intellectual task 
of considerably greater magnitude. Whether literally or figuratively, 
it is a whole lot easier to take an egg apart than to put it back • 
together again. 

There is a lot to be learned about things by taking them 
apart. The search for the ultimate parts and pieces of things and 
stuff was a driving force in transforming alchemy into chemistry. 
When earth, air, fire, and water were replaced by scores of chemical 
Clements, chemistry had become an analytical science. And when 
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l\'iRC Two 

the elements had been taken apart and found to be made of elec- 
trons, protons, and neutrons, part of chemistry had almost be- 
come physics. An enticing goal for many chemists has been to 
talk about the world in terms of simpler and simpler models, 
and in the process to cut up chemical reality into smaller and 
smaller pieces. 

Yet the rise of chemistry during the past hundred years as 
a significant force in our society has not appeared to the lay 
observer to be due to progress in identifying elements, nor to 
the insights of quantum mechanics, nor to the ability to deter- 
mine molecular structures, nor to the ability to make refined 
chemical analyses. Layman, physician, businessman, homemaker, 
soldier, and politician ha^'^e seen and felt the impact of chemistry 
in the twentieth century in terms of a host of new materials put 
together by the chemists: fabrics, drugs, plastics, explosives, 
fertilizers, and construction materials. There is private fun 
in taking things apart, but there is public reward in putting 
things together, especially for building — that is, in synthe- 
sizing — things new and useful. 

As we look to the future, we Cjan see new chemical challenges • 
Design and synthesis of new materials will continue to be important. 
Society will demand from chemistry new fuels, new fabrics, and new 
foods at low cost and in abundance, without regard to constraints 
imposed by thermodynamics. Society will ask chemistry for ways 
to perform miracles in waste disposal, air and water purification, 
and instant drug therapy for cancer, heart diseases, and olf age. 
In a world rapidly moving toward a confrontation with the limits 
of readily available natural resources, many of the material 
expectations of society may be found to be unrealistic or ^.nattain- 
able, and chemistry may soon be called upon to help delineate the 
' imits of the chemically-possible exploitation of this planet. 
Increasingly, chemists may expect to be dealing with more compli- 
cated and less well defined chemical systems: the various eco 
systems, for example, and living organisms, for another. Chemists 
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Page Three 

will also I'ind themselves talking to oconoinisr'? , politicians, 
voters, and consumers, trying to explain to this public tln« 
limitations placed on chemical systems by the realities oi" basic 
physical laws. 

In the relationships of , . ;y with society, the future 
will require increasing unders eluding of complex chemical systems. 
Chemists will need to know what to expect when lots of interacting 
molecules get put together. Chemical education should begin now 
to provide students with the intellectual tools required for 
coping with complex systems. 

TAKING IT APART/ VERSUS PUTTING IT TOGETHER: 
IMPLICATIONS FOR THE CHEMISTRY CURRICULUM 

The word system occupies a central place in the teaching of 
physical. chemistry , and the word is commonly used in all areas 
of chemistry to denote the mixture of chemicals being studied in 
an experiment. The need is to make the concept of system an 
open-ended concept, capable of being extended by the student to 
systems of greater and greater complexity. There needs to be 
explicit consideration of the ways in which sub-systems combine , 
for seldom are the properties of a chemical system the simple sum 
of its parts. There needs to be exploration of the advantages 
and disadvantages of the phenomenological description of a systcr, 
and of the relationships between macroscopic and microscopic 
variables and models ii. describing systems. Finally, there nccv. ; 
to be an awareness throughout chemistry that students should be 
reminded that the real world is made of chemical systems, few of 
which are understood, many of which need to be understood. 

We began to make substantial progress in getting students 
involved, at the freshman level, with systems chemistry when an 
IBM 1130 computer became available on campus. For our teaching 
of solution equilibria and kinetics, access to a computer with 
a plotter revolutionized our approach. The mathematics has 
shifted from calculus to arithmetic and simple algebra. The lab 
and classroom have come closer together, as the computer has 
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Page Four 

facilitated a teaching strategy in which students confront their 
experimental data with predictions o£ their models. By giving 
the students the powerful tool of computer simulation of a model, 
the role of the teacher has changed; since the students have the 
tools to challenge ideas and models presented by the instructor, 
it is possible to have students and instructor joining forces in 
an active pursuit of models that are faithful to good experimental 
data. Finally, the use o£ simulation techniques means that the 
students can discover and verify certain properties of systems, 
without having to accept textbook dogma. We find that realistic 
small research- type projects are quite appropriate at the college 
freshman level, at least in small classes with well -prepared 
students . 

Systems chemistry, when taught with the support of computer 
facilities and a well-staffed laboratory, can enter the che -xstry 
curriculum in various ways. We are struggling with the appropriate 
places where these concepts ought to be taught. They can be 
taught very early in a student *s career, at the point where we 
usually teach stoichiometry and titrations. It can come at the 
point where we usually talk about entropy and free energy. It 
can wait until senior year, and be used as an integrating course 
to help put it all together for a student. For the remainder of 
this discussion, I wish to dwell on the substance or the concepts 
that need to be taught, without necessarily implying the courses 
where the concepts might be introduced or reinforced. 

CONCiiPTS OF SYSTEMS CHEMISTRY 

LEVELS OF CONSTRAINT ON A DYNAMIC CHEMICAL SYSTEM . Every 
reacting system, given a constant environment, is predestined to 
a particular evolution and a particular equilibrium state. In 
the absence of an adequate reaction mechanism, few predictions 
can be made about the initial time-course of the syste^f. The 
passage of time, though, places constraints on any reacting sy^tam. 
As time goes by, less information is needed to formulate a 
description of the further evolution of the system. By the same 
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tckon, less information about the system is available from cluta 
obtained from a system close to equilibrium. 

A generally -reacting system can be described from time- zero 
to time-infinity by postulating a reaction mechanism, and writing 
a mathematical model that includes a set of mass-action rate 
equations. If values of the rh!:e constants are known, the system 
can be simulated by approximating the derivatives in the rate 
equations by ratios of small differences, and the time-dependence 
calculated by repetitive arithmetic via a process that is quite 
appropriate for computer computations. Experimental data can be 
used to evaluate rate constants-, but experiments must be carefully 
designed for that purpose. 

For a complicated mechanism, the early tim? phase may be 
categorized as "transient behavior*' and ignorec, so that attention 
may be focused on the later stages of the reaction. For any 
syhtem that can be characterized by mass -action rate equations, 
eventually all transient behavior passes by, and the final approach 
to equilibrium is a simple, exponential, first-order relaxation . 
With some systems, the reaction all (or virtually all) the way 
from time zero proceeds by one or more first-order relaxations. 
The time-constant for a relaxation (a macroscopic rate constant) 
is a function of all the microscopic rate constants of the 
mechanism. For a complicated system, the approach to equilibr i 
in a final relaxation is described by a single number. This is 
a great simplification for macroscopic or phenomenological desc r i pt i . i . 
The situation renders hopeless, however, the prospect of gaining 
detailed mechanistic information about a complicated reaction from 
data obtained near equilibrium. 

Often, but not always, a system may react in such a way so 
that after an initial transient phase, some of the components may 
be in a steady-state . Then just a few chemical species are bein^ 
transformed from reactant to product, while other species have 
concentrations that are time-invariant. Steady-state systems can 
be described with fewer variables, but again steady-state rate dat;i 
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Six 

yiolil loss information, than a generally react ing system. Most 
bioiojiical systems are in steady-states. It would seem that 
general analysis of steady-state systems must be an important^ 
part of a systematic study of systems chemistry. 

Finally, there is the equilibrium state# A system that is 
actually at equilibrium has reacted as far as it can react. Such 
a system can be described with the minimum number of variables. 
However, the cost of obtaining this simplicity of macroscopic 
description is that equilibrium experimental data has essentially 
no information content regarding mechanisms of reaction. 

STABILITY OF A REACTING SYSTEM, Many thermodynamicis ts speak 
with assurance when they say that no reacting system overshoots 
equilibrium, that entropy always increases during a reaction, 
and that the driving force for achievement of the equilibrium state 
is the derivative of the free energy with respect to the degree 
of advancement. Each of these statements may well be true, if 
properly applied and if properly qualified. Yet apparent contra- 
dictions abound. The concentration of an intermediate often 
overshoots its equilibrium value during the transient phase of a 
reaciion. Sustained oscillations have been observed in several 
reacting systems. 

Thus it is especially instructive to explore some of these 
general statements about the one-way nature of the approach to 
equilibrium and the staoility of the equilibrium state. A computi r 
simulation is an ideal method for the student to conduct this 
exploration. Along with the output of concentration versus time, 
one can calculate entropy and free energy versus time. One can in 
fact look at entropy and free energy of individual species, as well 
as for the system as a whole. A student should then formulate 
some of the general criteria for reacting systems in terms that 
can be defended on the basis of simulations of reaction mechanic 

Most elementary chemical reactions have inherent stabilitv 
because of negative feedback built into the form of the mass-action 
equations. Some chain reactions appear to lack this stability, 
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Page Seven 

and it is instructive to pinpoint wh-ere stability is lost, aiui how 
it can be restored. Such considerations of stability are very 
important in talking about control in a living organism, wIkmc 
steady-state concentrations are regulated by molecular control 
processes that need to be understood. Using simulation methods, 
an undergraduate student can subject some of the qualitative ideas 
that appear in the biochemical literature to some critical tests. 
The concept of self -regulated dynamic stability is really quite 
apoealing. 

STRUCTURE WITHIN A REACTING SYSTEM. Homogeneous solutions 
are convenient media for performing reaction-rate studies in the 
laboratory. Homogeneous systems are also particularly amenable 
to modelling with mass-action rate equations. However, there is 
evidence that suggests that -some types of interesting reactions 
cannot occur in a strictly homogeneous medium. Certainly, a living 
cell is a structured medium, and many of the special properties 
of a cell derive from this structure. There is continuing interest 
in ways in which an initially homogeneous solution, or gas, or gel 
can achieve structure. It is not enough to divide the woi .d into 
homogeneous and heterogeneous media, for the criteria for structure 
are probably quite subtle. The dynamic stability of non-uniform 
structures such as flames or living cells is intimately involved 
with reaction processes. Prototype models can be cast into 
mathematical form, and simulated, setting up open-ended investi^j 
tions for students". 

SOME FUNDAMENTAL CONCEPTS. There art a few ideas that are 
central to systems chemistry that need to be stressed to students. 
I shall merely list a few of the ones that I feel are significant. 

1. The macroscopic behavior of a multi-component equilibrium 
system, or of a multi-reaction reacting system, can be described 
with fewer numbers than the microscopic description of individu.'.l 
molecular processes. 

2. A system is the superposition of individual sub-systeiii . > 
but the properties of the whole are not the additive sum of the 
properties of the parts. 
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Page Eight 

3. There is room for subjectivity and creativity in making 
models of chemical systems. Various models may be equally useful, 
and equally faithful to experimental data. 

4* The best that we can sometimes accomplish is to make 
statements about the forms, not the details, of equations that 
govern a system. But sometimes it is the structure of the dynamic 
system, not necessarily all the details, that is significant. 
Thus, although probably every candle flame is different, there arc 
common features of candle flames that arise from the structure 
of the mathematical model, features that yield stability and 
integrity to a candle flame. 

SUMMARY 

It has been asserted that chemistry makes its impact on th^ 
world by processes of synthesis, by putting things together. 
It has been predicted that this putting- together activity will 
increasingly involve chemists in dealing with complex systems, 
chemical in nature, that are of concern to a wide range of non- 
chemists. To deal with such research concerns in the future, as 
well as to make current instruction relate to issues of sccicty 
as they impinge on chemistry, it seens appropriate to introduce 
some explicit systems chemistry into the chemistry curriculum. 
It is suggested that computer simulations of equilibrium syst< r. . 
and kinetic systems is an effective way to accomplish this, wii ^ » 
the present structures of chemical curricula. 
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INTRODUCTION TO CASE STUDIES IN SYSTEMS CHEMISTRY 

The focus of each of these case studies is a chemical system. 
We shall use investigations of chemical systems to gain experience 
with some useful experimental techniques and theoretical methods. 
We shall also use these investigations to explore certain central 
philosophical issues, including 

inter-relationships between facts and theories in science; 

analytical versus synthetic approaches to learning about nature; 

inter-relationships between macroscopic and microscopic descriptions; 

the role of time in describing nature; and 

the possibilities for subjective, creative expression in science. 
The systems themselves are worth studying for their own sake. 
By pointing directly to significant and interesting chemical systems, 
we shall try to keep tangible hunks of matter — a living cell, a 
lake, a beaker of solution, a jug of wine — in front of our eyes and 
at the front of our minds. And then, to develop realistic theories 
that are clearly and directly related to the real chemical systems, 
we shall adopt a strategy of con^Kontation oi data invm chmicjodL zxpMMmntA 
uUXh thz p^zdaitiom oi chmaial thzoKy. Over and over again, we shall 
set up confrontations by carefully designing experiments and by 
carefully designing chemical models so that both experiment and 
model yields a number that should be the same li experiment and model 
are consistent. We shall exploit the use of computoA ^AjnuZdUjon meXhodt 
for obtaining numerical predictions from our chemical models • 
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SOME COMMENTS ABOUT THE USE OF COMPUTERS 
IN THE INVESTIGATION OF SOLUTION EQUILIBRIA 

Much of the material that follows is oriented toward 
computer calculations. There is enough computer output presented 
so that a student can read about the chemistry, with understand- 
ing, without having access to a computer. However, the case 
study is primarily intended to involve each reader with calcula- 
tions and with interpretation of his or her own data, and for 
these purposes a computer is of great utility. 

This case study attempts to present the computer as a tool 
that can liberate the student from dependence on the instructor 
and on the textbook. The case study is also viewed as a resource 
for the student in developing a mature understanding of the 
relationships among hypotheses, experimental data, and interpre- 
tation in the special field of solution equilibria. Because the 
computer frees both instructor and student from the limitations 
of paper-and-pencil calculations, it is possible to explore in 
detail and with assurance the consequences of various chemical 
issumptions about the systems being investigated. 

The computer is not viewed as a substitute for laboratory 
experience. On the contrary, the author has found that extensive 
use of the computer has made possible a chemistry course in which 
laboratory experimentation and the use of chemical theory are 
very closely integrated. Throughout the case study, there is an 
emphasis on the confrontation of the student* s own data and the 
quantitative predictions of the student •s own model. 

In the main body of the case study, use is made of FORTRAN 




programming language, and the illustrative computer output has 
been obtained from an IBM 1130 computer system which includes 
an IBM 1627 strip-chart plotter. The discussion assumes that 
input of programs is via decks of punched cards. Our experience 
in chemistry instruction at Smith College has been with hands-on 
student use of the computer system. Some very useful plotting 
subroutines, available on our computer, are described in Appendix 
I. A single 15-minute orientation session, given to groups of 
about ten students at a time by an undergraduate computer-center 
assistant, is sufficient to provide all the information needed 
for running a program. 

In writing this case study, as well as in teaching this 
material to entering college freshmen, it is assumed that the 
student has no previous knowledge of computers or of programming 
languages. Working programs are presented that can be modified 
easily by the changing of just a few statements. As the case 
ftudy (or the semester) proceeds, the programs become somewhat 
more complicated, with the result that most students gain a 
substantial degree of facility with the computer while learning 
about the chemistry of equilibria in solution. Increasingly, 
students are entering college with some experience in FORTRAN or 
BASIC, and these students are encouraged to write their own 
programs and to attempt special computer projects. 

At the end of each chapter is a section discussing the use 
of the BASIC programming language. The illustrative computer 
output has been obtained from a time-sharing terminal of the 
University of Massachusetts computer system. We use a Cartertone 



S15C data terminal which is a modified IBM Selectric typewriter 
coupled to the University of Massachusetts computer center via 
a Bell System DATA •Phone. 

All FORTRAN programs in this case study are identified by 
number, beginning with FORTRAN Program 1. All BASIC programs are 
identified by letter, beginning with BASIC Program A. 

All Wang programs in this case study are identified by lower 
case letter, beginning with Wang Program a. We have attempted 
to use a Wang 700 programmable calculator with output presented 
on a Wang flat-bed plotter, but have found that the most con- 
venient usage by students in studying, solution equilibria is 
by numerical output. We have therefore given a program in this 
case study only for use with a Wang unit with just nixie-light 
output. 
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PERFORMANCE OBJECTIVES 
FOR THE CASE STUDY 

If f ^ 

GCMERAL Ol'JECTJVliS 

K You should be able to shov; how tlic predictions of ;i 
postulated set of chemical equilibria can be confronted vjitlj 
the quantitative facts from chemical titration experiments in 
such a way that the cher^iical model miQUt be found inconsistent 
with the experimental data, 

II. For a postulated set of chemical equilibria involving 
an acid with at least two dissociable protons, you should be 
able to show how numerical values for the equilibrium constants 
for that acid can be obtained from titration data of pH versus 
volume of added titrant. 

III. For a particular acid, you should be able to obtain an 
experimental titration curve, analyze that data, propose a set 
of equilibria to serve as a chemical model, perform a computer 
calculation of the titration curve predicted by that model, 
and finally compare the results of experiment and the predic* 
tions of theory. 

^ ¥ ¥ ^ V 

SPECIFIC OSJECnVBS 

To achieve the three general objectives, you should be 
able to start with a particular set of equilibria, postulated 
for a particular chemical system, and deal with that chemical 
model and the chemical system in a stepwise fashion as out- 
lined in the list of specific objectives that follows. Given 
the particular diprotic acid, you should be able to: 

1. STATE A CHEMICAL MODEL. Make a list of all the 
chemical species that you decide should be assumed to exist in 
the solution. V/rit3 a balanced chemical equation for each 
equilibrium reaction that you decide to assume. Check to be 
sure that there is at least one xvay, via the assumed chemical 
reactions, for each assumed species to be produced from the 
chep.icals that are mixed together to prepare the solution. 
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2. STATE THE MATHEMATICAL MODEL. Write a set of simultaneous 
algebraic equations as follows :^ 

a. For each chemical reaction, write an algebraic equilib- 
rium-constant equation in terms of molar concentrations. 

b. Write the electroneutrality equation. 

c. Vrite algebraic conservation equations in terms of 
molar concentrations, 

3. DESIGN A TITRATION EXPERIMENT. ProposeVan experiment for 
obtaining experimental data sufficient to pi t a titration curve. 
Decide on appropriate concentrations. VJhich solution will you 
use as titrant? 

^. OBTAIN AN ALGEBRAIC PREDICTION FOR THE TITRATION CURVE. 
The simultaneous algebraic equations that constitute the 
mathematical model should be combined to give a single equation 
that involves (H*) , the equilibrium constants, and the known 
total conczntJtation6 (the macroscopic concentrations or analytical 
concentrations) of acid and base. Then change to volume variables, 
obtaining an equation that relates the two experimental variables 
(H*) and V^^^^^^^. 

5. ESTIMATE VALUES OF THE EQUILIBRIUM CONSTANTS. Calculation 
of the theoretical titration curve predicted by your model 
requires that you have numerical values for each of the equilib- 
rium constai!ts. You may use a literature value for the ion 
product of water. Probably the most convenient assumption for 
the various acid dissociation constants involves use of the 

half -equivalence method. But however you proceed, you must put 
some numbers into the calculations as starting values, numbers 
that later can be verified or modified as calculations are 
compared with data. 

6. PERFORM A COMPUTER SIMULATION OF THE SYSTEM. Using the 
algebraic equation that you obtained from the mathematical model, 
and usinj the estimated numerical values of the various eqtiilib- 
rium constants, use a computer to calculate a predicted titration 
curve. Do your data agree with the Calculations? If not, then 
attempt to get a good fit by making appropriate changes in the 
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values of the equilibrium constants used in the simulation. By 
a series of successive approximations, you may be able to get a 
2ood fit in the regions of the titration curve where the shape of 
the curve is highly dependent on the values of the equilibrium 
constants. It may be that your data are not accurate or precise 
or reliable enough for this comparison; if this is the case, you 
must return to the laboratory and get some better data. Your 
results should include a set of values for the equilibrium 
constants, an estimate of the reliability of each constant, and a 
judgment of the extent to which the data are consistent with the 
original chemical model. 



HELPS ^OR ACHIEt/IMG THE OBJECTIVES 

General objective I, the confrontation of the chemical model 
with experimental titration data, is achieved in specific objective 
6. General objective II, evaluation of the equilibrium constants, 
is achieved in specific objectives 5 and 6. 

Key to Textbook: 

G. M, Fleck, Eqailibfila In Solution, -Jew York: Holt, Rinehart 
and Wins con, Inc., 1966. Objectives 1 and 2: chapter 3. Objectives 
3 and 4: chapters 4 and 5. 



AN EXAMPLE OF THE CONFRONTATION OF THEORY AND EXPERIMENT: 

NUMERICAL EVALUATION OF THE EQUILIBRIUM CONSTANT 

FOR THE DISSOCIATION OF A MONOPROTIC ACID IN SOLUTION 

During the course of the titration of an aqueous solution of 
a monoprotic acid with a solution of sodium hydroxide, t?ie pH value 
(measured with a glass-electrode — calomel-reference-electrode 
assembly immersed in the solution) increases. Initially, before 
any sodium hydroxide titrant has been added, the pH is the pH of 
a solution of the pure acid in water. At the equivalence point, 
when equal numbers of moles of acid and base have been added to 
the titration vessel, the pH value is the same as the pH of a 
solution prepared: by dissolving the sodium salt of the acid in 
water. Past the equivalence point, the pH continues to rise as 
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additional titrant is added, approaching the pH value of the 
sodium hydroxide solution used as titrant. 

The titration ca^ve is a plot of pH versus Vj^^^y, the volume 
of added sodium hydroxide titrant solution. We shall examine a 
way to use these titration data to obtain a numerical value of 
the equilibrium constant for the dissociation of a moncpKOtic 
acid, an acid with just one dissociable proton. We shall do this 
by suggesting a chemical model to describe (and perhaps even 
partially to explain) the molecular equilibrium within the 
solution. Then we shall perform a computer simulation set up in 
such a way as to confront the model and the data. Demolition of 
the model is possible as a result of this confrontation. Also 
possible is the exposure of bad data. In this discussion, we 
shall follow an or^^anization suggested by the acid-ba6^ zqultib^la 
specific objectives 1-6 for a diprotic acid. 

Let us consider the titration of 0.1000 molur acetic acid 
solution with 0.1000 molar sodium hydroxide solution at room 
temperature. 

OBJECTIVE 1. STATE A CHEMICAL MODEL. 



OH' 
H2O 

Equilibria assumed: 

HA # A" ♦ H* 

H2O 4? OH' ♦ H* 
OBJECTIVE 2. STATE THE MATHEMATICAL MODEL, 
a. Equilibrium Constant Equations: 



List of chemical species assumed: 




(A')(H*) 



K » 

(HA) 



(i) 



- (OH')(H*) 



(ii) 
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where (A"), (HA), (H*) , and (OH") are individual-species concen- 
trations (microscopic concentrations) in moles/ liter. 

b. Electroneutrality Equation* 

(Na*) + (H*) » (OH") + (A") (iii) 

c. Conservation Equations: 

[A] = (HA) + (A") (iv) 

[Na] » (Na*) (v) 
where [A] and [Na] are total concentrations of all acetate-con- 
taininf species and of all sodium-containing species. 

OBJECTIVE 3. DESIGN A TITRATION EXPERIMENT. The goal of this 
experiment is to obtain a set of ordered pairs {(H), Vjj^qjj} for 
representative points in the titration, all the way from the start 
to past the equivalence point. Note that the experimental 
variables, as directly observed in the laboratory, are pH and 
V 

NaOH* 

OBJECTIVE OBTAIN AN ALGEBRAIC EQUATION FOR THE TITRATION 
CURVE. Equations i through v, a set of five simultaneous equa- 
tions in five unknowns, constitute the mathematical model for 
this acid-base titration. V/e seek to eliminate four of these 
unknowns and thereby obtain a single equation relating known 
quantities in [ I's to the measurable (H*). Equations i and iv 
can be combined to eliminate (HA) , {jiving 



[A] 



f(H*) 
+ 1 



(A") (vi). 
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Equations iii and v together yield 

(A*) - (H*) - (OH") + [Na] (vii) 

W 

!'/e now put everything together. Equations vi, vii, and also ii, 
are combined, resulting in the following cubic equation in (H ) : 

(H*)» + {[Na] * K}(H*)2 

* {K[Na] - K[A] - K^1(H*) - K^K = 0 ^^"^^ 

Our titration is performed by mixing two solutions: a 



solution of acid having concentration [A] « and volume Y^; 
and a solution of base havinf concentration [Na] = M^^ and volume 
added at any particular point in the titration equal to V^^, If 
there is additivity of volumes upon mixin,^ (that is, if ^goiution 
is equal to the um + ^Na^* then we can vrrite the two macro- 
scopic (total) concentrations for the solution being titrated as 

[A] . and [Na] « (ix and x) 

During the course of the titration, only (If*) and Vjj^ are 
variables. I'e shall nov/ substitute Equations ix and x into 
Equation viii co give 

%a - -^-^ (") 

where the fraction F is the following ratio of two polynomials in 
(H*): 

(H*)' * IC(H*)^ - {K^ + K-H^}(H*) - K^K 



F = 
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(H*)' * {K + Mjj^}(I!*)^ + {K-Mjj^ - K^Hn*) - K^K 
OBJECTIVE 5. ESTIMATE A VALUE FOR THE EQUILIBRIUM CONSTANT. 
Try the half -equivalence method for getting a :rial value of the 
equilibrium constant K to include in the calculations that depena 
on Bquation xi. The half -equivalence method asserts that 
(H*^) ^ K (and therefore that pH = pK) at the point in the titration 

at which v., has a value equal to just hati its value at the 

Ha -I ^ it 

equivalence point. Thus, half-way between the beginning of the 
titration and the equivalence point, the pH has a value that 
permits a good guess for the value of K. 

OBJECTIVE 6. PERFORM A COMPUTER SIMULATION OF THE TITRATION 
CURVE. Now use Equation xi to calculate a titration curve, the 
curve predicted ioK tfoaK ovon zxpe,Kime.ntaZ conditions t ^ox tjOiiK 
ac^cd, jJoA. tjOiiK thtimattd \)alut jjo/t K, expressed as a plot of pK 
versus Vj^^. The procedure is to assume a value for (H ) and then 
to use Equation xi to calculate the corresponding value of V,^^. 
If this volume is negative, the calculation is rejected as not 
corresponding to a point in the domain of chemical reality, and 



another value of (H*) is "picked. This process is repeated over 

and over, and the positive values of volumes are recorded. 

These renetitive calculations can readily be performed by 

a digital computer or a proorammahle calculator. A FORTRAN 

program that v;ill accomplish our task is jiven below: 

// JOB 
// FOR 

* LIST SOURCE PROGRAM 

*IOCSCCARD, 1132 PRINTER) 

C TITRATION CURVE ACETIC ACID 

REAL K, KW, NUM, MA, MNA 

VA » 25.00 

MA » 0.1000 

MNA » 0.1000 

K » 1.75E-5 

KW « 1.008E-ltf 

PH « 0.00 

1 PH * PH + 0.10 
H = 10**C-PH) 

NUM « H**3 + K*CH**2) - CKW+ (K*nA))*H-KV'*K 

DEN » H**3 + CK+MNA)*(H**2) + (K*MNA-KV;)*H -• KW*K 

VNA - -VA*NUM/DEN 

IF CVNA) 2, 20, 20 

2 IF (PH - 7.00) 1, 1, 30 
20 WRITE 0,10) VNA, PH 

10 FORMAT (F10.3, F10.3) 

GO TO 1 
30 CALL EXIT 

END 

// XEQ 

Code for Translating Program Into Algebra 

FORTRAN Symbol Algebraic Symbol 



K 



"a 

Na 
pH 

(H*) 

raised to the power of 



K 

K\-' 

VA 

VNA 

MA 

MNA 

PH 

H 
It 

** 
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The output of this program is a lonn tabulation of ordered pairs 

of numbers, where the first number is the volume of titrant added 

and the second number is the resulting pH value. As a sample. 

the first few pairs in the tabulation are j^iven below: 

0.027 2.899 
0.178 2.999 
0.337 3.099 ^ pH 
0.513 3.199 
0.71^ 3.299 
0.9«f9 3.399 

Of course, there are many other v;ays to display the results of 
the calculations, and you may have your own ideas about the best 
way to set up your confrontation of model and data. 

Compare calculated numbers with numbers from your own lab 
data book. If they don't match, try the calculations again with 
a different value of K. Maybe you should also change K^. Can 
you match your data and the calculated curve by choice of K? If 
so, what is the best value of K? How reliable do think the value 
of your best K is? Express K as some number, ± some enror 
estimate. If you can't match data and calculations, then report 
what, in your judgment, is the problem. If you suspect the data, 
then return to the laboratory and get additional data which are 
more reliable. 



volume 

of 
titrant 




16 



Chapter I 



Titration of a 
Monoprotic Acid 



CARBOXYLIC ACIDS: 

A rnMMnN class of compounds 

Some would say that the simplest of the carboxylic acids is 
CARBONIC ACID, a chemical compound that has never been prepared 
as a pure substance. Carbonic acid (HaCOs) appears to exist only 
in solution. A solution that contains carbonic acid can be pre- 
pared by bubbling carbon dioxide gas (CO2) through water. This 
formation of carbonic acid can be described by the chemical equation 

CO2 + H2O t H2COS (1) 
Equation 1 is a symbolic chemical model, written to describe some 
observations . 

Does carbonic acid really exist? Is it a real compound? 
Perhaps. But real or not, H2CO3 is a useful chemical species to 
use in describing some very real reactions of carbonates, bicarbon- 
ates, and carbon dioxide with acids, bases, and water. A common 
model of these reactions consists of the reactions of carbonic 
acid, bicarbonate ion, and carbonate ion, written as the pair of 
chemical equations 

H2CO3 # H* + HCO," (2) 
bicarbonate 

HCO,' H* + CO,'" (3) 
carbonate 

In many respects. Equations 2 and 3 have been found to describe 
faithfully the observed behavior of solutions. This chemical model 
has been found to be consistent with experimental observation. 

Carbonic acid has only a few atoms per molecule, and it may 
appear to be a simple species. However, description of a system 
containing carbonic acid may require Equations 1, 2, and 3, as well as 
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the equation for the dissociation of the solvent, 

H2O it H"" + OH" 

A solution containing carbonic acid can be described as a system 



of coupled and competing simultaneous chemical reactions. Such a 
Solution may be a rather complicated chemical system. We shall 
find that many aspects of carboxylic acid equilibria can be discovered 
in a study of a monoprotic acid (an acid that has only one proton 



other acids whose structures are given in Table 1. We shall think 
about the system at equilibrium, ignoring the rate and the manner 
by which the system achieves its equilibrium state. Then we shall 
look at additional features of carboxylic acid systems by examining 
some diprotic acids in which each proton dissociates from a separate 
site (proton binding sites are not unambiguously definable for the 
various carbonic acid species). We shall return at last to carbonic 
acid, but this return will be in the form of a literature-search 
and literature- interpretation project. 

A molecule of a carboxylic acid listed in Table 1 contains the 
structural unit 



called the CARBOXYL GROUP. The proton can dissociate, giving the 
carboxylate anion and the hydrogen ion: 



We shall now focus on compounds that contain just one of these 
carboxyl groups. 



that dissociates) such as formic acid, acetic acid, or one of the 





and H**" 



EMC 
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Table 1 

STRUCTURES OF SOME MONOPROTIC CARBOXYLIC ACIDS 



n 



'H 

H 



c£-c-cr 



H H 
( ( 



( ( ^ 
H-C-C-(X 

H H 



H H 
H-C-C-CT 



H-Cs. formic acid 

H 

H-C-C. acetic acid 

H 
H 



monochloroacetic acid 



trichloroacetic acid 



propionic acid 



lactic acid 



CHEMICAL MODcL FOR A MONOPROTIC CARBOXYLIC ACID 



Let us imagine a monoprotic acid that can exist in solution 
both as the associated acid HA and as the dissociated anion A . 
Such a situation can exist with any of the compounds given in Table 
1. Solutions can be prepared by mixing quantities of the acid, 
the sodium salt, sodium hydroxide, and water. Such solutions may 
have pH values that range from the acidic region to the basic region. 
We shall construct a model for this system that allows the prediction 
of pH values from statements of the amounts of acid, base, and 
solvent used for preparing the solution. In order to make this 
discussion specific, we shall focus on the titration of a carboxylic 
acid with sodium hydroxide. 

The species present in solution are (we shall assume) as follows: 



The chemical equilibria required to establish equilibrium in this 
system are (we shall asume) as follows: 



This listing of the chemical species assumed, together with this 
listing of the chemical equilibria assumed, constitutes the chmltal 
modzl {^on, thz 6y6tzm. 

Associated with each chemical reaction is an £.qulZlb^lum- constant 
equation which we shall always write in a standard form, with the 
molar concentration of a chemical species symbolized by enclosing the 
formula of the species with ( ), and with the equilibrium constant 
symbolized by a capital letter K with some identifying subscript. 
Thus we write, to accompany Reactions 5 and 6, the algebraic equations 



H2O, HA, A', Na*, H*, OH 



(4) 



H2O ? H* + OH 
KA ? H"^ + A' 



(5) 



(6) 
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= (H*)(OH") 
(H*)(A'*) 



(7) 



HA 



(8) 



(HA) 



The numerator in each equation is a product of concentrations of 
^zaction pKoduct^^ and the denominator in each case is a product 
of concentrations of n,<L(i(itcint6 ^ with the terms "product" and 
"reactant" defined in terms of the chemical reaction equation as 
it is written in the chemical model. When the solvent occurs in 
a chemical reaction, its concentration aoes not appear in the 
equilibrium- constant equation; instead, we write the number one 
in its place. These equilibrium constants may be expected to have 
constant numerical values under constant conditions, but they will 
in general be functions of temperature, pressure, solvent, and 
concentration. Hopefully, an equilibrium constant will not change 
its value very mucli with concentration within the concentration 
range which is encountered in the experiments under consideration. 

The requirement of electrical neutrality of the solution 
requires that the number of positive charges in the solution be 
equal to the number of negative charges in solution. Written in 
terms of molar concentrations, the zlzttKltat no^utKatlty equation 
is 



We shall also require that there be conservation of mass and 
non-transmutation of elements, and that chemical reactiors that are 
not included in the model do not occur. Such requirements can be 
stated algebraically in terms of a set of ton^z^vatlon zqaation6^ 
of which we shall require two here: 



(A') + (OH") = (H"") + (Na*) 



(9) 



[A] = (HA) + (A^) 



(10) 
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[Na] = (Na"") 



(11) 



Rquation 10 states that all the carboxylic acid orijiinally placed 
in solution is distributed between the species IIA and A . The 
quantity [A] represents the total of all A-containing species, in 
moles per liter. Equation 10 is an algebraic statement that the 
acid, when added to water, cannot disappear, but can be transformed 
into any of the species of the model if there is a chemical equation 
for the transformation, ""nee there is only one sodium-containing 
species in the model, Equa ^on 11 is particularly simple. 

It is helpful to distinguish between two sorts of concentrations: 
we shall speak of macfio people conczntKatloM ^ symbolized by letters 
enclosed by [ ]*s, and microscopic concentrations^ symbolized by 
letters enclosed by ( )*s. Macroscopic concentrations can generally 
be known by the experimenter from the manner in which the solutions 
were prepared; macroscopic concentrations are thus experimentally 
observable quantities. Microscopic concentrations are concentrations 
of individual chemical species, and often these concentrations can 
be found only by interpreting experimental data in terms of a model. 
One goal ot our manipulations of this model of an acid is to 
obtain ways to relate observable macroscopic concentrations and 
the microscopic concentrations of the species listed in the chemical 
model . 



Equations 7,8,9 10, and 11 comprise a set of simultaneous 
algebraic equations that describe in quantative terms the chemical 
nodel of a carboxylic acid. This set of equations is the mathmatlcal 
mtdtt ioK tht systtm. These five equations can be combined to 
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eliminate all the microscopic concentrations except (H*) , giving 
the cubic equation 

(H*)' + {[M] + Kjj^lCH*)^ ^^^^ 

This equation is written in terms of a rather special microscopic 

concentration, (H^) , which is used here as a ma^ttn. van.iablz. We 

shall use (H^) to describe the chemistry of the system, largely 

because of the relationship between (H^) and the experimentally 

measurable quantity, pH, We shall assume the equation 

pH = -ilogio(H^) (13) 

where pH is the reading on a properly-standardized pH meter, and 

(H*) is the molar concentration of the species called the hydrated 

hydrogen ion, or the hydronium ion. 

It appears, from a casual inspection of Equation 12, that 

there are two macroscopic variables: [A] and [M] . In an ordinary 

titration, however, [A] and [M] are not independent variables. It 

all becomes clearer if we change variables so as to make the actual 

experimental variable — the volume of sodium hydroxide titrant added 

to the acid solution, called V^^^^jj — explicit and prominent in the 

equation. If the two solutions used in a titration contain only the 

solutes HA and NaOH respectively, if we call the volumes of the two 

solutions that have been mixed V„. and V„ and if the volumes of 

HA NaOH 

the two solutions are additive (that is, if ^^^^^^ = Vjj^ + V^^^^jj) , 
then 

V M 

[A] = «^_HA 

V + V 
HA NaOH 

and 
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V M 

NaOH NaOH . c^ 

[M] = (15) 

V + V 

HA NaOH 

where V^^^^^^ is the volume of sodium hydroxide titrant added, 

^HA original volume of carboxylic acid solution, 

^NaOH molarity of the sodium hydroxide solution, and 

^HA molarity of the carboxylic acid solution. 

Introduction of Equations 13 and 14 into Equacion 12 permits the 
elimination of [A] and [M] , and gives an equation that relates 
the two experimental variables (H*^) and V^^^^j^. The equation is 

NaOH ' ]iA[dznomlnaton.j 

denominator = (H^)' . {K,^ - M^^^^} (H^) =^ . {KHA^NaOH " ^HH^) ^ 
EXERCISES 

1. Combine the appropriate algebraic equations and obtain 
Equation 12. 

2. Learn about the use and limitations of Equation 13. One 
ruthori itive source of information is R. G. Bates 9^V^t^n.mination 
0^ pH, New York: John Wiley § Sons, Inc., 1964. 

3. Perform the algebra that takes you from Equations 12, 14, 
and 15 to Equation 16. 




PREDICTIONS OF THE MODEL FOR A TITRATION EXPERIMENT 



In a titration experiment, quantities of the titrant are 
dispensed from a burette into a solution of acid, with the volume 
of added titrant being recorded after each addition. The pH of the 
resulting solution is measured, after each addition of titrant. 
There results a pair of two numbers, V^^^^^ and pH, which characterizes 
the solution at each point in the titration. More titrant is added, 
and a new set of numbers, V^^^^^ and pH, are recorded. When this 
process is continued throughout the titration, a set of ordered pairs 
^%aOH* obtained that can be plotted as a ^itfiation cunvz oi 

thz canboxytic acid. If required, we can use the operational defini- 
tion 

pH = -aogioCH'') 

to convert the ordered pairs fV^j^^jj, pH} into the ordered pairs 

f\aOH' f^^^H- Our task, in using our model to simulate this 

titration experiment, is to obtain a predicted titration curve in 

a form suitable for comparing with actual experimental data. 

Numerical substitution of a value of (H"*^) into Equation 16 

yields a value of V^^^^^, provided that numerical values for K^^^ and 

are already known, or are assumed. 

Methods for making reasonable estimates of the numerical values of 
these two equilibrium constants will be discussed in a later section. 
For purposes of understanding the following computer programs, just 
assume that the appropriate numbers have somehow been revealed. If 
you have your own data, and want to make a prediction of the titra- 
tion curve, you should Jump ahead to the later section , and find out 
how to make an informed guess. 

We shall program the IBM 1130 computer to perform this calcu- 
lation for us for many values of (H"^) . A program that will accomplish 

erJc '^5 



this task is FORTRAN Program 1, listed on page 28, with sample output 
given on pages 28 and 29. Let us examine this program, since it is 
typical of programs that can do the drudgery of repetitive calculations 
for us, leaving us free to think about interpretations and meaning 
of the numbers. 

Each typewritten line in this program is called a 6tatzfmnt. 
In each case, the statement was given to the computer in the form 
of a single punched card. The first four cards are conViot ciaKd^; 
they must be present and in the proper order for the program to 
function. (Actually, the * LIST SOURCE PROGRAM statement is not 
needed for the program to function; when it is included, the printer 
will type out the program as it appears on page 28. Such a print- 
out is helpful in spotting mistakes when you are running the program 
for the first time; later on, you may find that such a listing of 
the program takes too much time, and is not helpful, and you will 
then want to remove this control card.) The C on each of the next 
three cards indicates a commznt cafid. A comment card does not carry 
instructions to the computer for the program, but is included only 
to remind the user of special features of the program. 

A special feature of the FORTRAN language that is often helpful, 

but is annoying to us at this point, is the use of the letters I, J, 

K, L, M,andN (the letters from I to N, the first two letters of 

INteger) as first letters of variables used as integers. Integers 

are numbers without decimal points. Since we shall wish to use the 

variable names KHA , KW, NUM, MHA , and MMOH (because they remind us 

of K„^ , K , numerator, M„. , ana M„ ^„)» we tell the computer that 
HA' w 'HA' NaOH *^ 

we intend to violate these rules of FORTRAN grammar; we include the 
statement REAL'KHA, KW , NUM, MHA, MMOH. In FORTRAN, the word ''real'' 
is used to describe numbers with decimal points. 

A series of cards follows next to define various variables and 
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// JOB S71007tS5l 



// FOR 

» LIST SOURCE PROGRAM 
»IOCS(CARDt 1132 PRINTER) 

C TITRATION -CURVE 

C PREDICTIONS OF THE CHEMICAL MODEL 

C MONOPROTIC ACIDt ASSUMED VALUE OF KHA 

REAL KHAt KWt NUMt MHAt MMOH 

KHA = l,75E-5 

VHA = 25.00 

MHA = 0.1000 

MMOH = 0.1000 

KW = 1.008E-14 

H = 1. 

1 H = .5»H 

PH = -(ALOG(HJ ) /2.303 

NUM = H»»3 + KHA»(H»*2J - (KW + {KHA»MHA))»H - KW»KHA 

DEN = H»»3 + (KHA + MM0H)»(H»»2) + ( KHA«KMCH - KW)*H - KW»KHA 

VMOH = -VHA»NUM/DEN 

IF (VMOH) 2t 20t 20 

2 IF (PH - 7.00) It It 30 
20 WRITE (3tlO) VMOHt FH 
10 FORMAT (F10.3t FIG. 3) 

GO TO 1 
30 CALL EXIT 
END 

// XEQ 



PRINTER OUTPUT: 



0.194 


3.009 


0.739 


3.310 


1.607 


3.611 


3.100 


2.912 


5.552 


4.213 


9.100 


4.514 


13.349 


4.815 


17.406 


5.116 


20.524 


5.417 


22.542 


5.718 


23.707 


6.019 


24.336 


6.320 


24.663 


6.621 


O 
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24.830 


6. 


.922 


24.915 


7, 


.223 


24.957 


7« 


.524 


24,979 


7< 


.825 


24,99C 


til 


.126 


24.996 


8i 


.427 


25.000 


8< 


.728 


25.004 


9, 


.029 


25.010 


9, 


.330 


25.021 


y. 


.631 


25.043 


9. 


.932 


25.086 


10. 


.233 


25.173 


10. 


.534 


25.348 


10. 


.835 


25.702 


11. 


.136 


26.424 


11. 


.437 


27.933 


11. 


.738 


31.232 


12. 


.039 


39.239 


12 


i340 


64.818 


12 


• 640 


416.096 


12. 


.941 



constants, and to state the initial conditions of the titration. 

Since only capital letters (and no subscripts) arc available for 

communicating with the computer, we use such expressions as KHA and 

VHA to represent K^.^ and V^^^. The writer of this program thought 

that Mj^^Qj^ could be most clearly represented by MMOH rather than by 

MNAOH; you can use whatever combination of letters you wish, with 

a maximum of five letters per variable name. Large and small numbers 

are written in txpomntlat notation, which means that two of the 

FORTRAN statements translate as follows: 

KHA = 1.75E-5 => Kj^^ = 1 . 75 x 10" = 

KW = 1.008E-14 =o K = 1.008 x 10"^'* 

w 

Most of the remaining cards contain statements that constitute 
a computation loov, and we shall look at this loop in detail later* 
Statement 30 and the following card provide the necessary mechanism 
for ending the program. The last card in the deck is the instruction 
to execute the program. 

Unless there are instructions to the contrary, the computer 
proceeds sequentially through the statements. After one statement 
has been executed, the statement immediately following it is then 
executed. However, you as the programmer may change the sequence 
of operations by using the GO TO or IF statements: 

tIThe statement 

GO TO 1 

will be followed by execution of the statement numbered 1. Then 
the program follows sequentially after statement 1. 
^The statement 

IF (A) 1, 2, 3 
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provides for transfer of the program to one of three different 
statements, depending on the numerical value of the variable A. 
If A is negative, then the next statement executed will be the 
statement numbered 1. If A is zero, then the next statement that 
is executed will be statement 1. If A is positive, then the next 
statement executed will be statement 3. 

Note that only some of the statements are numbered. You may- 
number any statement with any number, remembering to never use the 
same number twice. 

You can add spaces within a statement, almost at your own 
convenience, in order to make the typewritten statement more 
readable to you. If there is doubt in your mind about the order 
in which the computer will perform the indicated arithmetic within 
a statement, then add some clarifying parentheses; parentheses have 
the same meaning in FORTRAN as in arithmetic and algebra. You must . 
always use a closing parenthesis with each opening parenthesis. . 
Other arithmetic symbols used are: 
+ dddJitJion 6jign 
6ubt/iactZon 6lQn 

* muttA.plA.aatA.on 6tgn (do not indicate multiplication 

with just parentheses) 
** means n^aihtd to thz powen. o^ 

= zqual6 6lgn 

The printer instructions to print ordered pairs of values of 

VMOH and PH can be written as 

WRITE (3,10) VMOH PH 
10 FORMAT (F10.3, FlO.3) 

where the 3 in the WRITE statement is the code number for the printer, 



and the 10 refers to FORMAT statement 10, which gives some details 
about how the printer is to print the output. The list of variable 
names in the WRITE statement indicates the variables whose numerical 
Values will be printed by the printer. The notation 10.3 means 
that a total of ten spaces has been reserved for typing the number, 
with three places following the decimal point. The F i^o^ ''floating 
point^^) means that the number will be printed as an ordinary number 
with a decimal point (that is, not in exponential notation, and not 
as an integer) . 

The heart of this particular FORTRAN program is the computation 
loop. A large value (larger than will be encountered in the actual 
titration) ^or (H^) is assumed; the concentration of is symbolized 
by the letter H. Then this value for (H^) is reduced to 501 of 
its previous value in statement 1, the pH value is calculated for 
future reference, and then (H^) is substituted into Equation 16. 
A value of titrant volume is calculated, and it is tested in two 
IF statements to see if it is negative, positive, or zero; and if 
it is negative, whether the titration is in the acidic or basic 
region. Only positive values of volume are acceptable as having any 
chemical meaning. Negative values are encountered at the beginning 
of the calculations, and here the computation loop cycles without 
any output being printed. When (H^) gets small enough to produce 
a positive volume, the cycle begins to include the WRITE statement 
and output results. Finally, when (H^) becomes smaller than its 
value in the titrant itself, the calculations pass outside the 
realm of chemical reality, the volume becomes negative, and the 
program is terminated. 
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// JOB 



S71007.S51 



// FOR 

» LIST SOURCE PROGRAM 
*IOCS{CARD» 1132 PRINTER) 



C 

c 
c 



TITRATION CURVE 

PREDICTIONS OF THE CHEMICAL MODEL 
MONOPROTIC ACID* ASSUMED VALUE OF KHA 



REAL KKA» KW« NUM» MHAt MMOH 

KHA » l,75E-5 

VHA • 25.00 

MHA « 0.1000 

MMOH = 0.1000 

KW « 1.008E-1^ 

H » 1. 

1 H « •8*H 

PH « -(ALOG(H) )/2.303 

NUM » H»»3 + KHA»(H»»2) - (KW + (KHA»MHA))»H - KW»KHA 

DEN « H»»3 + (KHA ♦ MM0H)»(H»»2) + (KHA»MMOH - KW)»H - KW»KHA 

VMOH = -VHA»NUM/DEN 

IF (VMOH) 2* 20* 20 

2 IF (PH - 7.00) 1» 1» 30 
20 WRITE (3»10) VMOH* PH 
10 FORMAT (F10.3t FiO.3) 

60 TO 1 
30 CALL EXIT 



END 



// XEQ 



PRINTER OUTPUT; 



0.038 
0.184 
0.339 
0.510 
0.703 
0.928 
1.194 
1.509 
1.886 
2.336 
2.871 
3.502 
4.240 



2.906 
3.003 
3.100 
3.197 
3.294 
3.391 
3.488 
3.585 
3.681 
3.778 
3.875 
3.972 
4.069 




33 



ERIC 



D ( 




A 
H 4 


. 1 A A 

1 loo 


1 

• 




I U3U 


0 
7 4 


"7 Q A 
I f 00 


i 




*T ( 


i 263 


C, D i 


1 vy .90 


0 ^ 
7 fl 


> ft ft 


7, 


l1 SI 

1 X 7 X 


/ft. i 


. "^A 0 

1 yj 




^ V 4 


. HAft 
J UHO 


7 fl 


> 7 f 7 


w ( 


> ^ •♦o 


( 


1 ^S7 






> HAO 


1 0 - 

X U fl 


»n7A 




f W ^ £m 


4 1 


9 J ^ J 




C 7 4 


n7R 


X U 4 


. 1 7*^ 




i984 


4i 


i650 




^ ^ 4 


^ no A 

1 U 7H 


1 U 4 


. 5 7n 




i372 


4 1 






2 ^ 4 


I X Ic 


iU 4 


Q AT 




,754 


4 1 


i ft44 




^ ? 4 


. 1 A7 
1 X** f 


1 0 - 

X V# 4 


1 HOH 








1 7** X 




2 ^ 4 


I X o4 


iU 4 


^^A 1 


16 




^ 4 


i 03ft 

1 w ^ O 


j 


^ 4 




1 V 4 


AKft 
1 0!>0 


17. 


1 w ^ £ 


^ 4 


1 X ^ ^ 




2^4 


O Q O 
I ^07 




» r DD 


1 ft 




^ 1 






Z 0 4 


^ A O 


iU4 


Q & 1 
I 0 51 


IQ 




7 4 


1 ^£ 7 




x5 4 


A & A 

1 45H 


xU 4 


1 74c 




1 ^ ^ 


^ 1 


1^2 S 




^ 0 4 


K A Q 
I DOO 


X 1 4 


. r\ A** 
1 UH5 






^ 1 


. S22 


• 


2 D 1 


"7 10 

1 71 3 


1 1 

1X4 


1 /- 0 
1 14^ 


? 1 1 


• ^ o ^ 


7 4 


A1 O 

I O X 7 




£ 0 4 


Q OA 
1 O7H 


1 1 4 


5^0 

1 A .37 


C £ i 


I ^ D 


7 4 


. 71 A 
1 r X O 




0 A 

4 


* lc.D 


1 i 4 


1 5 50 


^ C i 


1 7 O W 


7 4 


. ft 1 
1 O X ^ 




O A - 
4 


» A 1 0 
I HX ^ 


1 1 . 
1 1 4 


IH55 


£^ D 1 


t ^ O £ 


^ 1 


1 7 X U 




5 A - 
4 


« 777 
1 r r f 


1 1 - 
X X 4 


1 5 5U 






O 1 


k 007 




57 . 
^ r 4 


. 9 a5 
1 ^ H^ 


1 1 - 
X X 4 


.A 57 
1 0 A f 




iQ26 

1 7 ^ V 


A i 


k 1 04 




57 - 
^ r 4 


1 0 ^H 


1 1 - 

X X 4 


k72^ 

1 / AH 




i 1 34 


6 i 


> 20 1 

1 C W X 


• 


5 ft - 

^ O 4 


1 3 7H 


JL X 4 


. ft 5n 
1 0 aU 




i302 


6 1 


i298 




^ 7 4 


. S7S 


1 1 i 

X X 4 


.017 

» 7 X » 


24< 


»438 


6 1 


i394 


J 


D \J 4 


. ft S 5 
1 0 5 ^ 


1 5 - 

X ^ 4 


lU XH 


24i 


i649 

P ^ "T 7 


6 i 


k^Q 1 

> •♦7 * 






¥^ DO 


1 2- 

X £ 4 


kill 
1 X X X 


24, 


•637 


6i 


ft 6ft8 


1 


4 


k 7ftQ 

1 f 0 7 


1 0 - 

X C, 4 


. 50ft 
1 A v/O 


24. 


■ 709 

P V W 7 


6 i 


ft685 

t V u ^ 




"^7 - 
D 1 4 


. ft A A 
1 OOD 


1 0 - 

X C 4 




24i 


i767 


6 i 


k 7ft 2 




AO . 
4 


. 1 ft ft 
1 X Oo 


1 7 . 

X<^ 4 


1 HU< 


24< 


i813 


6 i 


i879 


• 


Aft - 
HO 4 


. SHA 


1 5 - 

X ^ 4 


. AQQ 
1 H 7 7 






A . 
O 1 


. OTA 
1 7 f O 




c o 

584 


i296 


124 


i59o 


5 A . 


^ Q Q r\ 


7 
r 1 


t\J f ^ 


j 


744 


i933 


12 4 






• 7 U*f 


7 , 
r 1 


I X r U 




1084 


lot 

»185 


124 


*T 0 
i789 




. Q ? ^ 


7 ^ 


1 cOO 




2034 


»028 


124 


1886 


£ H I 


1 7 ^ O 


7 , 






20504 


i702 


124 


i983 




t 7 ^ w 


7 , 


k^AO 

1 HO U 














1 7 Q V 


7 , 

f 4 


» SS7 


• 










5Zi . 
^ H 4 


1 7 O C 


7 

1 4 


AS^Zft. 












2^1 


1 7 f ^ 


7 ^ 


» 7S 1 
1 r ? X 












^ H 4 


1 7 O U 


7 ^ 


^ ft Aft 

1 OHO 












^ H 4 


1 7 O** 


7 ^ 


^ OA ^; 

1 7H^ 












2£i 1 
£ H 4 


1 7 O f 


ft < 














2£ii 


> 7 7 W 


ft « 

O 4 


k 1 "^ft 
1 X ^ o 












C.H i 


• 7 7 C 


ft . 
O 4 




• 












> 7 7 


ft < 
O 4 


» 9 
1 ^ ^ ^ 














> 7 7 O 


ft i 
9 4 


»£l2Q 

1 H 7 












24. 


P 7 7 f 


ft 1 
O 4 


i S26 

1 V b C7 












24. 


■ 999 

P 7 7 U 


ft 1 


ft623 












24. 


■ 999 

P 7 7 7 


ft i 


k 7? 0 












2 5 1 


.001 


ft « 
O 4 


kft 1 7 
1 O X f 












C 7 4 


\\J\J c 


ft - 


» Q 1 A 
1 7 X H 


• 










C, D 4 




o . 

7 1 


» Ul X 












^ ? 4 




O - 

7 1 


1 X U r 












254 


i007 


9fl 


P «• w ~ 












25i 


»009 


9i 


.301 












25< 


»012 


9i 


.398 












25< 


»015 


9i 


.495 












25< 


»019 


9i 


.592 












25i 


»024 


9i 


.689 











Perhaps you may feel that the output listed on pages 2 8 and 29 
inadequate to give a realistic prediction of a titration curve 
that is, in fact, a continuous function. Some of the output pairs 
are rather widely spaced. It is easy to get the output more-closely 
spaced; we simply make the decrements in (H^) smaller in Statement 
1* Such a revised program, together with the resultant output, is 
given as FORTRAN Program 2 on pages 33 and 34. 

A GRAPHICAL PREDICTION OF THE MODEL FOR A TITRATION EXPERIMENT 

The numerical output of the previous computer calculations 

is a prediction with greater precision than the corresponding 

experimental data. However, most people find it rather difficult 

to visualize the shape of the predicted function by looking at the 

long listing of ordered pairs of numbers. The computer program 

v/hich we shall now examine yields a graphical prediction. The 

principal readout from the computer is a graph of the titration 

curve, drawn and labelled by the plotter. The program, FORTRAN 

Program 3, is presented on page 37 and the plotter output is given 

on page 38. We shall examine some aspects of this program which 

pertain to the use of the plotter as the output device. 

PLOTTER The plotter preparation statement in the program is 

PREPARATION 

SUBROUTINE of the form 



where 



CALL PREP (XS, YS, XO , YD) 

XS is the scale chosen for the x-axis, written in inches 
per unit, using a number with a decimal point. 

YS is the scale chosen for the y-axis, written in inches 
per unit, using a number with a decimal point. 
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XO is the distance in inches from the left-hand edge of 
the graph to the origin. In'rliulo decimal point. 

YO is the distance in inches from the bottom of the graph 
to the origin. Include decimal point. 



tells the computer to write, using the plotter (the plotter h 5 the 

code number 7), the numerical value of K„^ according to the format 

HA ° 

given in tp. iiient 15. The fovmat information contained in state- 
ment 15 tells the computer to write the characters typed on irf 
card between the two ' marks. This statement is used to provide a 
label for the graph. Similar statements are used to label the 
X-axis and the y-axis. 

Following the CALL PREP statement, a WRITE (7, ) statement 
results in a label placed at the top of the graph. Following each 
of the CALL XAXIS and CALL YAXIS statements, a WRITE (7, ) state- 
ment results in a label placed along the corresponding axis. Other 
words and characters can be placed anywhere on the graph; see the 
IBM plotter instruction manual for detailed instructions. 

SUBROUTINE The y-axis statement in the program is of the form 



where 

X and Y are the coordinates of the origin of the axes. 
U is the distance between tick marks on the y-axis. 
N is the number of tick marks; thus N x U is the length of 
the y-axis. 

NLAB is the number of tick ryiarks that are to be labelled 
with che value of the scale at that point. 



WRITE 
STATEMEl.T 



The statement in the program 



WRITE (7,15) KHA 



YAXIS 

SUBROUTINE 
XAXIS 



CALL YAXIS (X,Y,U,N,NLAB) 



The corresponding statement is used for the x-axis. 
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// FOR 

♦ONE WORD INTEGERS 
• LIST SOURCE PROGRAM 

♦lOCSCCARD. il32 PRINTER* TYPEWRITER* PLOTTER) 
C TITRATION CURVE 

C PREDICTIONS OF THE CHEMICAL MODEL 

C MONOPftOTIC ACID* ASSUMED VALUE OF KHA 

REAL KHA* KW* NUM* MHA* MMOH 

KHA « 1.75E-5 

VHA - 25.00 

MHA = 0.1000 

MMOH = 0.1000 

KW « 1.008E-1^ 

CALL PREP (.2* .55* 1.* l.J 

WRITE (7*15) KHA 

15 FORMAT ('PREDICTION OF A CHEMICAL MODEL WITH KHA ««*E11.3) 

C CALCULATION AND PLOTTING OF THE TITRATION CURVE 

H « 1. 

1 H « .5*H 

PH « -(ALOG(H) )/2.303 

NUM « H»»3.+ KHA»(H»»2) - (KW + (KHA»MHA))»H - KW»KHA 

DEN « H»»3 + (KHA + MM0H)»(H»»2) + ( KHA«MM0H - KW)»H - KW»KHA 

VMOH » -VHA»NUM/DEN 

IF (VMOH) 1*2*2 

2 CALL FPLOT (-2*VM0H*PH) 
5 H = H - .1»H 

PH = -(ALOG(H) )/2.303 

DEN = H»»3 + (KHA + MM0H)»(H»»2) + (KHA»MMOH - KW)»H - KW»KHA 
NUM = H»»3 + KHA»(H»»2) - ( KW + (KHA»MHA))»H - KW»KHA 
VMOH = -VHA»!\UM/Dfc.N 
IF(50.-VMOH) 10*3*3 

3 CALL FPLOT (0*VMOH*PH) 
GO TO 5 

C DRAWING AND LABELLING THE AXES 

10 CALL YAXIS (0.* 0.* 1.* 15* 5) 

WRITE (7*17) 
17 FORMAT ( 'PH' ) 

CALL XAXIS (0.*0.*1.*50*10) 

WRITS (7*16) 

16 FORMAT ('MILLILITERS OF SODIUM HYDROXIDE') 
CALL FINPL 

CALL EXIT 
END 
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PREDICTION OF A MMICAL MODEL WITH KHA = 0-175f:-04 




EVALUATION OF EQUILIBRIUM CONSTANTS 
BY COMPARING PREDICTIONS Or A MODEL 
WITH THE DATA FROM A TITRATION EXPERIMENT 



If we proceed to interpret experimental data according to the 
chemical model described on page 21 we shall have reduced all the 
differences between various monoprotic acids simply to the differences 
in the values of the equilibrium constant K„^. Since we intend to 
describe the particular acid by specifying the numerical value of 
its equilibrium constant, it is important that we have a reliable 
method for evaluating that equilibrium constant. 

The shape of the titration curve depends on the value of the 
equilibrium constants K^^^ and K^. Thus it would seem that one way 
of obtaining a good value for the equilibrium constants would be to 
find the values for K^^^ and that yield a titration curve that 
gives the best fit to experimental data from an actual titration 
experiment. This is, in fact, what we shall do. 

FORTRAN Programs 1, 2, and 3 each give predictions of the 
titration curve for a monoprotic acid. In each case, however, it 
is necessary to supply a numerical value for both of the equilibrium 
constants before any calculations can be made. Thus we need a way 
to make an informed guess for such numbers. Then we need a convenient 
way to plot the predictions of the model on the same graph as the 
experimental data; or a convenient way to list the predictions of the 
model in the same table as the experimental data. We shall proceed 
first to develop two ways — one numerical and one graphical — to 
use the computer to facilitate a comparison of data and model for 
various values of the equilibrium constants. 
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FORTRAN Program 4 has been written to accept (as input) your 
experimental data from a titration, and to yield output as a tabular 
comparison of observed and calculated values of titrant volume for 
each of the experimental pll values. The input data are in the form 
of a set of ordered pairs {V^^^jj, pH}. Each pair of numbers is 
punched on a separate card, and the computer reads these cards one 

by one. From each card, the card-reader reads the values of V„ ^„ 

NaOH 

and pH. The computer converts each pH value to a value of (H^) , 
and calculates a value of V^^^^^^ which we shall call VMOH (CAL) . Then 
the computer prints out VMOH (CAL) and the experimental value of 
\aOH ^^^^^^ iEX?) originally read from the card. 

Statements 28 and 25 give instructions to the card reader (the 
number 2 is the code number for the card reader) to read the next 
card in the hopper, to interpret any number punched in the first 20 
spaces as V, and to interpret any number in the next 20 spaces as 
PH. The computer (because of FORMAT statement 25) expects i lat each 
number will have a decimal point; no matter how many decimal places 
are in the punched number, the number will be read with two decimal 
places and additional digits will be ignored. These data cards are 
placed immediately following the // XEQ card. 

Every program requires a mechanism for ending itself. Here we 
use the device o. placing a negative number in the first 20 spaces. 
The IF statement following statement 25 transfers the program to 
tne next statement for each data card that actually contains 
experimental data (all experimental volumes are zero or positive), 
but terminates the program when V is negative. Thus, by punching 
the last data card with a negative value of volume, we end the program 
after all the experimental data have been read and processed. 
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// FOR 

♦ LIST SOURCE PROGRAM 
»IOCS(CARO» 1132 PRINTER) 

C CONFRONTATION OF DATA AND MODEL. A NUMERICAL COMPARISON 

C TITRATION CURVE 

C MONOPROTIC ACID» ASSUMED VALUE OF KHA 

C DATA CARDS HAVE VMOH PUNCHED WITH IN THE FIRST 20 SPACES AND 

C PH PUNCHED WITHIN THE SECOND 20 SPACES 

C INCLUDE DECIMAL POINT IN EACH VALUt OF VMOH AND PH 

C LAST DATA CARD MUST HAVE A NEGATIVh VALUE OF VMOH 

REAL ICHA* KW» NUM» MHA» MMOH 

KHA = 1.75E-5 

VHA = 25.00 

MHA = 0.1000 

MMOH = 0.1000 

KW « 1.008E-14 

28 READ 12.25) V» PH 

25 FORMAT (F20.2. F20.2) 
IF (V) 30» 29. 29 

29 H = 10**(-PH) 

NUM = H**3 + i;hA*(H**2) - (KW + {KHA»MHA))*H - KW*KHA 
DEN = H**3 + (KHA ♦ MM0H)»(H**2) ♦ (KHA*MM0H - ICW)»H - ICW»KHA 
VMOH = -VHA»NUM/DEN 
WRITE (3»10) PH» VMOH. V 
10 FORMAT (IX.'PH » • »F6 .2 »4X » • VMOH (CAL) = ' »F7.2 »4X » • VMOH (EXP) 
•.F7.2) 
GO TO 28 

30 CALL EXIT 
END 

// XEV 



PRINTER OUTPUT: 



PH 


S 


3.42 


VMOH 


(CAL) 


1.00 


VMOH 


(EXP) 


= 1.00 


PH 


Jl 


4.11 


VMOH 


(CAL) 


= 4.57 


VMOH 


(EXP) 


= 4.47 


PH 


z 


4.51 


VMOH 


(CAL) 


9.02 


VMOH 


(EXP) 


= 8.89 


PH 


s 


4.58 


VMOH 


(CAL) 


9.97 


VMOH 


(EXP) 


= 10.26 


PH 




4.85 


VMOH 


(CAL) 


' 13.82 


VMOH 


(EXP) 


= 14.53 


PH 




4.96 


VMOH 


(CAL) 


15.36 


VMOH 


(EXP) 


= 15.96 


PH 


m 


5.32 


VMOH 


(CAL) 


' 19.62 


VMOH 


(EXP) 


« 19.42 


F 


<5 


5.54 


VMOH 


(CAL) 


' 21.46 


VMOH 


(EXP) 


= 21.19 



In actual practice, we would run this program several times, 
each time with a difrorcMit value of KHA. The ••best'* value of the 
equilibrium constant would bo the value that gave the best lit of 
experimental and calculated values of titrant volume. 
EXERCISES 

1. Certain regions of the simulated titration curve, calculated 

with programs such as.l, 2, 3, or 4, are markedly more sensitive to 

the particular value of the equilibrium constants than other regions. 

Determine, by running such a program with a range of equilibrium 

constant values, what experimental data will be most useful in 

determining the best value of K„,. 

HA 

2. The computer can be helpful in making a decision about which, 
of several sets of output volumes in program 4, represents the best 
fit of data and model. One criterion for best fit might be that, in 
the range where the calculated curve is most sensitive to changes in 
the value of the equilibrium constant, the differences between calcu- 
lated and experimental values of volumes should be as small as possibl 
Thus we might find a program that made such a subtraction, and 
printed out the difference, to be helpful. Even more. helpful would 

be a program that summed the absolute values (or the squares) of 
these differences, and printed this sum at the end of the calculations 
Write such a program, and execute it with your data, and with several 
values of an equilibrium constant. 

A GRAPHICAL FORTRAN Program 5 takes the experimental data from a 
COMPARISON OF 

DATA AND THE titration experiment, punched on cards in the same 
PREDICTIONS 

OF A MODEL format as required by program 4, and plots these data 
on the same graph as a calculated plot of the predic- 
tions of the chemical model. When such comparisons are made for 
O 42 
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several different values of an equilibrium constant, inspection of 
the computer -drawn graphs may be sufficient to permit an informed 
estimate of the best equilibrium constant value. 

The graphical output from this program gives the prediction of 
the chemical model as a smooth curve, and the experimental data as 
points plotted on the same graph. Such output, using the same data 
cards as used for demonstrating program 4, is presented on page 46. 
The model and the data are consistent if the points fit the curve 
within the limits of the experimental uncertainty of the data. 

• What l6 an dxptanatlon zach thz iottovolng dz\jJiatlon^ 
bztu)zzn pZottzd points and tattutatzd cuKvz? 

II The data points lie above the curve for volumes between 
5 and 20 ml for a titration curve that has a sharp inflection point 
at about 25 ml. 

II The data points lie below the curve for volumes between 
5 and 20 ml. 

H The data points lie to the right of the curve forpH values 
from 7 to 10. 

II The data points lie below the curve for volumes greater 
than 30 ml. 

H The data points lie randomly above and below the curve 
for volumes between 5 and 20 ml. 

SOME NOTES ON The smooth solid curves drawn by the plotter are 
THE COMPUTER 

PLOTTING PROGRAM in fact constructed by instructing the plotter 

to draw a series of very short connected straight 
line segments. The procedure is to bring the pen to one end of the 
desired line, with the pen up (lifted off the paper) so that the 
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// FOK 

»ONE WORD iNTEGtKS 

» LIST SOUKCE PROGRAM 

»IOCS(CARDt 1132 PKINTERt TYPEWRITEKt PLOTTER) 
C TITRATION CURVE 



C 

c 



COMPARISON OF EXPERIMENTAL DATA WITH THE PREDICTIONS OF A 
CHEMICAL MODEL 



C 



MCNOPROTIC ACIDt ASSUMED VALUE OF KHA 



C 

c 
c 
c 



DATA CARDS HAVE VMOH PUNCHED WiTHiN THE FIRST 20 SPACES 
DATA CARDS HAVE PH PUNCHED WITHIN THE SECOND 20 SPACES 
INCLUDE DECIMAL POINT IN EACH VALUE OF VMOH AND PH 
LAST DATA CARD MUST HAVE A NEGATIVE VALUE OF VMOH 



REAL KHAt KWt NUMt MHAt MMOH 

KHA = l,75E-5 

VHA = 25.00 

MHA = 0*1000 

MMOH 3 0*1000 

KW = 1*0C8E-14 

CALL PREP (*2t •55t l.t 1*) 

WRITE (7»15) KHA 

15 FORMAT ('CONFRONTATION OF CHEMICAL MODEL WITH Tn;?ATlON DATAt KHA 
X =«tEll«3) 

C CALCULATION AND PLOTTING OF THE TITRATION CURVE 



1 H = .5»H 

PH = -(ALOG(H) )/2*303 

NUM = H»*3 + KHA*(H**2) - ( Krt + (KHA*MHA))*H - KW*KHA 

DEN = H»»3 + (KHA + MM0H)»(H«»2) + (KHA»MMOH - KW)»H - KW»KHA 

VMOH = -VHA»NUM/DEN 

IF (VMOH) It2t2 

2 CALL FPLOT (-2tVM0HtPH) 
5 H = H - •1»H 

PH = -(Al.OG(H) )/2*303 

DEN = H:r»3 + (KHA + MM0H)»(H*»2) + (KhA*MMOH - KW)»H - KW*KHA 
NUM = H»»3 + KHA»(H»»2) - ( KW + (KHA»MHA))»H - KW»KHA 
VMOH = -VHA»NUM/DEN 
IF(50.-VM0H) 10t3f3 

3 CALL FPLOT (OtVMOHtPH) 
GO TO 5 

C PLOTTING DATA POINTS FROM THE TITRATION 

10 CONTINUE 

28 READ (2t25) VtP 

25 FORMAT (F20*2tF20*2) 
IF (V) 30t29t29 

29 CALL FPLOT (ItVtP) 
O CALL FPLOT (-2»VtP) ' L^O 44 



H 



= 1* 
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CALL POINT (1) 
GO TO 28 

C DRAWING AND LABELLING THE AXES 

30 CALL YAXI5 (0«t 0«* !*• 15* 3) 

WRITE C7tl7) 
17 FORMAT ( 'PH* ) 

CALL XAXI5 (0«tO«tI« t50tlO) 

WRITE (7tl6) 
16 FORMAT ('MILLILITERS OF SODIUM HYDROXIDE') 

CALL FINPL 

CALL EXIT 

END 



// XEQ 



ERIC 



64 



45 



COraTATION OF [MM. KEEL WITH TITRATION DATA, M - 




pen is not writing* The pen is then lowered. When the pen is then 
moved in the down position to the other end of the desired lino, the 
straight line is drawn. An instruction to move to another position 
yields a continuation of the line if the pen is left in the down 
position. The necessary information to control these pen movements 
is included in the computer program in the form of FPLOT statements 
of the form 

CALL FPLOT CI ,X,Y) 
where I is an integer that controls the up-down pen position as 
follows : 

= 0 => no change 

= positive ^ control pen before movement 
= negative ^ control pen after movement 
= odd =«> raise pen away from paper 

= even lower pen to writing position 

The variables X and Y are the coordinates of the graph position to 
which the pen is moved. 

If you want to place a printed label on the graph, use the 
FPLOT statement to bring the pen to the position where you wish the 
lower left corner of the first character to appear. The pen should 
be in the up position. Then use the statement 
WRITE (7, J) 

where J is the number of a FORMAT statement of the form 

J FORMAV ») 
The desired letters and numerals for the label go between the ' ' 
marks. 

If you want to plot a point on the graph, you may use the 
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statement 

CALL POINT (1) 

This statement causes the plotter to draw an x mark at the current 
pen position. It expects to find the pen down and it leaves the pen 
down when finished. Other numbers may be used in this statement, 
with other marks resulting; details are given in the IBM plotter 
manual . 

At the end of a graph, when all plotting and graphing is finished, 
it is convenient to use the statement 
CALL FINPL 

This subroutine plots a cross at the lower left hand corner on top 
of a plus sign put there by PREP. This is a check on machine errors 
during the plotting. If the two symbols do not land exactly on top 
of each other, the graph should not be trusted. FINPL then raises 
the pen and moves it to an appropriate position for beginning a new 
graph. 

USE OF THE MACROSCOPIC HALF-EQUIVALENCE METHOD 
FOR MAKING AN INFORMED GUESS OF THE VALUE OF THE 
EQUILIBRIUM CONSTANT FOR A MONOPROTIC ACID 

In order to confront the data from a titration experiment with 

the predictions of a particular chemical model, it is necessary to 

assume some specific numerical value for the equilibrium constant 

Kjj^. Only then can a theoretical titration curve be calculated for 

direct comparison with the experimental curve. If the agreement is 

not very good, then a different value for K„. can be tried. This 

trial -and-error approach — the method of successive approximations — 

will usually yield a •-best value'' of after three or four such 

HA 

48 



ERIC 67 



calculations if the first trial value of K^^^ is reasonably close to 
the value that finally gives the best fit between chemical model and 
experimental data. With a computer being used for these calculations, 
almost all the effort goes into the first calculation. Succeeding 
calculations, with different trial values of Kj^^, require only that 
you change a single card in the program deck, and each calculation 
takes only about five minutes of IBM 1130 computer time. We have 
already seen how to perform the calculations and get a readout, on 
the plotter, of the simulated titration curve plotted together with 
experimental data points obtained from an actual laboratory titration. 
What we need now is a simple way to make a good guess for a value 
of Kjj^ to begin the successive approximations process. 

During the titration of a monoprotic acid, the pH value of the 
solution will be equal numerically to the pK^^ of the acid when the 
concentrations of species HA and A" are equal. The criterion of 
equal concentrations of these two species is the m^ccAo^ copied ddijiYuLt^on 
0^ thd hali-zqiilvatzndz point oi tkz tltn.atJion. A ma(iH.o6coplc 
ha^^-equlvaZzncz point is often defined as the point in the titration 
at which a volume of base has been added which is equal to half that 
volume required for reaching the equivalence point. 

Microscopic and macroscopic half -equivalence points are not 

necessarily identical. If they are the same, within the limits of 

experimental errors of measurement of volume and pH, then the pH 

value at the macroscopic half-equivalence point provides the datum 

for a convenient determination of pK,.^ and therefore of K„,. In 

^ hA HA 

any event, the determination of the pH at the macroscopic half- 
equivalence point provides a convenient and straight-forward way to 
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get an approximate value of K^^^ to serve as a starting point for 
the successive approximations calculat ions. The K^^ so ohtainoU 
will be an excellent approximation for cases in wliich the actual 
value of the dissociation constant lies between 10*** a;id 10"^^ — 
where the actual ^alue of pK^^^ lies between 4 and 10, Even when 
^HA order of 10"^, the value given by the half -equivalence 

method is a useful starting point for a series of successive approx- 
imations . 

FORTRAN Program 6 enables us to evaluate the half -equivalence 
method by simulating the titrations of acids with Kj^^ \^alues ranging 
from 10'^ to 10'^^. The resulting computer-produced plots are 
presented on pages 53 to 58. For each of these simulations, a 
straight line is drawn with pH = pK„., and another with V,,^„ = 

tiA MOH 

^^^^MOH, equivalence* intersection of these two lines lies on 

the calculated titration curve, the macroscopic and microscopic 
half-equivalence points coincide. All these calculations are for 
titrations in which both acid and titrant have concentrations of 0,1 
moles per liter. The half -equivalence method becomes less reliable 
at lower concentrations,, but the loss of accuracy at pK^^^ = 4 is not 
ver^' dramatic at concentrations of 0,01 molar, for instance, 

SOME NOTES ON THE The titration curve for an acid with a large 
COMPUTER PROGRAM 

value of the dissociation constant — that is, for 
a relatively strong acid ~ begins with a slope close to zero. An 
aqueous solution of a strorg acid is a good buffer solution. This 
fact introduces a strange feature into these computer plots that is 
entirely a quirk of our programming: the calculated titration curve 
doesn't necessarily begin at V = 0! The problem is that in the 
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// FOR 

♦ONE WORD INTEGERS 

♦ LIST SOURCE PROGRAM 

»IOCS(CARD» 1132 PRINTER* TYPEWRITER! PLOTTER) 

C TEST OF THE VALIDITY OF THE HALF-EQUIVALENCE METHOD 

C FOR EVALUATING THE EQUILIBRIUM CONSTANT 

C MONOPROTIC ACID 

REAL ICHAt KW» NUMt MHAt MMOH 
KHA = l.E-^ 
VHA = 25.00 
MHA = 0.1000 
MMOH = 0.1000 
KW = l.OOSE-l^ 
CALL PREP (.2t .55* lot 1.) 
WRITE (7*15) KHA 
15 FORMAT ('TEST JF THE HALF-EQUIVALENCE METHOD* KHA « 'tEg.Z) 

C MARKING OF THE HALF-EQUIVALENCE VOLUME* AND PH = PKHA 

VMOH = .5*VHA*MHA/MM0H 
PH = -(AL0G(KHA5 )/2.303 
VMOH 1 = VMOH - 3. 
VMOH 2 = VMOH +3. 
PHI = PH - 1. 
PH2 = PH + 1. 
CALL FPLOTdtVMOHtPHl) 
CALL FPL0T(-2»VM0H»PH1) 
CALL FPL0T(-2»VM0H»PH2) 
CALL FPL0T(1»VM0H1»PH2) 
(.ALL FPL0T(-2»Vy0Hl»Ph) 
CALL FPLOT (-1 »VMCH2»PH) 

C CALCULATION AND PLOTTING OF THE TITRATION CURVE 

H = 1. 

1 H = .5*H 

PH = -(ALOG(H) )/2.303 

NUM = H»*3 + KHA*(H*»2) - ( KW (KHA»MHA))»H - KW*KHA 

DEN = H»»3 + (KHA + MMOH ) ♦ ( H»»i: ' + ( KHA«MMOH - KW)*H - KW»KHA 

VMOH = -VHA*NUM/DE.\ 

IF (VMOH) 1»2»2 

2 CALL FPLOT (-ZtVMOHtPH) 
5 H = H - .1»H 

PH = -(ALOG(H) )/2.303 

DEN = H*»3 + (KHA + MM0H)»(H»»2) + (KHa'^.MMOH - KW)*H - <W*KhA 
NUM = H**3 + KHA*(H**2) - ( KW + (KHA»MHA))*H - KW*KHA 
VMOH ' -VHA»NUM/DEN 
IF(50.-VM0H) 10»3»3 

3 CALL FPLOT (0»VM0H»PH) 
GO TO 5 

C DRAWING AND LABELLING THE AXES 7^ 51 



PAGE 2 



10 CALL YAXIS ( 0* tO* t 1* f 15 f 5 ) 

WRITE (7tl7) 
17 FORMAT { 'PH' ) 

CALL XAXIS (0««0**1* *50*10) 

WRITE {7tl6) 
16 FORMAT {'MILLILITERS OF SODIUM HYDROXIDE') 

CALL FINPL 

CALL EXIT 

END 



// XEQ 
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searching process at the beg inning of the calculat ions , when the 
program is looping until V^^^^ becomes positive, the pH increment is 
too large. It is possible for an overshoot of almost 0.5 pll units 
to occur. This is of no consequence for a weak acid, where the 
titration curve begins with a rather steep slope. The simple remedy 
for this problem (if you consider it to be a problem) is to reduce 
the pH increment to some smaller value by modifying one statement in 
the program. 

The section of the program for drawing the two intersecting lines 
contains some strange instructions ft. the form of the FPLOT state- 
ments. The procedure is to bring the pen to one end of the desired 
line, with the pen up so that the pen is not writing. The pen is 
then lowered, and moved while in the down position to the other end 
of the desired line. This process yields a straight-line segment. 
The necessary information for decoding these FPLOT instructions is 
given on page 47. 

A METHOD OF ESTIMATING THE VALUE OF K^^ 
FROM AN INSPECTION OF THE TITRATION CURVE 

Confrontation of data and model also requires that you assume 

some specific numerical value for the equilibrium constant K^. Then 

a series of successive approximations can be used to find the value 

that gives the best fit between calculations and titration data. 

A straightforward method of estimating the value of K involves 

w 

use of Equation 7: 

= (H^)(OH") 

A place on the titration curve is needed where you can know the 
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values of both (H*) and (OH'). 

If you continue the titration well past the equivalence 
point, you will find that the pH approaches a limiting value. 
This limiting value is the pH of the titrant. Eventually, if you 
were to add an extremely large volume of titrant, the resulting 
solution would be essentially pure titrant. If you can make a 
reliable estimate of the limiting pH value, then you can use 
Equation 7 to calculate K . This pH value is converted into a 
concentration, giving the value of (H"*") that results when (OH') 
is equal to M^^^^. 

Rather than to rely on guessing the value of this limiting 
pH, you may find that it is more convenient to measure the pH 
of the pure titrant solution. 
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WHERE TO FIND EQUILIBRIUM-CONSTANT VALUES IN THE LIBRARY 



Much quantitative information about chemical reactions has 
been reported in research articles in the form of numerical values 
of equilibrium constants* Such numbers for acid-base reactions 
have been collected in several useful compilations that are listed 
and described in this bibliography. Obviously, research results 
reported after a set of tables has been compiled will not appear 
in those tables. For recent research, consult Chemical Abstracts 
or Chemical Titles for references to the original research papers. 

GENERAL TABULATIONS 

L. G. Sillen and A. E. Martell, Stability Constants of 
Metal-Ion Complexes , 2nd ed., London: The Chemical Society, 
196^ • Supplement, 1971- 

An excellent and reliable source of equilibrium conr.tHiil. 
values. For each entry, there is information about the solvcMit. ^ 
ionic strength, temperature, and experimental method, togothcr* wil.li 
a roferenco to the original research report. This is a virtuiUly 
complete r.urvey of the chemical literature for acid-base reuctloiir. , 
metaJ-ligand reactions , inorganic oxidation-reduction react i otir. , 
and precipitation reactions. 



ACID-BASE REACTIONS 

D. D- Perrin, Dissociation Constants of Inorganic Acids 
and Bas es i n Aqueous Solution , London: Butt erworthr, , 1969* 

D. D. Perrin, Dissociation Constants of Organi c Bases in 
Aqueous Solution , London: But terworths , 19o5 

G. Kortflm, W. Vogel, and K. Andrussow, D isG^ciat j on Con;; t:tn tr; 
of Organic Acids in Aqueous Solution ^ London: Butierwor Ihr. , 
1961. 

Each of these volumes was published under the uur.pi***:: ni' 
the International Union of Pure and Applied Chemistry. Compl«-l.o 
information regarding solvent, temperature, experimental moth*Ml, 
data treatment, and a reference to the primary research publ j cai, i chi 
is given for every entry. 
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USE OF A VISUAL INDICATOR TO DETECT 
THE EQUIVALENCE POINT IN A TITRATION 

Thus far in this case study, we have been considering the 
use of titration data to learn about the chemical nature of a 
carboxylic acid in water solution. The most common use of the 
burette is in the titration of solutions to determine the quantities 
of some material in those solutions • Such titrations are among the 
most convenient methods of quantitative. volatmtn,ic anatyii6^ and 
depend for their reliability upon methods of detecting the titration 
equivalence point. For determining the amount of a carboxylic acid 
in a water solution, the color change of an acid-base indicator is 
especially convenient* We shall use a computer simulation to 
investigate the utility and the limitations of visual equivalence- 
point detection in the titration of a monoprotic acid. 

Let us examine a solution prepared by mixing aqueous solutions 

of the weak acid HA, the strong base NaOH, and the indicator HI. 

Visual detection of the equivalence point depends on the species HI 

and I having markedly different colors , with at least one of those 

species being highly colored so that only trace quantities of the 

indicator are needed in solution. 

Often it is necessary to dissolve the indicator in an alcohol-wat^r 
mixture in order to get enough indicator into a drop of indicator 
solution. However, in describing the chemistry of a titration, it 
is customary to ignore the drop of alcohol that is introduced into 
the titration mixture when the indicator solution is added* 

We shall examine the conditions under which an abrupt color change 

would be expected to coincide with the equivalence point, thus 

signalling the point to halt the titration and read the burette. 
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Consider the special case of a titration in which the indicator 
concentration is so low that throughout the titration (but in 
particular in the vicinity of the equivalence point) 

(!') « (A") (15) 

If 

the inequality 15 will be satisfied if 

[I] « [A] (17) 
It is often quite easy to arrange experimental conditions so that 
conditions 15, 16, and 17 are satisfied. 

The relevant equilibria and associated equilibrium constant 
equations that we need for constructing a model for this system are 
H2O # H"" + OH' K = (H"") (OH') (18) 



(H )(A") 



HA i» H + A K„ = (19) 

(HA) 



(H )(I ) 

HI ii H + I K„ = (20) 

(HI) 



The algebraic equations for conservation of mass are 

[A] = (HA) + (A') (21) 

[Na] = (Na*) (22) 

[I] = (HI) + (D (23) 

The electroneutrality relationship is 

(H*) + (Na*) = (A") + (OH') + (l') (24) 

The volume variables that tie our description direcMy to burette 

readings in a titration can be introduced into the model via the 

three equations 
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[A] = 



Vha^ha 



V + V + V 
"^HI ^HA ^NaOH 



(25") 



[Na] = 



V M 
NaOH NaOH 



V + V + V 
HI "ha "NaOH 



(26) 



[I] = 



^HI^^HI 



V + V + V 
"^HI ^HA ^NaOH 



(27) 



where V^^^ is the volume of acid solution (of molarity M„^) that 



HA' 



was originally placed in the titration vessel, V„^ is the volume 

of indicator solution (of molarity M^^^) added, and V^^^^^^ is the 

volume of titrant (of molarity Mj^^^^^) which has been delivere^^ from 

the burette at a given point in the titration. Only V„ ^„ is a 

' NaOH 

variable during a particular titration. 

With approximation 15 introduced. Equation 24 becomes 

(Jfj + (Na"") = (A") + (OH") (28) 
Equation^- 18, 19, 21, 22, 25, 26, 27, and 28 can be combined to 
give 



V 



w HA 



- V 



HI 



- K (H*) \ . 



w 



V HA 



(H'^V + {K^., + M„ ^„}(H'')2 + {K„.M„ ^„ - K }(H''') - K K, 
^ ^ hA NaOH^ ^ ^ ' ''■>. NaOH w I 



NaOH 

NaOH ^ w 'HA 

Equation 29 contains no reference to either HI or I as chemical 
species; the volume of indicator solution appears, but this solution 
functions in Eouation 29 only to increase the total volume of the 
titrated solution, thus slightly diluting the solution. The re- 
quirement tha*- [I] be small has resulted ip separation of variables. 
The indicator plays no role In e.6 tab Hi king the value of (H*) for 



(29) 
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any value of added titrant. The distribution of total I between 
the species HI and l' depends solely on the value of (H**") . Thus we 
see that this concentration variable (H**") functions as a ma^taK 
vaKiabtz^ providing the chemical coupling between equilibria 18, 
19, and 20, and providing the algebraic coupling between Equation 
20 and the rest of the mathematical model. 

We shall now proceed to make a quantitative prediction for the 
color change of the indicator in a titration. Quantitative measure- 
ment of color changes in solution is best made in terms of the 
optical ab60Kbo,ncz, A, of the solution at some appropriate wavelength, 
ADsorption of visible light by either of the two indicator species 
occurs when a photon of light is absorbed by electrons in the anion 
or the molecule. Association or dissociation of the proton changes 
the permitted wavelengths of light that can be absorbed by the-; 
indicator species. The rw;ult is a change in the color of the^ 
solution containing the i^icator. We are interoste^ in hou) much 
.light is* absorbeV by V\ie*; lA? utiQn\ This chang^\n tl^ amount of 
iignt of each wavelength ^^^bso^bed by the solution as ille titration 
proceeds is directly relatW \o the variation in concentrax:ions. of 
the chemical species in solution. 

A beam of monochromatic light will be reduced in intensity 
while passing through a solution that contains chemical species 
which absorb photons of that wavelength. \he decrease in intensity 
depends on the number of molecules of the absorbing species in the 
path of the light beam. If the beam has an intensity lo before 
entering the solution, emerging from the solution with an intensity 
I, then the Boug'3r-Beer Law states that, for a solution containing 



8-1 



only a single absorbing species X, 

I = lolC^^x^'^^ , (30) 

where I is the optical pathlcngth, the distance travelled through 
the solution by the light beam; (X) is the molar concentration of 
chemical species X; ar;d is a proportionality constant, called the 
molar absorptivity of species X. 

If there are two different absorbing species in solution, then 
Equation 30 is modified to become 



I = I( 



10'^^Xif^^))(lO-^^X.f^^)) (31) 



Species X and X are here assumed to absorb light independently of 
one another. Equation 31 can be written in the equivalent form 

I = lolO-^^'^X,^^^) " ^Xa^^^^^ (32) . 

To get from Equation 31 to Equation 32, use was made of the relation- 
ship 

{io^}{io^} = 10^^*^^ •"■ " > . i 

Quantitative measurements of light absorption in ^plution| |?re 
conveniently made in the laboratory with an optical absorption ^ 
spectrophotometer, and the measurements are usually repoiibed in -^rms 
of either t ransmittance , T, or absorbance, A. Transmitt ance is 
defined by 

I 

T = — * 

; 

Absorbance is defined by 

lo] 

7 (33) 



A = ilogio 



There are both advantages a;id limitations to the use of either T or 
A* Transmitt ance values are contained within the interval from zero 
to unity, whereas absorbance values can range from zero to infinity. 

r 
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Transmittance is related to concentration in an awkward way, but 
absorbance is a convenient linear function of concentration. 

A derivation and discussion of the Bouger-Beer law is given 
in M. G. Mellon, Analyt ical Absorption Spectroscopy, New York: John 
Wiley & Sons 5 Inc., 1950, pp. 90-101. This book remains an excellent 
source of information about spectrophotometry, although continuing 
rapid advances in instrumentation make it necessary to supplement 
such a book with the latest information from instrument manufacturers. 
Alternative derivations of the Bouger-Beer Law are given by H. A. 
Liebhafskey and H. G. Pfeiffer, J. Chem. Educ. , 30, 1*50 (1953) , and 
by J. H. Goldstein and R. A. Day, J. Chem. Educ. , 31, 1*17 (195^*). 

We shall now investigate how this linear dependence of A on 

concentration comes about. Equations 30 and 33 are combined, 

taking advantage of two fundamental properties of logarithms: 

AogiolO^ = X 



The result is 



10 



10^ 



\ 

for a^solut 



ion V 



A = le^iX) 



ontaining the single fbi?Qrbing Species X. 



\ 7 



(34) 
T% 



corresponding equation relating absorbance to concentrations for a 
» ». 

solution witlj. tw* absorbing chemical species is obtained by com- 

t 'Kit 



(35) 



bi)tiing Equatilo 3,1 and 33. . The result is 

All important chemical research tool is the optical aKsorption * 

{ 5 

spectrophotometer, an instrument used to measure the amount 9f s 
monochromatic light absorbed by a sample. White light combsjfrom J: 
an incandescent tungsten-filament electric light bulb for .-jieasure- 
ments in the visible portion of the spectrum, and this light is 
passed through a prism or grating monochromator . Monochromatic 
light comes from the monochromator as a narrow band of wavelengths 
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distributed around a nominal value shown on the monochromator 
dial. Numerical values of absorbance are obtained by determining, 
with a photosensitive detector, the intensity of a light beam 
that passes through a transparent cell containing the solution. 

If, at the wavelength chosen, the only significant light 

absorption is due to the two indicator species HI and l"^, then 

the absorbance of the solution is given by 

A = AejjjCHI) + Aej-CD (36) 

where e^^ and e^- are the molar absorptivities of species HI and 
I . We can conveniently measure A in the laboratory, but there 
is no direct method of measuring (HI) or (I), so we seek now 
to eliminate both (HI) and (I) from Equation 36* We seek an 
equation relating only experimentally-measurable variables. 
Combination of Equations 20 and 23 gives both 



[I] = (HI) 



1 + 



HI 



(H^) 



(37) 



[I] = (D 



(H") 



HI 



(38) 



Then introduction of both Equations 3i 
36 gives 



and 38 into Equation 



A - 



1 + 



'^HI 
(H*) 



HI 



(39) 



Rearranging, we get 
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A = 



(H*) + Kjjj (H*) + K 



HI 



(H") + K 



HI 



(40) 
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Introduction of volume variables by means of Equations 25, 26, and 
27 transforms Equation 40 into 



A = 



Vhi^hi 



V + V + V 
*HI ''ha ^NaOH 



HI 



^HI^" ^^'^I-'^HI 



(41) 



If, at some particular point in the titration, paired values of 
(H*) and V^^^^jj are known, then Equation 41 can be used to 
calcu'.ate the value of A at that point. Our computer calculation 
of A versus Vj^^^jj will use this procedure. Equations 29 and 
41 will be used together for the calculations. 

Computer calculations could give us insight into several 
questions a.jout the use and usefulness of indicators in the 
titration of an acid with a base: -y-'^Ss. 

• feoto macfe do tka vatau the. two ab^o^ptlvltlu have 
to dt^^Q./i In ofiddfi that a ^hoAp ab^o^bance. change can 
occu^? > 

• What afic the ^utfilctlonh on the value o^ K^^ con6l6tent 
i^lth an abrupt abco^bance change? 

• Voe6 an abrupt ab6 0Kbance change neceihafilty coincide 
u)lth the equivalence point? 

• l6 It po6Alble to have an abrupt change In ab6CKbance 
without a coKKe^pondlngly abrupt change In (H*)? 

• Can an ab^o^bance change occu^ be^o^e the equivalence 
point? a^tcK the equivalence pclnt? 
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FORTRAN FORTRAN Program 7 begins with a trial value of (H ) , 
COMPUTER 

PROGRAM calculates a value of Vj^^^oH* ^^^^8 Equation 29. 

If the calculated value of V^^^^j^ is negative, a 
smaller value of (H*) is chosen, and the calculation is made 
again. This procedure continues until a positive volume is 
obtained, and then that value of (H*) is also used to calculate 
the value of A by using Equation 41. A graph of A versus Vj^^^j^ 
is drawn from the points obtained from the series of further 
calculations, each calculation made from a smaller value of (H^) . 
Finally (H*) becomes so small that the calculated volume'"bf 
titrant again becomes negative, and then this phase of the 
program is completed. Then on the same coordinates, a plot of 
pH versus V„ ^„ is drawn with numbers obtained from a new set 
of calculations. The statement following statement 22 adjusts 
the scale axes so that pH, ranging from 4 to 13, can be plotted 
on the same axes as absorbance, ranging from 0 to 1. The 
quantities B, C, D, E, and F were introduced for convenience 
in punching cards, on the asiumptiol} ♦ttj^ t^e chance <^ eijrors 
in punching was thereby reduce|^. | 



conclusions: some The results of spme; representative \ 

GENERALIZATIONS ABOUT ^ \ 

THE PRACTICAL USE calculations are presented on the 

OF AN INDICATOR 

. lowing pages. It is clear that an 
abrupt change in color of the solution (an abrupt increase io 
the absorbanct for the cases presented here) can be achieved i^ 
K„- has a numerical value equal to a value of (H*) within a 
nearly-vertical portion of the pH-volume cu^ve near the equivalence 
point. The sharpest endpoint is found when K^^^ has a value equal 
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// FOR 

♦ LIST SOURCE PROGRAM 
»ONE WORD INTEGERS 

«IOCS(CARD» 1132 PRINTER* TYPEWRITER* PLOTTER) 
C COLOR CHANGE OF AN INDICATOR 

REAL KHA* KW* KHI* MHA* MMOH* MHI • L 

VHA . 25.00 

VHI s 0«1 

MHA s 0*1000 

MMOH s OalOOO 

MHI = 0«03 

KHA = l,E-5 

KHI = l.E-12 

KW = 1«008E-1<» 

L = 1«00 

EI s l,5E/» 

EHI = le5E2 

CALL PREP (.lA* 8«* 1«* !•) 

WRITE {7»15) 

15 FORMAT ('COLOR CHANGE OF AN ACID-BASE INDICATOR*) 
CALL YAXIS (0«* 0«* .l* 10* 2) 

WRITE (7»17) 
17 FORMAT CABSORBANCE •05PH') 
CALL XAXIS (0«* 0«* 5«* 10* 5) 
WRITE (7»16) 

16 FORMAT ('MILLILITERS OF SODIUM HYDROXIDE') 
H = 1. 

1 H = •1»H , 
E » VHA»(H»»3 + KHA»(H»»2) - (KHA»MHA + KW)»H - KW»KHA*^ 



,F - VHI»(H^»3 + KHA»(H»»2) - KW»tt - »W*>KHA) .^S A. « 

"^.DEM = H*»3 + (KHA + MMOH ) » ( M»»i S ♦ {-KHA^MMOH Kl^)i»«-5 
VM9V1 =1 {.ftE + F)>D£N « • I .-ill 




KW«KHA 
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8 =fi VHI»MHI/(VHI + VMCH + VHA) ' \ 

C « L'(H KHI ) 

D » EH^»H + EI^KHI . 
A = B»q»D A 
IP (VMC^) lt.2»2, . \\ 

2 CALL FPL'OT t-2»VM0H»A) \ '\ 

5 H « H - •1»H V ^ l\ ; 

E » VHA»(H»»3 + KHA»(H»»2) - {KHA»MHA + KW)»H - KW»KMA) 
F » VHI»(H»»3 + KHA»(H»»2) - KW»H - KW»KHA) 

DEN ■ H»»3 + (KHA + MM0Hl»(H»»2) ♦ (KHA«MMOH - KW)«H - KW«KHA 

VMOH = -(E 4 F)/DEN 

B = VHI»MHI/(VHI + VMOH + VHA) 

C = L/(H + KHI ) 

D = EHI»H + EI»KHI 

A = B»C»D 

IF(50«-VM0H) 10»3t3 

3 CAL.'. FPLOT (OfVMOHtA) 
GO TO 5 

10 CONTINUE 

CALL FPLOT (1»VM0H»A) 71 



PAGE 2 

H = 1. 

21 H « •1»H 

E ~ VHA»(H*»3 + KHA*(H»»2) - (K,HA»MHA * K.W ) »h - K.W*K.HA) 
F * VHI»{H»»3 + K,HA»{H»»2) - <W«H - KW»lCHA) 

DEN * H««3 + (KHA ♦ MM0H)»(H»«2) + {KHA«MMOH - K.W)«H - KW»ICHA 
VMOH « -JE ♦ F>/DEN 
IF (VMOH) 21f22f22 

22 PH » -(ALOGiH) )/2«3C3 
P « PH/20. 

C LL FPLOT (-2tVM0H»P) 
25 ' « H - .l^H 

c - VHA»(H»»3 + KHA»(H»»2) - {KHA«MHA + KW)»H - KW«KHAJ 
F « VHI»(H»»3 + KHA»(H»»21 - KW»H - KW»KHA} 

DEN » H»»3 + (KHA + MM0H)»(H»»2) + {KHA«MMOH - KW)»H - ICW«KHA 
VMOH « -(E + F)/DEN 
PH = -(AL06{H) )/2«303 
P « PH/20. 

IF(50«-VMOH) 30»23»23 

23 CALL FPLOT (0»VMOH»P) 
GO TO 25 

30 CALL FINPL 
CALL EXIT 
END 
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mm mm. [f m acio-base micatgr 



o 
o 
a 



VHA = 25.00 
VHI = 0.1 
MHA = 0.1000 
MMOH = 0.1000 
MHI = 0.03 
KHA = l,£-5 
KMI = l.E-6 
<W = 1.008£-i^f 
L = 1.00 
EI = 1.5E4 
EHI = 1.5e2 



A 0 




Um CHANE Df AN ACID-BASE INDICATOR 



VHA = 25.00 
VHI = 0.1 
MHA = 0.1000 
MMOH = 0.1000 
MHI = 0.03 
O I KHA = l.E-5 

^ KHI = l.E-8 

KW = .1.008E-14 
L = 1.00 
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COLOR CHANGF CF AN ACID -BASE INDICATOR 



O 

o 



VHA = 25.00 
VHI = 0.1 
MHA = 0.1000 
WViOH = 0.1000 
MHI = 0.03 
KHA = l,E-5 
KHI = l.E-10 

= l.OOeE-14 
L = 1.00 
EI = 1.5E4 
EHI = 1,5E2 



O 

o 

HI 



♦ 



X 

a 



10 



8 



LiJ 
LJ 

o 

CD 
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iO'OO 20-00 30-00 

millilit[:rs t sodium hydroxide 
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WSR mANCE OF AN ACIO-BASE IMjICATCR 



o 
o 
o 



VHA •- 25.00 
VHI = 0,1 
MHA = 0,1000 
KMOH = 0.1000 
MHI = 0.03 
KHA = l,E-5 
KHI = l.E-12 
KW = 1.008E-14 
L = 1,00 
EI = 1.5E4 
EMI = 1,5E2 
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to (H*) at the basic side of this equivalence zone. 

If the indicator has an equilibrium constant 100 times too 
large or too small, then no abrupt change of colcr will be 
observed, and the color change that does occur may not take place 
at the equivalence point. 

These calculations assumed a particular set of values of 
^HA* ^HA* ^NaOH* ^ different acid, and a different set of 
concentrations, would be expected to affect the results. What 
does your chemical intuition say? How can you check your 
intuition? 

^ t ^ ^ 7^ 
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CHANGES IN THE TWO DISTRIBUTION FRACTIONS 
DURING THE TITRATION OF A MONOPROTIC ACID 

A fruitful way of thinking about the shifts in chemical 

equilibria that occur continuously during a titration is to look 

at the behavior of the ^pzdizi dl6tn,lbutlon {^xaction^ during 

that titration. The species distribution fractions aie functions 

that tell quantitatively, for each value of (H*) , how the protons 

are distributed between HA ?nd A", and how total A is distributed 

between these two chemical species. For the case of the mono- 

protic acid HA being titrated with NaOH, the two fractions are 

(HA) (A") 
and 



[A] [A] 
Note that (HA) /[A] is the fraction of all A- containing species 
present in solution as HA, and that (A") /[A] is the fraction 
present as A . To obtain algebraic equations relating these 
fractions to (H*) , we shall combine the defining eq ation for 
the equilibrium constant Yi^^ with the cons er^'^at ion equation for 
all A-containing species. Thus, Equations 8 and 10 can be 
readily combined to give 



[A] = (HA) 



1 + HA 



(H*) 



whereupon rearrangement yields 

(HA) 1 (H*) 



(H*) 



(42) 



In the same manner, we can again start with Equations 8 and 
10, this time obtaining 
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[A] = (A-) 1 + 



K 



HA 



which can be rearranged to give 



(A ) 



1 




lA) 



(43) 



1 + 



HA 



K. 



HA 



• Show that the. 6um oi the. tMo dl&tKlbtitA,on ^KCic.tlon& l6 
e.qaal to 1 . 

• Show that when (H*) Ju 6ub6tantlally gfieatex than K„., the 
pfiQ^domlnoLYit 6pz(ilz6 In solution l6 HA. 



predominant 6pz<ilz6. 

• Wfeen the tu)o distribution fractions equal one another land 



To illustrate the vie o£ distribution fractions, let us 
investigate the chemical changes involved during the titration o£ 
a monoprotic acid with a dissociation constant equal to 10"^ (that 
is, with pKjj^ = 5). FORTRAN Program 8 can be used to plot the 
titration curve and the two distribution fractions (see the output 
on pages 85 to 88). The fraction (HA)/ [A] is almost equal to 
1.0 at the beginning of the titration, meaning that in this 0.1 
molar solution of acid in water, most of the acid is undissociated. 
Then during the titration, the fraction decreases essentially 
linearly with addition of NaOH, reaching the value CO almost 
exactly at the equivalence point. The fraction (A')/ [A] begins 
almost at zero, and increaijes linearly to 1.0 at the equivalence 



• Show that when (H*) Is much less than K 



HA' 



a' Is the 



thus each Is equal to Va), what Is the numerical value o^ (H^)? 
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point. From the beginning of the titration to the equivalence 
point, each aliquot of sodium hydroxide solution has the same 
incremental effect on each of the fractions. 

Then, at thz zqalvatzncz point, thz situation markedly 
change,6. Further addition of titrant has no observable effect 
on either fraction. There is a qualitative change in the 
dependence of the fractions on volume of added titrant. Each 
fraction ceases its linear change, and remains essentially equal 
to either zero or to one. 

The fraction (HA)/ [A] reaches essentially zero on the linear 
scale; on a logarithmic scale, ve would observe this function 
to continue its decrease. In fact, for (H+) < ^ha * Equation 
1*2 assumes the limiting form (HA)/ [A] ^H+)/Kjj^, and the 
value of the fraction plunges toward zero together with the value 



This qualitative change signals a change in the whole character 
of the significant equilibria; the chemistry is different. The 
entire first portion of the titration from start to equivalence 
point can be described accurately as the reaction of NaOH with 
HA. However, at the equivalence point the solution is unable to 
furnish unreacted HA to continue the reaction. The equivalence 
point in this titration marks a sharp dividing line between two 
very different chemical situations. The most dramatic occurrence 
seen on these plots is at the equivalence point on the plot of 
pH versus Vj^^^qu* where the quantity (H*) changes by a factor of 
a million with the addition of just a few milliliters of titrant; 
this change is abrupt even on the logarithmic pH scale. It is 
this abrupt change that permits the reliable detection of the 



of (H+). 
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equivalence point in the titration. 

Note that both distribution fractions arc equal to V2 at 
the macroscopic half-equivalence point, and that the pH at that 
point is numerically equal to the value of pK^j^. 

For a substantially stronger acid, the picture drawn by 
the distribution-fraction curves is different in several signif- 
icant respects. A stronger acid (see the curves on pages 85 and 
86) is dissociated to a greater degree at the beginning of the 
titration. With an acid having pKj^^ = 5, we could say with 
justification that equimolar amounts of OH (from the NaOH solu- 
tion) and HA react following each addition of titrant prior to 
the equivalence point. With the stronger acids (acids with a 
larger value of K^^, and a smaller value of pK^^) , the slope of 
the (HA)/ [A] versus Vj^^^g curve is less 5teep, and this fact can 
be interpreted as meaning that OH" from the titrant reacts only 
partly with HA and also partly with H^. 

With very weak acids, the distribution fractions do not 
change linearly during the titration (see page 88) , noticeably 
just before the equivalence point. Some HA remains in solution 
even after the equivalence point has been passed. Endpoint detec- 
tion is then difficult and less precise, primarily because there 
is no longer a qualitative difference between the chemistry just 
before and just after the equivalence point. Note that this 
plot suggests a requirement for an abrupt pH change at the equiva- 
lence point: tho. 6topz 0^ thz distribution' {^^action ca/tve, plotted 
on a limaK 6cate. vq./l6u6 vo£ame 0^ added titrant, mast cfeange 
abfiaptty at the equivalence point. 
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The arguments of the preceding section are based on graphs 
in which (H*) is presented logarithmically, and each of the 
distribution fractions is presented linearly. You may wish to 
consider whether such presentations are appropriate for a 
faithful and informative presentation of the relevant facts. It 
is possible that this presentation distorts some relationships. 
Feel free to make your own computer-drawn graphs, plotted in 
other ways, and to draw your own conclusions about the chemistry. 

COMPUTER CALCULATION Calculation of the fractions (A") /[A] and 
OF THE DISTRIBUTION 

FRACTIONS FOR A (HA) /[A] as functions of V proceeds by 

MONOPROTIC ACID * + 

considering the concentration (H ) to be 

a master variable that connects Equations 16, 42, and 43. The 
chemical reason for this algebraic convenience is that the chemi- 
cal spevies is common to the two chemical equilibria, providing 
a chemical coupling between the two reactions in solution. For 
a particular value of (H*) , use of Equation 16 yields the paired 
value of Vjj^Qjj, Equation 42 gives (HA)/[A], and Equation 43 gives 
(A")/[A]. 

This program introduces two new FORTRAN features: writing 
and use of a stored subroutine, and storage of variables in COMMON. 
The COMMON statement has the general form 

COMMON A, B, C, ... 
where A, B, C, ... are variable names. Five variables are assigned 
to COMMON by our main program, and the values of these variables 
are available for use by the subroutine without being redefined 
in the subroutine. 

The subroutine appears first, so that it can be used later 
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// FOR 

♦ LIST SOURCE PROGRAM 

»I0CSJCARD» 1132 PRINTER* TYPEWRITER* PLOTTER) 
♦ONt WORD INTEGERS 

REAL NUMi KHA* MMOH* KW* MHA 

COMMON V»HH»KHA»MMOH»KW#MHA»VHA 
C SUBROUTINE PL0T2 

C CALCULATION AND PLOTTING OF THE DISTRIBUTION FRACTIONS 

CALL PREP (.l^* 1,, 1.) 

CALL YAXIS 0., ,1, IQ, 2) 

WRITE I7f37' 
37 FORMAT COI JN FRACTION*) 

CALL XAXIS • 1., 50» 10) 

WRITE (7#54 

56 FORMAT ( 'MI I LITERS OF SODIUM HYDROXIDE*) 
VMOH « V 
H = HH 

FHA « H/(H + KHA) 
CALL FPLOT (-2 »Vr10H»FHA) 
35 H » H - .l^H 

DCN » H»»3 •♦■ (KHA MM0H)»(H»»2) (KHA*MM0H - KW)«H - KW»KHA 

NUM « H»»3 KHA»(H»»2) - (KW (KHA»MHA))»H - KW«KHA 

VMOH « -VHA»NUM/DEN 

FHA = H/(H ■♦• KHA) 

IF(50«-VM0H) 30*33*33 
33 CALL FPLOT (0»VMOH»FHA) 

GO TO 35 
30 VMOH = V 

H = HH 

FA = KHA/(H ■♦• KHA) 
CALL FPLOT (1»VM0H»FA) 
CALL FPLOT (-2fVM0H»FA) 
U5 H ' H - »l*H 

DEN « H»»3 ■♦• (KHA ♦ MM0H)»(H»»2) (KHA«MM0H - KW)«H - KW«KHA 

NUM = H**3 + KHA»(H»»2) - (KW -f (KHA»MHAJ)»H - KW»KHA 

VMOH = -VHA»NUM/DEN 

FA = KHA/(H ■♦• KHA) 

IF (30«-VM0H} 40*43*43 
43 CALL FPLOT (0»VMOH»FA) 

60 TO 45 
40 CALL FINPL 

CALL EXIT 

END 



// DUP 

♦STORE WS UA PL0T2 
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PAGE 2 



// FOR 

»IOCS(CARD» 1132 PRINTERt TYPEWRlTERt PLOTTER) 
• LIST SOURCE PROGRAM 
•ONE WORD INTEGERS 

REAL KHA* KW* NUM* MHA* MMOH 

COMMON V»HH»KHA»MMOH.KW»MHA»VHA 
C CHANGES IN THE TWO DISTRIBUTION FRACTIONS DURING A TITRATION 

VHA ' 25.00 

MHA « 0.1000 

MMOH * 0*1000 

KHA « 1,E-10 

KW * 1«008E-14 

CALL PREP (.lA* .25» 5.5) 
WRITE 17»15) 

15 FORMAT I 'CHANGE OF THE DISTRIBUTION FRACTIONS DURING A TITRATION*) 

CALL YAXIS (0«*0«*1«*15*5) 

WRITE 17»17) 
17 FORMAT CPH') 

CALL XAXIS (0«*0«*1«*50*10) 

C CALCULATION OF THE FIRST POINT OF THE TITRATION CURVE 

H « 1, 

1 H « H - .l^H 

PH « -(ALOG(H) )/2«303 

NUM - H»»3 ♦ ICHA»IH»»2) - (KW (KHA*MHA))*H - KW»KHA 

DEN » H»»3 ♦ IKHA ♦ MM0H)»IH«»2) ♦ IKHA»MM0H - KW)«H - KW»KHA 

VMOH = -VHA»NUM/DEN 

IF IVMOH) 1»2»2 

2 CALL FPLOT (-2»VM0H»PH) 
V = VMOH 

C CALCULATION AND PLOTTING OF THE TITRATION CURVE 

HH « H 
5 H « H - .l^H 

PH = -lALOGlH) )/2«303 

DEN X H»»3 IKHA MM0H)«IH#»2) + IKHA»MM0H - KW)»H - KW»KHA 
NUM X H»»: + KHA»IH»»2) - IKW {ICHA»MHA) )»H - |CW»KHA 
VMOH = -VHA»NUM/DEN 
IF150.-VMOH) 10»3»3 

3 CALL FPLOT (0»VM0H»PH) 
60 TO 5 

10 CALL FINPL 

CALL LINKIPL0T2) 
END 



// XEO 
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when needed in the main program. 7\e subroutine is written as an 
ordinary program, with a specification of the variables the 
subroutine expects to find in COMMON. The last rtatement is END. 
The programmer decides on some name for the subroutine, and then 
uses the // DUP statement and the *STORE statement to transfer 
the subroutine into a temporary storage location on the disk. 

We then proceed to the main program, written in the usual 
manner* The COMMON statement instructs the computer to store in 
COMMON each of the five variables whenever they are defined, and 
to replace these values with new values whenever the variables 
are redefined. 

The CALL LINK statement calls the subroutine from disk 
storage into the working core memory, and transfers control to 
the first statement of the subroutine. The general form of the 
statement is 

CALL LINK (NAME) 

where NAME is the name of a subroutine already stored on the disk. 

There are two distinct plotting programs, and in each case 
there is a pause in the exacution before the plotter begins 
plotting to allow you to set the pen in the proper location on 
the paper. Careful readjustment of the pen position is necessary 
to assure that the second plot will be located on the paper 
directly under the first. 

^ ^ ^ ^ ^ 



COMPUTER CALCULATIONS USING THE 
BASIC PROGRAMMING LANGUAGE 
AND A TIME-SHARING TERMINAL 



Many of the computer calcv"* illustrated by FORTRAN 

programs can be readily performe-u using a remote terminal for 
input and output. There are some advantages, as well as some 
limitations, to use of a simple teletype or teletypewriter 
terminal. The principal advantage is the ease with which the 
user can interact with the program. Disadvantages of the simplest 
terminals include the necessity of relying on key' >ard input of 
programs and data (punched cards, while more difficult to prepare, 
are more reliable for repeat runs), and the lack of high-resolution 
plotting capability. In essence, batch processing of cards with 
an IBM 1130 system that includes a high-resolution plotter is 
rf^jj{e>ten^ than time-sharing computations from a simple terminal. 
Some illustrative programs, all written in BASIC programming 
language, are presented on the following pages to show how the 
inherent advantages of a terminal may be utilized for learning more 
about the chemistry of r onoprotic carboxylic acids. 

BASIC Program A is written to provide simulated data from a 
titration of the acid HA with a solution of NaOH. The variables 
are given short names (for example, C, ...) so that the 

program can be typed with a minimum of effort. The output is 
printed as a listing of paired values of V^^^^^ and pH. 

Let us examine this program in detail. The program is based 
on Equation 16, and this equation appears in a somewhat disguised 
form in statements 100, 110, and 1^0. Rather than use the symbols 
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Table 2 

ARITHMETIC SYMBOLS USED IN THE BASIC PROGRAMMING LANGUAGE 

+ addition 6ign 
- AubViaction 6ign 
I division 6ign 

" multiplication 6ign (do not indicate multiplication 

with just parentheses) 

t means Kai^zd to tht pou)z/L o^ 

= cqual6 6ign 

RELATIONAL SYMBOLS 

< le.66 than 

<= lz66 than OK equal to 

= equal to 

> QKzatzK than 

>= QfizatzK than ok equal to 

<> not equal to 
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directly from Equation 16, we have used certain arithmetic 
symbols from the listing in Table 2, and we have jjiven the various 
variables the following short names: 

Code for Translating Program into Algebra 

BASIC Symbol Algebraic Symbol 

" "ha 

D M 

NaOH 

E K 

w 

I change in (H^) 

H (H^) 
P pH 

Dl dtnomlnatoK 
V V 

NaOH 



Large and small numbers are written in zxpontntlal notation^ 

which means that two of the BASIC statements translate as follows 

A=1.75E-5 => K„. = 1,75 x lO"*^ 

HA 

E=1.01E-"lil => K = 1.01 X lO*"^** 

w 

Each typewritten line in the program is called a ^tattmtnt^ 
and each statement has a number. The first statements define 
various variables and constants, and state the initial conditions 
of the titration. Unless there are instructions to the contrary, 
the computer proceeds sequentially through the statements in 
numerical order according to their statement numbers. After one 

92 



10 A=1.7'jl:-'J 
20 B=2'j 



30 C=.l 

'50 0=.l 

50 E=1.0lE-li| 

60 I=.5 

70 H=l 

80 H=H"I 

90 P=-(LOGCiO)/2.303 

100 N=Ht3 + A"CHt2) - CE+CA«C))"H - CE"A) 
110 Dl=Ht3 + CA+D)"CHt2) + CA"D-E)5«H -CE-A? 
120 V=CC-ii)"N)/Dl 
130 IF V>=0 THEN 170 
li|0 X=P-7 

150 IF X>0 THEN 190 
160 GO TO 80 
170 PRINT V, TABC30),F 
180 GO TO 30 

190 PRINT "END OF CALCULATIONS" 
200 END 
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RUN TITRl 
8K 

.iy'l0772O) 

.7393181'j6 

1.607182 

3. 100 318611 

5.55210606 

9.10084005 

13.3'193'131 

17.'l06655i» 

20.52110076 

22.51121029 

23.7075563 

2I1.3366I18 

24.6638325 

2I1.830836 

24. 9151943 

24.957652 

24.979062 

24.9900347 

24.9960329 

25.0000496 

25.0040919 

25. 0101815 

25.0213661 

25.0432503 
25.0868258 
25.1740792 
25.3^94378 
25.7038263 
26.4277668 
27.9394798 
31.2461742 
39.2757244 
64.9608534 

423.050333 
END OF CALCULATIONS 

200, HORHAL EXIT 

time: 0.318 SEC. 



3.0097576;' 

3.3107333''- 

3.61170914 

3.9126849 

4.21366067 

4.51463643 

4.8I561219 

5.11658795 

5.41756372 

5.71853948 

6.01951524 

6.320491 

6.62146677 

6.92244253 

7.22341829 

7.52439405 

7.8253698I 

8.12634558 

8.42732134 

8.7282971 

9.02927286 

9.33024863 

9.63122439 

9.93220015 

10.2331759 

10.5341516 

IO.B351274 

II.136IC32 

11.4370789 

11.7330547 

12.0390304 

12.3400062 

12.640982 

12.9419577 
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statement has been executed, the statement with the next-higher 
number is then executed. However, you as the programmer may 
change the sequence of operations by using the GO TO or IF ••• THEN 
statements: 

HThe statement 

GO TO 80 

will be followed by execution of the statement numbered 80. Then 
the program will follow sequentially the statements that follow 
statement 80. 

HThe statement 

IF X>=0 THEN 190 

transfers the program to statement 190 the numerical value of 
the variable X is greater th^n or equal to zero. Otherwise, the 
program moves along in order to the statement next in numerical 
order following the IF ••• THEN statement. 

You can add spaces within a statement, almost at your own 
convenience, in order to make the typewritten statement more 
readable to you. If there is doubt in your mind about the order 
in which the computer will perform the indicated arithmetic within 
a statement, then add some clarifying parentheses; parentheses 
have the same meaning in BASIC as in arithmetic and algebra. You 
must always use a closing parenthesis with each opening parenthesis. 

Statements 80 - 180 constitute a computation loop. This 
loop is the heart of this particular BASIC program. A large value 
(larger than will be encountered in the actual titration) for (H^) 
is assumed in statement 70. Then this value is reduced to 501 of 
its previous value in statement 80, the pH value is calculated 
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for future reference, and then (H*) is substituted into Equation 
16. A value of the titrant volume is calculated, and this value 
is tested in two IF ••• THEN statements to sec if it is negative, 
or positive or zero; and if it is negative, whether the titration 
is in the acidic or basic region. Only positive values of volume 
are acceptable as having any chemical meaning. Negative values 
are encountered at the beginning of the calculations, and here 
the computation loop cycles without any output being printed. 
When (H*) gets small enough to produce a positive (or zero) value 
of the titrant volume, the cycle begins to include the PRINT 
statement, and output gets typed. Finally, when (H*) becomes 
smaller than its value in the titrant itself, the calculations 
pass outside the realm of chemical reality, the volume as cal- 
culated becomes negative, and the program is terminated. 

It is^ necessary that each program have a name. It was 
decided to call this program TITRl, and so the line 
NAME TITRl 

was typed prior to the typing of the statements of BASIC Program 
A. In order to get a listing of the program, one types 

LIST TITRl 
To execute the program, one types 

RUN TITRl 

The results of executing the program are presented on page 94. 
The entry 8k is printed automatically, indicating the amount of 
storage space that the computer has available for running this 
program. The output is a listing of paired values of V^^^^^ and 
pH for evenly-spaced values of (H*) . 
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After inspecting the typewritten output, you might decide 
that you want the simulated data to be more closely spaced, and 
that you have no need for data beyond V^^^^j^ = 23 mil. You can 
get closer spacing by changing the value of the variable I . This 
is easily accomplished by typing the line 
60 I=.7 

The result is that the former statement 60 is replaced by this 
new statement 60. One way to end the program at 23 mil is to add 
a new statement: 

125 IF V>=23 THEN 190 
Typing this statement adds it to the program. To check the pro- 
gram after making the change and the addition, we call for a 
listing of the program: 

LIST TITRl 

10 A=1.7'-jli-'j 

20 13=25 

30 C=.l 

^0 D=.l 

50 E = 1.0li£-lij 

60 I=.7 

70 H=l 

80 n=i^"I 

90 P=-(LOGClO)/2.303 

100 N=Ht3 + A"(Ht2) - (l£ + (A"C))"H - (f:"A) 
110 Dl=Ht3 + (A+D)"(Ht.?) + CA-p-::")"il -Cr."A-) 
120 V=(C-")"rO/Dl 
125 IF V>=23 TMliN 100 
130 IF V>=0 THEN 170 
IHO X=P-7 

150 IF X>0 THtifJ 190 
160 GO TO oO 
170 PRINT V, TAt;(30),P 
loO GO TO BO 

190 pr.IIiT "liUD OF CALCULATIOriS" 
?.00 EfJD 

We can then run the program by typing the instruction 
RUN TITRl 

which is the instruction to run the edited program. The output 
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of the program is as follows : 



RUN TITRl 



Ok 



3.'j(;:?10321 
3.71<''97727 
3. .17105132 

M.02672!337 
i4.lGl'j99^2 

^1.336^173^1- 
'I. '1 9 13 '1 70 3 

i\.?01005C'\ 

^1.9131596959 

5.11GC'l37^i 

5. 26 5 7 17 B 

5.^I20591G5 

5.575^!659 

5. 730 339 P5 



.33^133270,) 
.olb'loa7^l9 
.9O9'j02251 
1.430132132 
2.0^1068329 
2.8'l8J]0uuo 
3.9010i|676 

[3.23776234 
6.87i324774 
8.79210319 
10.9186577 
13.1 '109 3 87 
15.3226125 
17.3365122 
19.0926658 
20.5495796 
21.7090571 
22. 6016024 
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When a chemist is involved in a research project for the 

'j 

determination of the numerical value of an equilibrium constant, 
he or she must be concerned with the precision and the accuracy 
of that value. How many significant figures can be justified in 
reporting an equilibrium constant? The answer is not a simple 
one, but one important piece of information required for the 
answer is the answer to this second question: what is is smallest 
change in an assumed value of that causes an experimentally- 
detectable change in simulated values of Vjj^^jj and pH? We can 
easily get such information by repeating this simulation with 
different values of K^^; that is, we can rerun this BASIC program 
with a slightly-changed value of A, and check the output to see 
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if there are any significant changes in the predicted volumes for 
identical valufes of pH. The value of A is readily changed simply 
by retyping statement 10, and then the program is run simply by 
typing the RUN command. The commands and the results for a change 
of K„. from 1.75 x 10"^ to i.76 x 10"^ are as follows: 



10 A=1.7uL-'j 
l:UN TITfa 



8k 





?.9-'l260 7 


.337358772 


3. 09 7^1 01 J 5 


.620592737 


3. 2d;/ 35511 


.:)7533ii028 


3.m072:'-}1^; 


1.''13C17729 


3. 5 6 21;) 32] 


..'.05161339 
2.oG2'j6illl 


3.71697727 




3.9199 3179 


■I. 02672137 


5.2bl'i.'l359 


'■|.l."1599'i2 


6.90 37102^1 


•1.336^173'in 


.'■;.82n6l7'l7 
10.9537108 

13. 176 '19''; 9 


'';.'i:il3'i753 


^■|.G'!'''221';" 


'■; .r-oio )56''; 


15.3563779 


';.955');'j(-9 


17.366750''l 


5..''1 )0-'37'' 


19.113325 


5.2657178 


20.5703706 


5. '1205 )1 ^-5 


21.7253007 


5.575'i6-;:) 


22.6ll|00i|i4 


5. 730 33'- 9 5 



r::^D OF CALCULATIONS 

i:u(}, ;;oi:;iAL ^ixiT FRwfj Pnn . 
time: u*2oL) sue. 
Comparison of these numbers with the output listed on the previous 
page reveals that the predicted volumes differ generally only by 
t.02 m£, an amount comparable to the uncertainty in reading a 
burette. A closer look shows that the largest differences occur 
near the macroscopic half -equivalence point, suggesting that such 
mid- titration data are most useful in evaluating an equilibrium 
constant. 
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USE OF A BASIC The easiest way to use a time-sharing terminal 
PROGRAM FROM 

PERMANENT STORAGE is to call a program that someone else has 

written and placed in permanent storage. You 
don't need to know anything about programming language, and you 
don't need to know what the program is doing. Of course, there 
is the possibility that you therefore won't know what the output 
means! BASIC Program B is such a program. It performs the same 
calculations as BASIC Program A. In order to obtain the printout 
of simulated data, the user types the following two commands: 

FETCH PROGB 

RUN PROGB 

The results of this series of commands are given on page 102. The 
FETCH command obtains the program from storage, and places it in 
the user's active working area, ready for use. (Of course, this 
procedure only works if BASIC Program B has already been stored 
in your particular computer system, and if it has been given the 
name PROGB.) This program is written so that it always has some 
reasonable numbers supplied for each of the variables and constants 
However, if you wish to simulate data for your own special set of 
conditions, you will want to enter your own numbers. To do this, 
you merely retype the statements (30, 40, 50, 60, etc.) that are 
required to define your own experimental conditions. You need a 
listing of the program in order to identify these statement numbers 

Here you have a program in storage that serves as a starting 
point for your calculations. It has an output that includes some 
conventional English- language sentences that describe the set of 
experimental conditions assumed for the simulation. You can then 




10 PRINT "THIS PROGRAM SIMULATES DATA FROM A 

20 pitiirr ".'louoPROTic Acin tithation" 

;'l PkiNT 

PRIiJT 

iO A=1.7')l.-'» 

'10 fi=?'j 

•1)0 C=.l 

60 D=.l 



6S 


f-: = l. OK: 


-li! 






70 


I = .5 








80 


'1=1 
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I'RINT 
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PRINT 


"13 — 


THE 


VOLU.'IE OF 


101 


PRIUT 








110 


PRINT 


"C— 


THE 
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111 


PKINT 








120 


PRINT 


"U- 


THE 


VALUE FOR 


121 


I'RINT 








130 


PRINT 


"1: — 


THE 


VALUE FOI": 


131 


PIUKT 








li|0 


PRINT 


"I — 


THE 


DECnul'.ENT 


I'll 


PRINT 








li>o 


PKliiT 


".1— 


THE 


VALUE FOP 




PRffJT 








r;2 


PRINT 








lti3 


PRINT 









IbO PRINT "VOLUnE OF NAOH", T AB (ifO ) /'PH" 

161 PRIiJT 

162 PRIIJT 
170 M=H5M 

180 P=-(LOG(H))/2.303 

190 N=Ht>A5J(Ht2 )-(£ + ( A55C))"H-rF"A) 

200 ri = !n3+(A+D)"CHtL')+(A"D-E)"H-Cn"A) 

210 v=r(-i5)"ro/i>i 

220 IF V>=0 THEN HoO 

230 y=p-v 

:."io IF :->f. THEN nno 

2'jO GO TO 170 

260 PRINT V,TAl;(^IO),P 

270 CO TO 170 

J80 PRIfJT "END OF CALCULATIONS" 

^!90 END 
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THIS PnOGR>\M SIf-5ULATES DATA FROf A f'ONOPROTlC ACID TITRATION 



A— rm. vAi.ui: i-on (h+)(a-j/cma)= i.7'.,i:-', 

U— THE VOLUML- OF THE ORIGINAL COOH SOLUTION= 2') 

C— THL VALUE FOR THE MOLARITY OF THE COOH SOLUTIO.'s'= .1 

D- THE VALUE FOR THE MOLARITY OF THE NAOH 50LUTI0N= .1 

E--THE VALUE FOR CH+)C0H-)= 1.0lE-li| 

I— THE DECREMENT= .5 

H— THE VALUE FOR (H+)= 1 

VOLUME OF NAOH PU 




.19^f077263 

.739 318196 

1.607132 

3.1003186^1 

5.55210606 

9. 100 ail 00 5 

13.3193^31 

17.1066551 

20.52I0076 

22.5121029 

23.7075563 

2I.3366I8 

21.6638825 

21.830836 

2I.91519I3 

21.957652 

2iL979062 

21.9900 31 7 

2I.990O329 

25.0000^196 

25.00iJ0919 
25.0101815 
25.0213661 
25.0^132503 
25.OR68258 
25.17I0792 
25.3191378 
25.703O2C3 
26.1277668 
27.9391798 
31.2^161712 
39.2757211 



6.621il6677 
6.922II253 
7.223II829 
7.52I39I05 
7.825^.6931 

(1.126 3ii55''^ 
P.I273213I 
8.7282971 
9.029272-16 

9.33021863 
9.63122I39 
9.93220015 
13.2331759 
10.5311516 

10.3351271 

11.1361032 

11.1370739 

31.73805i'7 

12.0390301 

12.3100062 

12.610982 

1'''. 91.19577 



3.00975762 
3.31073338 
3.6II7091I 
3.91268U9 
1.21366067 
I.51I636I3 
1.81561219 
5.11658795 
5.11756372 
5.718539I8 
6.01951521 
6.320I91 
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modify the program to suit your own purposes by adding, changing, 
or deleting statements. Happily, the original program in storage 
remains unaffected by anything that you do to its copy. 

We have used the PRINT command in several ways. The state- 
ment 

mmm PRINT "MONOPROTIC ACID" 
results in the printing of the words MONOPROTIC ACID, followed by 
advance of the paper to the next line. The statement 
PRINT 

advancer, the paper one line. The statement 
PRINT X 

results in the printing of the numerical value of the variable x. 
The number is printed with nine significant figures, except that 
trailing zeros are not printed. The empty spaces at the extreme 
right end of the numbers printed on page 102 represent zeros. If 
a numbeV happened to have only zeros following the decimal point, 
the decimal point would be omitted, and the number would appear 
as an integer. The typewriter's tabulation capability can also be 
used by inserting the command TABCN) within the PRINT statement; 
this results in the typewriter "tabbing" across the line. The 
computer considers the typewriter set with the left-hand margin at 
column 0. The command TAB(N) causes the next character to be 
typed in column N. N can be either a number or an arithmetic 
expression. 

PLOTTING WITH AN The listings of numbers that result from BASIC 
OUTPUT PRINTER 

USING BASIC Programs A and B have many more significant 

figures than are needed for any purpose involving a comparison 

ERIC 



with experimental data. Despite this superabundance of digits, 
however, the lists do not necessarily convey the full signifi- 
cance of the calculations. An overall picture of the calcula- 
tions is often seen better if the data are graphed. A graphic 
presentation can be produced by a conventional typewriter or 
teletype printer. Such a graph is rather crude, since the coord- 
inates of each point must of necessity be positions on the paper 
where a character can be typed. Notwithstanding this limitation, 
such plots can be informative, and we shall present two programs 
that produce graphical output of a titration simulation. 

BASIC Program C illustrates a simple plotting program. Very 
similar to BASIC Program A, this program gives output via the 
statement 

170 PLOT 10+V="''", 10="+" 
The PLOT statement is executed from right to left, so that the 
first command is .to plot 10="+". This means that the typewriter 
is to type a + sign in column 10. And then, without advancing 
the paper, the command to plot 10+V=""" is executed. This means 
that tie typewriter is to plot the symbol " in column 10+V. 
Finally, when all the characters have been typed, the paper is 
advanced to the next line. In order to have pH plotted along the 
vertical axis, the c^>lculation loop has been instructed to change 
pH by increments of 0.5 pH units, thus making each print line on 
the graph correspond to 0.5 pH units. The horizontal axis is the 
^NaOH ^® used 50 spaces (50 character columns) for this 

graph, and the volume range was from 0 to about 50 mil, so we 
were able to let each space equal one volume unit. Otherwise, we 
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LIST PROGC 

10 A=1.75E-5 
20 n=?.5 
30 c=.l 

iiO D=.l 
60 I=.5 

70 p=l 
80 p = P+I 

90 H=10t(-P) 

100 N=Ht3 + A"CHt2) - (E+(A':C))"H - fE»A) 
110 lJl=Ht3 + (A+lj)5!(Ht2) + (A"D-l£)"H - CF"A) 

120 v=((-a)«N)/ni 

130 IF V>=0 THcN 170 
mo X = P-7 
150 IF X>0 THtN 190 
160 GO TO 80 

170 PLOT 10+V=""'», 10="+" 
130 GO TO 80 

190 PRINT "EiND OF CALCULATIONS" 
200 END 

RUN PROGC 



8K 

+ 

+ " 

+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 

FMD OF CALCULATIOiNS 

^00^ IWKMAL cXIT FROM PROG. 
time: 0.218 SEC. 
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would have been required to multiply the variable by an appropri- 
ate scaling factor so that the graph would fit into the allotted 
space on the paper. 

BASIC Program D is written for permanent storage, and has 
several useful features • The output begins with a title, and a 
listing of the numerical values of the various constants that 
characterize the particular acid being titrated and the experi- 
mental conditions of the titration being simulated. The actual 
plot has labelled axes, and an indication (with the printing of 
the symbol X) of the half -equivalence point. This program can 
be called from storage and used without understanding the various 
plotting statements, modifying computation statements as desired 
in the same manner as described on the previous pages. However, 
there are some additional programming features which may be of 
interest. 

Use of a 6iibKoatjint is illustrated in statements 424 to 430. 
These statements are a miniprogram to plot an X at the half-equiv 
alence point. The final statement is RETURN, Execution gets 
transfered to the subroutine by the command 
GOSUB k2k 

After completion of the calculations or instructions of the sub- 
routine, the execution returns to the statement immediately fol- 
lowing the GOSUB statement. 

For this graph, it was decided to use a scaling factor to 
expand the horizontal axis. Can you find where the scaling facto 
is inserted? 

In order to get the vertical axis labelled with pH values, 
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ID 



10 A=1.7'JU-'.> 
20 11=2 'j 
30 C=.l 

'10 n=.i 

50 E=l.008t-in 
60 l=.5 
65 H=l 

70 PRINT "PLOT OF A SIMULATED MONOPPOTIC ACID TITRATION CURVE 

71 PRINT "V^ITH (m-)(A-)/(HA)=";A 

72 PRINT 

73 PRINT 

80 PRINT "A— THE VALUE FOR (H+)(A-)/ (HA)="; A 

81 PRINT 

90 PRINT "li— THE VOLUME OF THE ORIGINAL COOH SOLUTION=";B 

91 PRINT 

100 PRINT "C— THE VALUE FOR THE MOLARITY OF THE COOH SOLUTION=";C 

101 PRINT 

110 PRINT "D— THE VALUE FOR THE MOLARITY OF THE NAOH SOLUTION=";D 

111 PRINT 

120 PRINT "E— THE VALUE FOR (H+) (0H-) ="; E 

121 PRINT 

130 PRINT "I— THE decrements"; I 

131 PRINT 

110 PRINT "II— THE VALUE FOR (H+)=";H 
mi PRINT 
150 PRINT 
160 PRINT 

170 PRINT TAb (15),"MILLI LITERS OF SODIUM HYDROXIDE" 
180 PRINT "PH" 

190 PRINT TAB C10),"0 5 10 15 20 

200 PRINT TAB (10)/'+ + + + + + + + + 

210 FOR V.'= 1 TO 10 
220 PRINT TAB(10)/'+" 
230 NEXT W 

235 X=0 

210 H=H5JI 

215 IF H>A THEN 250 

216 GOSUo 12 1 

250 P=-(LOGCH))/2.303 
260 P1=INT(10"P)/10 
270 P2=INT(P) 
280P=P/2 

290 ll=Ht3+A"(Ht2)-(E+(A"C))"H-(E"A) 
300 Dl=Ht3+(A+D)"(Ht2)+(A"D-fc:)"H-rii"A) 

31U ■•/= a-ny-: 0/ui 

3201 F V>=0 THEi.36U 

330 >;=P-Y 

3IU IF A>0 THEN 12J 
350 GO TO 210 

360 IF Plol'2 THuN loo 

361 IF PI <=J TilEU370 

362 PRINT PI; 



303 PLOT u+C2"V)+.'j=""'M)="+" 

3Cl\ GO TO ?.i\0 

370 PRINT Fi; 

330 PLOT 7+C2"V)+.i5=""'Sr)="+" 

390 GO TO P.Ho 

nOO PLOT 10 +C2"V)+.'j=""",9="+" 

^110 GO TO 2 '10 

■^20 PRINT "fZnn OF CALCULATIONS" 

^121 GO TO iJ3I 

i|2i| IF X=l TMErj i|30 

il25 X=l 

H2G P=- CLOG (A) )/,"'. 303 

^J27 Vl=.5"«"C/f> 

iJ28 PLOT 9+C2"Vl)="X",9="+" 

A 30 RHTURN 

i|31 fi4D 

OUTPUT 

PLOT OF A SIMULATED MONOPROTIC ACID TITRATION CURVE 
WITH (H+)CA-)/CHA)= 1.75E-5 

A— THE VALUE FOR CH4-)CA-)/(HA)= n.75t£-5 

B— T!IE VOLUME OF THE ORIGINAL COOH SOLUTION^; 25 

C— THE VALUE FOR THE ."OLARITY OF THI£ COOH SOLUTION= .1 

D— THE VALUE FOR THE fOLARITY OF TilE NAOH SOLUTIOH= .1 

E— THE VALUE FOR CH+)COH-)= l.OOSE-lH 

I— THE DECREM':rjT= .5 

H— THE VALUE FOR CH+)= 1 



OUTPUT CONTINUED ON NEXT PAGE 
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f:iLLiLiTLr:r> of .>oi)ium 'iydroaI i.r 
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+ 

+ Y 
+ X 
+ X 
l:UD OF CALCULATIONS 

^131/ NOKMAL FXIT FROM PROG, 

f^i: o.ti98 sec. ;n:8 
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it was necessary to have integral values of pH available. To do 
this, use was made of the INT(X) command to convert a variable into 
its next- lowest integer (the number rounded off to an integer by 
simply discarding all digits following the decimal point). This 
program is pre* ^nted to illustrate the sorts of programming 
niceties that can be used to dress up the output of a printer- 
generated graph. You may have other ideas. 

CONFRONTATION OF BASIC Program E provides a point-by-point 
SIMULATED DATA AND 

EXPERIMENTAL DATA comparison of experimental data from a 
FOR A TITRATION 

titration experiment with the predictions of 
a mathematical model for that experiment. The program is written 
for permanent storage. When run, the program provides for a 
dialogue between computer and user, with instructions coming from 
the computer in the form of English- language questions and state- 
ments. In the interest of simplicity, and in the interest of 
getting the user involved in the actual program, this interactive 
dialogue has been kept to a rather minimal level, but you may 
wish to write a more extensive program with more conversational 
interaction. 

The program begin, with a dimension statement 
10 DIM A(50),B(50) 
that reserves storage space for 50 paired numbers A^,B^; that is, 
for 50 subscripted variables that will be identified with the 
experimental values of pH and V^^^^jj. Within the program, the 
subscripted variables are written as A(X) and B(X), where identical 
values of X means that the two variables are paired. Statement 
60 prints a question that asks the user to type a number that 



tells the computer how many sets of pairs {pH,Vjj^Qjj} are to be 
entered. The following statement 
70 INPUT N 

orders a pause in the execution of the program until a number 
has been entered by the user at the terminal. That number becomes 
the value of the variable N, and the statements numbered 90 through 
140 then ask for values of pH and Vj^^^jj a total of N times, storing 
the entered numbers as an array of the form 

ACl), BCD 

AC2), BC2) 

AC3), BC3) 

• • • 

ACN), BCN) 

The program now proceeds to take each of these pairs in turn, 
calculating a predicted value of the volume that corresponds to 
the experimental pH value, and then printing the predicted and 
experimental values of the titrant volume for comparison. 

In order to perform the calculations, it is necessary to have 
values for the two equilibrium constants, the two molarities, and 
the initial volume of acid. Tfce^e numbe.u have. aiAzady been 
enie^ed by 6tatQ.mQ.nt 160. You might wish to alter the program by 
deleting statement 160 from the program, allowing yourself the 
option of entering these numbers during the execution of the 
program. You would need a statement prior to statement 150 that 
printed out an appropriate question to indicate the form in which 
these constants are to be entered from the keyboard. 



vnr 




10 DIfl AC:>0'),liC >i'>') 

;H) i'iunt "fins pi;oc.i'>A/. coriHi-.orjT-.; i.Mi'iinif.iuTAi. ivma" 

]{) I'UIiiT "'.JiTil Tin: PRin^ICTIOIiS OF A t'ATIIKHAT I CAl. r'OfMI " 

3'j iMiIfjT "Hoi; tui-. TIinATloN OF A ;'or«!oiM'wTl(: AC in" 

'ID VtAtn 

i>o Puii;T 

60 PRINT ":IOU |;ANY SETS OF lIXI'uR IJilinTAL DATA '..'ILL YOlJ .JF USIIIC"; 
70 INPUT U 

80 pnI;^lT 

JO FOR X=l TO t' 
100 PRir4T "i'H="; 
110 INPUT A(X) 

120 PIUiiT "VOLUf!E OF 1!A0H="; 

130 INPUT iKX) 

1^10 liEXT X 

130 READ A,J,C,D,E 

160 DATA 1.75,25/0.1,0.1,1.008 

170 A=A"10t(-j) 

180 r:=E"10tC-liJ) 

190 PRirjT "THE VALUE FOR (H+)CA-)/(HA)="; A 



200 
210 
220 
230 
2^10 
250 
260 



"THE VOLUP.E OF THE ORIGINAL HA SOLUT I0f4="; 15 

r.OLARITY OF THE MA nOl.tlTI ON=";C 



"THE VALUE FOu TiiE 
"THti VALUE FOR THE 



r.OLAniTY OF THE HAOH 50LUTlOri=";f^ 



THE VALUE FOR CMCO\i-')=" }il 



PRINT 
PRINT 
PRINT 
PRINT 
PRIi4T 
PRirjT 
PRINT 
270 PRINT 
280 PRINT 
290 PRINT 

300 PRINT TABC20)/'CALCULATED",TAiK'lO)/'EXPERHlENTAL" 

310 PRINT "PH",TAB(20>/'VOLU:iE OF I1AO.<",TAUC^O), "VOLUME OF NAOH" 

320 PRINT 

3 JO PI; I NT 

3^»0 PRINT 

350 FOR X= 1 TO il 

360 T=AC'0 

370 ;l=10tC-T) 

380 H=-(LOGOO)/2.303 

390 N=Ht3+A"(Ht2)-(E+CA"C))"H-(E5JA) 

i|00 Dl=Ht3+CA+D)«CHt2)+CA"D-i;)'-H-CE"A) 

'!10 V=((-o)''iO/Ul 

'120 IF V>=0 THEN HGO 

'130 X=P-7 

H1\0 IF X>0 THEN ^180 
'J 50 no TO i|70 

l\60 PRINT A(X),TAlK20),V,TAB(i|0),ii(X) 
i»70 NEXT X 

i|80 PRINT "END OF CALCULATIONS" 
^190 END 
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THIS PROGRAM CONFRONTS EXPERIMIZNTAL DATA 
WITH THE PREDICTIONS OF A l ATilEflATICAL I'ODEL 
FOR THE TITRATION OF A HONOPROTIC ACID 



HOW J'ANY SETS OF EXPEP.IHEMTAL DATA V/ILL YO!' OH unlHG? 



PH=? 

VOLUME OF flAOH=? 
.19 

PH = ? 

3.31 

VOLUME OF imOii=? 

PH=? 
3.62 

VOLUME OF I-MOH=? 

1.73 

PH=? 

3.89 



VOLUME OF UAOils? 

3.12 

PH=? 

n.34 

VOLUME OF fJAOH=? 
5.55 

THE VALUE FOR CM+) CA-:)/(HA)= 1.75E-5 

THE VOLUME OF THE ORIGINAL HA SOLUTIONS 25 

THE VALUE FOR THE MOLARITY OF THE HA SOLUTION= .1 

THE VALUE FOR THE MOLARITY OF THE fJAOH SOLUTION= .1 

THE VALUE FOR (H+)COM-)= 1.008E-li| 



PH 



CALCULATED 
VOLUME OF NAOH 



EXPERIflE.'JTAL 
VOLUME OF NAOH 



3. Oil 

3.31 
3.62 
3.89 
H,3^ 



.2^104817^17 
.736393^1^11 
1,635901^6 
2.95390902 
6.90691832 



.19 
.7i| 

1.73 
3.12 

5.55 
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You can change just the value of A (that is, just the value 
of K„.) by adding a statement of the form 
l6l A=:1.76 

Note that the exponent portion of each equilibrium constant is 
added later, in statements 170 and 180* This procedure for 
entering the constants is just one of many possibilities, and 
you should write the program in the way that is most convenient 
for you. \ 

A sample of the user — computer dialogue, as well as the 
calculated comparison of model and data, is presented on page 
113. 

Certain data points , taken from the region of the titra- 
tion curve in the vicinity of the half -equivalence point, are 
most sensitive to variations in the value of K^^^, and these are 
the only points that need to be considered when searching for a 
value of that equilibrium constant. A final comparison should 
include data from throughout the whole titration curve. 

You may wish to write a program that allows you to try a 
whole series of different values of K^^^, while keeping the exper 
imental data stored for repeated comparisons. Another worthwhil 
refinement would be to have the program find the difference be- 
tween calculated and experimental titrant volumes, printing out 
that difference as a part of the output. 

^ ^ ^ ^ ^ 
? ? ? ? ? 



USE OF A WANG 700 PROGRAMMABLE CALCULATOR 
IN EVALUATING AN EQUILIBRIUM CONSTANT 

The Wang 700 calculator is a miniature computer with a 
limited memory (a limitation) and with extremely simple input- 
output capability (an advantage). Such a calculator can be of 
great use to the chemist who is trying to find the numerical 
value of an equilibrium constant that brings the predictions 
of a mathematical model into agreement with a few experimental 
data points clustered about the half -equivalence point of a 
titration. Wang Program a performs such calculations. The 
user is closely involved with the actual calculations, having 
close control over the numbers being used in the calculations. 
The user A.ntzKact6 with the calculator and with the calculations 
in a meaningful and effective manner. 

The Wang 700 calculator permits the user to perform arith- 
metic in the same manner as with any electronic calculator, by 
entering numbers and depressing keys labelled +, -, ^, and 
so forth. Numbers at any stage of the calculations can be stored 
in one of 122 data registers of the calculator's core memory. 
The^ whole sequence of operations that constitutes a program can 
be stored in memory so that the calculator can then execute the 
steps automatically. Finally, this program can be stored for 
future use on a magnetic- tape cassette. 

This program is based on Equation 16. The strategy is 
quite simple: we shall store (in registers 00 to 08) the various 
terms that need to be added together to give name.KatoK and 
danomlnatoK^ then perform a division at the very end. Since the 
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a 



STEP 


KEY 


CODE 


STEP 


KEY 


CODE 


00 


MARK 


0408 


33 


X 


0602 


01 


1 


0701 


34 




0605 


02 


CLEAR X 


0715 


35 


STORE DIRECT 


0404 


03 


STOP 


0515 


36 


REGISTER 04 


0004 


04 


CHANGE SIGN 


0711 


37 


RECALL DIRECT 


0405 


05 


10* 


0613 


38 


REGISTER 03 


0003 


06 


STORE DIRECT 


0404 


39 


STORE DIRECT 


0404 


07 


REGISTER 00 


0000 


40 


REGISTER 05 


0005 


08 


STORE DIRECT 


0404 


41 


CLEAR X 


0715 


09 


REGISTER 02 


0002 


42 


STOP 


0515 


10 




0713 


43 


X DIRECT 


0402 


11 


STORE DIRECT 


0404 


44 


REGISTER 05 


0005 


12 


REGISTER 01 


0001 


45 


CLEAR X 


0715 


13 


STORE DIRECT 


0404 


46 


STOP 


0515 


in 


REGISTER 08 


0008 


. 47 


X DIRECT 


0402 


15 


X DIRECT 


0402 


48 


REGISTER 01 


0001 


16 


REGISTER 02 


0002 


49 


X DIRECT 


0402 


17 


CLEAR X 


0715 


50 


REGISTER 03 


0003 


18 


STOP 


0515 


51 


CLEAR X 


0715 


19 


X DIRECT 


0402 


52 


STOP 


0515 


20 


REGISTER 08 


0008 


53 


STORE DIRECT 


0404 


21 


STORE DIRECT 


0404 


54 


REGISTER 07 


0007 


22 


REGISTER 03 


0003 


55 


RECALL DIRECT 


0405 


23 


+ 


0604 


56 


REGISTER 02 


0002 


24 


RECALL DIRECT 


0405 


57 


+ 


0604 


25 


REGISTER 00 


0000 


58 


RECALL DIRECT 


0405 


26 


X DIRECT 


0402 


59 


REGISTER 08 


OOOB 


27 


REGISTER 03 


0003 


60 


+ 


0600 


28 


+ 


0600 


61 


RECALL DIRECT 


0405 


29 


1 


0701 


62 


REGISTER 04 


0004 


30 


4 


0704 


63 




0601 


31 


CHANGE SIGN 


0711 


64 




0605 


32 


10* 


0613 


65 


STORE DIRECT 


0404 



til 



STEP 


KEY 


CODE 


66 


REGISTER 06 


0006 


67 




060il 


68 


RECALL DIRECT 


0il05 


69 


REGISTER 01 


0001 


70 


+ 


0600 


71 


RECALL DIRECT 


0il05 


72 


REGISTER 03 


0003 


73 


+ 


0600 


7i* 




0605 


75 


STORE DIRECT 


OHOH 


76 


REGISTER 09 


0009 


77 


RECALL DIRECT 


0il05 


78 


REGISTER 06 


0006 


79 


i 


060il 


80 


RECALL DIRECT 


0il05 


81 


REGISTER 05 


0005 


82 


— 


0601 


83 


RECALL DIRECT 


0il05 


8i| 


REGISTER 09 


0009 


85 




0603 


86 


RECALL DIRECT 


0il05 


87 


REGISTER 07 


0007 


88 


CHANGE SIGN 


0711 


89 


X 


0602 


90 


STOP 


0515 


91 


END PROGRAM 


0512 



SAMPLE OUTPUT OF THE PROGRAM 

= 1.75 X 10-5, Mj,^Qg« 0.1000, = 0.1000, V^^ = 25.00 

pH = 3.009, Vjj^Qg = 0.192--' 
3.310 0.736 
3.611 1.602 
3.912 3.091 
4.213 5.538 
4.514 9.081 
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user's operating instructions 



PROGRAM STEP INVOLVED 



depress SEARCH, then 0 
depress GO 

enter value of pH, then GO 

enter value of K„., then GO 

enter value of M„., then GO 

enter value of M., ^„ , then GO 

NaOH ' 

enter value of V^^^, then GO 



00 

goes to 03 and stops 

goes to 18 and stops 

goes to 42 and stops 

goes to 46 and stops 

goes to 52 and stops 

goes to 89, stops, with 

displayed value of 

V 

NaOH* 



STORAGE LOCATIONS AT TIME OF FINAL CALCULATIONS 


Register 00 




Register 01 


NaOH*- ^ 


Register 02 




Register 03 




Register 04 




Register OS 




Register 06 




Register 07 




Register 08 
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calculations are designed for data in the region of the half- 
equivalence point where predicted values of V^j^^^ are insensi- 
tive to the assumed value of K^, we shall always use the value 
= 1.00 X io-^*» 

* 

Certain keys on the Wang keyboard need an explanation, for they 
are used to move numbers in and out of the storage registers. 

STORE DIRECT. This key is followed by a register key, 
and the result is ' it the number that is displayed in the 
x-register on the calculator (displayed in nixie lights) is 
stored in the indicated register. The number also remains in 
the x-register, available for further computation. 

RECALL DIRECT. When this key is followed by a register 
key, the contents of that register become displayed in the 
x-register, ready for further computation. The number remains 
stored, and can be recalled as many times as desired. 

X DIRECT. When followed by a register key, this command 
multiples the contents of that indicated register by the 
contants of the displayed x-register, and stores the product 
in the indicated register. There is now a new number in the 
storage register, the result of the multiplication. The 
displayed number in the x-register remains unchanged. 

+, -, T. Adds, subtracts, multiplies, or divides the 
numbers displayed in the x- and y-registers, and places the 
result in the y-register. In subtraction, the order is y-x; 
in division the operation is y/x. 

+. Places the contents of the x-register into the y-reg- 
ister, leaving the x-register unchanged. 
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+. Places the contents of the y-register into the x-reg- 
ister, leaving the y-register unchanged. 

Interchanges the contents of the x- and y-registers. 

If you want to become more familiar with programming and 
using the Wang calculator, try to modify Wang Program a to get 
some of the following features: 

• W^^tt a pfiOQfiam that pzfimlt^ znt^y a^^amed valuz6 
tha ^qultlb^4.um constant K^. 

• W^^ie a p/iog^am that Kttain^ tht valaz6 oi ^^^^^^^ 
^NaOH* that thz6z valuz6 do not have, to be entered each 

timn a neto valuz o^ pH ok Kjj^ cho^an. 

• ia^ltz a pfiOQ^am that attom the. znViy o^ KzpQ.atzd valae.6 
0^ pH (^Uth the. 6amz valuz6 o^ all tha constants, 60 that 6Q.vQ.Kal 
pizce6 oi zxp^Klmzntal data can be compared ialth the predictions 
0^ the model u)lthout re-entry OjJ constants. 

Output from the Wang is in the form of a display of nixie 
lights, so the user must write down these numbers before pro- 
ceeding to the next calculation. While this may appear to be 
more demanding than simply watching a high-speed printer or a 
Selectric typewriter produce numerical output, it does have an 
advantage: the user becomes more involved with the calculations, 
and has time to think about the meaning of the numbers. This 
can often be a blessing; mistakes can be caught earlier, and 
misconcenptions can be corrected earlier. 
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Chapter 2 

Titration of a 
Diprotic Acid 



The chemical model for the titration of a diprotic carbox- 
ylic acid involves an additional chemical equilibrium. An 
additional chemical equation appears in the chemical model which 
in turn introduces an additional algebraic equation into the 
mathematical model. The equation for the titration curve thus 
becomes more complicated. Once we have derived the equation for 
the titration curve, however, this complication need cause no 
problems or inconveniences in making computer calculations. 
Indeed, instead of complications, titration of a diprotic acid 
provides additional insights into the chemistry of coupled reac- 
tions and competitive equilibria. 

We shall investigate the analysis of experimental titration 
data for diprotic carboxylic acids by again confronting labora- 
tory titration data with the predictions of our model. We shall 
also investigate several aspects of the model by performing com- 
puter simulations of the model foravariety of different situa- 
tions . 

There are many common diprotic carboxylic acids, a few of 
which are presented in Table 3. Many are major constituents of 
certain plants or food products, or of significant biochemical 
importance. Most of the low-molecular- weight dicarboxylic acids 
are sufficiently soluble in water so that a titration curve can 
be obtained without difficulty. Some of these acids have been 
studied in great detail, but the equilibria for others have been 
investigated only in passing, and some research topics of interest 
can be found within this class of rather simple compounds. 

We shall begin this chapter with a statement of the chemical 
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Table 3 

STRUCTURES OF SOME DIPROTIC CARBOXYLIC ACIDS 




T-C'-Ci malonic acid 



H H 

Tc-c-c-cT 

H H 



H H 

^-c-c-cC 



H H H 



H-'O^ { I ( 
H H H 



succinic acid 



H 

I tAanA isomer: ■ fumaric acid 

~c=c-cf 

H^G^ I D'^H cl6 isomer: maleic acid 

H 



malic acid 



-C-C-C-CL glutaric acid 
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and mathematical models that we shall use for simulating diprotic 
acid equilibria in solution. 

THE CHEMICAL MODEL Any one of the representative diprotic car- 
FOR OUR SIMULATION 

OF THE TITRATION boxyJic acids represented by structural form- 
OF A DIPROTIC ACID 

ulas in Table 3 can be represented as an 
associated acid H2A which can lose one or two protons in solution 
to become anions HA" or A~. In a solution prepared by mixing 
quantities of the acid, the salts NaHA and Na2A, sodium hydroxide, 
and water, the species present (we shall assume) are 

H2O, H2A, HA", a", Na^, H'^, oh" 
We require a sufficient number of chemical equilibria in our model 
to permit the interconversion of all A-containing species, as well 
as the dissociation of the solvent water. We shall assume the 
following three chemical equilibria: 

H2O 4! H'^ + OH" (44) 
H2A i! HA" + H'^ (45) 
HA i! a"" + h'^ (46) 
This listing of the chemical species assumed, together with this 
listing of the chemical equilibria assumed, constitutes the cfeem- 
Icat modzl iofi thz ^y^tzm. 

THE MATHEMATICAL MODEL Associated with each of the chemical 
FOR OUR SIMULATION OF A 

DIPROTIC ACID TITRATION reactions of the chemical model is an 

equilibrium- constant equation. Thus we 
write the algebraic equations 

\ = (H^)(OH") (47) 

(H^)(HA") 
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(H*) (A ) 

K„,_ = — (49) 

(HA ) 

The requirement of electrical neutrality of the solution gives us 
the equation 

2(a') + (HA') + (OH') = (H*) (Na*) (50) 
Notice the coefficient "2" that appears so that two negative 
charges will be counted for each A" anion in solution* Finally, 
we include two conservation equations in our mathematical model: 
[A] = (H2A) + (HA~) + (a") (51) 
[Na] = (Na*) (52) 
Equations 47 — 52 suffice to permit elimination of all micro- 
scopic species concentrations except foi (H*) , giving us an 
algebraic equation relating [A] and [Na] to (H*) • For the 
particular experimental situation of the titration of the acid 
HjA with the base NaOH, it is convenient to change variables by 
introducing the two equations 

V„ . 

[A] = (53) 

V + V 
HaA NaOH 

[Na] = ^NaOH^aOH ,543 

V + V 

HzA NaOH 

Equations 47 — 54 constitute a set of simultaneous algebraic 
equations that can be combined to give an equation that relates 

the two experimental variables (H*) and V^^^^jj. The equation is 

Yiamo-fLatofL 



V = -V 
NaOH HA 



(55) 



dQ.n.omlnatOK 

where the quantities namzKatoK and de.n.omln.a.to'L are given by the 

relations ^ 
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Equation 55 is the heart of all the computer calculations in 
this chapter. 

• Vt^iom thz cilQzbfiCL that ton^^otldata/^ thz Iniofimntlon oi 
Equations 47 — 54 Into Equation 55. 

• Compare. Equations 16 and 55. On the. ba6l6 o^ thl6 com-' 
pa^l^on, pxe^dlct the. iom {and -6ome o^ all oi the. do^talli] oi 
thz co^^Q^pondlng equation ioK a t^lpn^otlt acid. Hou) miild 
you vtnJL^y youK pKzdlctXori^ 
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HKtUiCllUNS Uh IHh MuOtL hOR THt 
TITRATION OF A DIPROTIC ACID 

Again we shall present two different methods of using an 
IBM 1130 computer to simulate the titration of an acid with a 
base in aqueous solution. Comput^i programs will be presented 
with output both as a tabulation of numbers and as a graph of 
titrant volume versus pH. We shall assume, for purposes of 
discussing these programs, that values of the various equilibrium 
constants have been obtained in some way. In the next section, 
we shall show that t^ie *pa^t>ca£-eqa>cva£ence point mctfiad, an ex- 
tension of the half-equivalence point method, is an adequate 
method for estimating a trial value of dissociation constants 
for diprotic acids in most cases, and an excellent method in 
some cases. 

PREDICTIONS OF THE MODEL FORTRAN Program 9 utilizes Equation 

PRESENTED AS A TABLE OF _ ^ . , , . r 

pH AND VOLUME VALUES 55 for its calculation of the titrant 

volume that corresponds to various values of (H*) . The program 

begins with (H*) = 1.00, and then systematically reduces the 

value of (H*) with the statement 

H = .8*H 

a statement that replaces the value oc the variable H by 80% of 
its value. When the assumed value of (H'^) is larger than possible 
in the titration. Equation 55 yields a negative value of the 
titrant volume. Such a situation is chemically impossible, and 
we reject all such pairs of ^H^) ,V^^q^ values with the statement 
!F CVMOH) 2, 20, 20 
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// JOB 



SV1007#S51 



// FOR 

« LIST SOURCE PROGRAM 
«IOCS(CARD* 1132 PRINTER) 



C 

c 
c 



TITRATION CURVE 

PREDICTIONS OF THE CHEMICAL MODEL 
DIPROTIC ACID 



REAL KHA* KHHA* KW* NUM* MHA* MMOH 

VHA « 25.00 

MHA ' 0.0500 

MMOH > 0*1000 

KHHA « l.E-6 

KHA » l,E-7 

KW ■ l«008E-l<f 

H > 1« 

1 H » .a^H 

PH ■ -(ALOG(H) J/2.303 

NUM » H##4 ♦ KHHA»(H#»3) ♦ (KHHA»KHA - KHHA»MHA - KW)«(M»»2) 
»-(2»KHHA»KHA#MHA ♦ KW«KHHA)*M - KW#KMA#KHHA 

DEN ■ H##4 (KHHA MM0H)*(H#«3) (KHHA#KHA KHHA«MMOH - KW) 
##<H»»2) ♦ (KHHA#KHA»MMOM - KW»KHHA)»H - KW»KHA»KHHA 

VMOH ■ -VHA»NUM/DEN 

IF (VMOH) 2* 20* 20 

2 IF (PH - 7.00) 1» 1» 30 
20 WRITE (3*10) VMOH* PH 
10 FORMAT (F10«3* F10«3) 

GO TO 1 
30 CALL EXIT 
END 



// XEQ 
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which transfers the program back to statement 1 (by way of 
statement 2) if the volume is negative, and allows output of the 
calculated value of titrant volume only if that volume is zero 
or positive. Statement 2 checks to be sure that the negative 
volume has been encountered early in the titration; if the pH 
value has passed 7 and the titration has been proceeding, a 
negative volume signals that we have assumed a value of (H*^) 
that is less than chemically possible, that we have assumed a 
value of (H^) lower than the value in pure titrant. Such an 
event is the sign that all meaningful calculation has been 
completed, and that it is time to stop. Transfer of the program 
to statement 30 calls exit and end of the program. Output of 
the program, with the particular set of equilibrium constants 
and initial concentrations and volume give in the program listing, 
is presented on page 124. 

• W^itz a compute.^ p^og^am that pK(UQ.nt6 output <u evenly- 
spaced pH vatue6. 

• HzKZ^6 a chattznge! W^ite a computer pKogKam that glvz^ 
output a6 zvQ.nty'6pa(ie.d V^^^^^ vatue.6. 

Note that the algebraic expression for NUM is too long to 
fit on a single punch card, and so it was continued to a second 
card. A character punched in column 6 indicates that such a 
card contains a continuation of the previous card. The statement 
gets typed on the program listing as a two-line statement, but 
gets executed as if it were all printed on a single (extra- long) 
line. 

As expected, the equivalence point occurs at 25.00 ml when 
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25.00 mZ of a diprotic acid of molarity 0.0500 moles/liter is 
titrated with a sodium hydroxide solution of twice that molarity. 
The criterion for equivalence is that twlct a6 many molz6 oi 
ba^c ma6t be addzd to a givtn numbeJi OjJ mott6 oi acid. Why 
"twice"? 

ON SEEING TWO ENDPOINTS The number of endpoints that can be 
IN THE TITRATION OF A 

DIPROTIC ACID observed during the titration of a 

diprotic acid with a solution of sodium hydroxide depends on the 

numerical values of the two dissociation constants of the acid, 

as well as on the values of the molarities of the acid and 

titrant solutions. On the following pages, we shall examine the 

results of some numerical calculations simulating the titration 

of several hypothetical acids in order to gain insight into the 

effect of the values of K„ . and K„. on the shape of the titra- 

H2A HA" 

tion curve. We shall focus our attention on the question: "Hou) 
many tndpoint6 can be Ktcognlztd during tht titration?" These 
calculations are based on Equation 55, and are performed using 
FORTRAN Program 10. Output for various values of the two equi- 
librium constants is presented on pages 128 — 133. 

Inspection of the simulated titration curves leads one to 
propose some generalizations of the following sort: 

Uln order to get a sharp first endpoint, there must be a 
large separation of the numerical values of the K's, probably of 
the order of at least a factor of 10^. 

UTo achieve a sharp second endpoint, it is necessary to 

reduce the value of (H*) to a value much less than the numerical 

value of K„... 

HA 
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// JOB S71007tS51 



// FOR 

» LIST SOURCE PROGRAM 

»I0CS(CARD» 1132 PRINTER* TYPEWRITER* PLOTTER) 
»ONE WORD INTEGERS 

C ON SEtING TWO ENDPOINTS IN THE TITRATION OF A DIPROTIC ACID 

t.*- KHA* KHHA* KW* NUM * MHA* MMOH 
V^. 25.00 
MHA « 0*0500 
MMOH « 0.1000 
KHHA = l.E-A 
KHA » l.E-13 
KW ■ i v.l08E-I4 
Cr-L .cP (.2* .5* 1.* 1») 
WR*.t£ <7*15) 

15 FORMAT ('NUMBER OF ENDPOINTS IN THE TITRATION OF A Dlf.HOTIC ACIDM 
H >^ Is 

1 H « .5*H 

PH « -(ALOG(H) ) /2.303 

NUM « H**4 + KHHA*(H**3) + (KHHA»KHA - KHHA*MHA - KW)*(H**2) 
»-(2*KHHA*KHA*MHA + KW*KHHA)*H - KW*KriA*KHHA 

DEN « H**4 + (<HHA+MMOH)*( H**3) + (KHHA*KHA + KHHA*MMOH - KW) 
»»(H»»2) + {KHHA*KHA*Mi>10H - KW»KHHA)*H - KW*KHHA*KHA 

VMOH » -VHA*NUM/DEN 

IF (VMOH) 1*2*2 

2 CALL FPLOT (-2*VM0H*Ph; 
5 H = H - .1*H 

PH « -(ALOG(H) )/2.303 

NUM « H**4 +• KHHA*{H**3) + (KHHA»KHA - KHHA*MHA - KW)*(H**2) 
»-(2*KHHA*KHA*MHA + KW*KHHA)*H - KW*KhA*KHHA 

DEN « H**4 + (KHHA+MM0H)*{H**3) + {KHHA*KHA + KHHA*MMOH - KW) 
»*(H**2) + (KHHA*KHA*MMOH - KW»KHHA)*H - KW*KHHA*KHA 

VMOH = -VHA*NUM/DEN 

IF (50. -VMOH) 10*3*3 

3 CALL FPLOT (0*VMOH*PH) 
GO TO 5 

10 CALL YAXIS (0** 0** 1,* 15* 5) 

WRITE (7*17) 
17 FORMAT ( 'PH' ) 

CALL XAXIS (0** 0** 1*» 50* 10) 

WRITE (7*16) 
^6 FORMAT ('MILLILITERS OF SODIUM HYDROXIDE') 

CALL FINPL 

CALL EXIT 

END 

// XEQ 
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You can learn a great deal about diprotic acid titrations 
by studying and comparing and contrasting the graphs presented 
here, and by trying some different sets of conditions for your 
own calculations. You probably will find the following special 
cases helpful in reading some chemistry into the computer simu- 
lations . 

Q Case of K„ . = 4-K„.., This is a simulation of the titra- 
H2A HA 

tion of a diprotic carboxylic acid in which the two protons are 

bound by indistinguishable, non-interacting groups. 

For an algebraic discussion of the relationships between 
macroscopic and microscopic equilibrium constants for such an 
acid, and the reason for the factor U, see G. M. Fleck, Equilibria 
In Solution ^ New York: Holt, Rinehart and Winston, Inc., I966, 
pp. 91 - 93. 

A titration curve calculated with such a set of equilibrium con- 
stants is presented on page 128. This titration curve has exactly 
the same shape as the curve for a monoprotic acid present in 
twice the molar concentration as the diprotic acid. Interpreted 
as a monoprotic acid, what would be the numerical value of the 
single equilibrium constant? There is clearly just a single 
endpoint in the titration of this acid. 
(3 Case of K, 



< (H^) 



Under these 



HA" "lowest value in titration" 

conditions, only a small fraction of the second protons can be 

removed, and the diprotic acid behaves as if it were a monoprotic 
acid. The titration curve has the same shape as the correspond- 
ing curve for a monoprotic acid of the same concentration. See 
page 133 for one such case. Interpreted as a monoprotic acid, 
what would be the numerical value of the single equilibrium con- 
stant? 
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i3 Case of K„ . > 10'**K„... This large difference in equi- 
librium constant values implies that the groups that bind the 
two protons are chemically very different, or that they are 
strongly interacting* The titration curve reveals two endpoints 
if the pH at the conclusion of the titration is high enough. If 
the value of K^j^. is too small, or if M^^qjj is too small, then 
the removal of the second protons may be supressed, and the 
second endpoint may not be observed. Compare the titration curves 
presented on pages 129 — 132. 

• Fo^ one OjJ the. />lmutatlon6 that yle.td6 tm zndpoA.nt6 , ilnd 
a hzt numzfiicat vatuQ.6 (Jo^ the. ioafi microscopic equilibrium 
constants that is consistent u)ith the assumed values o^ the tu)o 
macroscopic equilibrium constants . Is your set o^ values ior the 
micro constants the only possible set? 

• Investigate the e^^ect oi concentrations on the existence 
ojj tu)o endpoints in the titration a diprotic acid. Is there 
any concentration e^^ect? Which {acid or titrant) concentration 
has the greater e^^ect? Why? 
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EXTENDING THE HALF-EQUIVALENCE METHOD TO THE 
ESTIMATION OF BOTH EQUILIBRIUM CONSTANTS FOR 
THE DISSOCIATION OF A DIPROTIC ACID 

A close look at several computer- simulated titration curves 
for monoprotic acids showed that the macroscopic half-equivalence 
method yields reliable values for equilibrium constants within 
the range 4 < pK < 10, and values that are useful as estimates 
for a successive-approximations procedure well outside that range. 
We were able to gain insight into the method, getting a feeling 
for both its usefulness and its limitations, by looking at the 
results of some numerical calculations. 

In order to interpret your titration data for a diprotic 
acid in terms of the chemical and mathematical models presented 
in this case study, you need a method for making informed guesses 
of the numerical values of two dissociation constants for that 
acid. It would be nice indeed if an extension of the half- 
equivalence method could be used. We shall explore the possibil- 
ities, testing the validity of such an extension by making some 
calculations with the IBM 1130 computer, using FORTRAN Program 
11. We shall call this extension tho, paKtlat-zqalvattnct point 
mzthod of estimating equilibrium constant values by inspection 
of a titration curve. 

USE OF A "do loop" IN A new programming feature is introduced 
A FORTRAN PROGRAM 

in this program. We have used a DO loop 
that enables us to do the same sequence of operations twice 
without punching a duplicate set of cards. Had we wanted, the 
DO loop could have been repeated many times. The statement 
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// JOB 



S71007»S51 



// FOR 

♦LIST SOURCE PROGRAM 

♦IOCS{CARD» 1132 PRINTER* TYPEWRITER* PLOTTER) 
♦ONE WORD INTEGERS 



C 
C 



TEST OF THE VALIDITY OF THE PARTIAL-EQUIVALENCE METHOD 
FOR EVALUATING THE EQUILIBRIUM CONSTANTS 



C 



DIPROTIC ACID 



REAL KHA* KHHA* KW* N IM* MHA* MMOH 

VHA = 25.00 

MHA = 0.0500 

MMOH = 0*1000 

KHHA » l.E-3 

KHA = l.E-5 

KW « l.OOSE-l'^ 

CALL PREP (.2* .55* 1** 1*) 

WRITE (7*15) 

15 FORMAT ('TEST OF THE HALF-EQUIVALENCE METHOD FOR A DIPROTIC ACID*) 

C MARKING OF THE HALF-EQUIVALENCE VOLUMES 

C MARKING PH = PKHA* AND PH ■ PKHHA 

VMOH « 1*50*VHA»MHA/MM0H 
. PH » -{ALOG(KHA) )/2.303 

DO 20 I » I* 2 

VMOH 1 = VMOH - 3. 

VM0H2 « VMOH +3* 

PHI = PH - 1. 

PH2 = PH + 1* 

CALL FPL0T(1*VM0H*PH1) 

CALL FPL0T(-2*VM0H*PH1) 

CALL FPLOT (-2*VM0H*PH2) 

CALL FPLOT (1*VM0H1*PH2) 

CALL FPL0T(-2*VM0H1*PH) 

CALL FPLOT (-1*VM0H2*PH) 

VMOH = •50»VHA»MHA/MMOH 
20 PH = -(ALOG(KHHA) )/2.303 

C CALCULATION AND PLOTTING OF THE TITRATION CURVE 



1 H « H - .1»H 

PH = -(ALOG(H) )/2.303 

NUM = H**4 + KHHA*(H»*3) + (KHHA*KHA - KHHA*MHA - KW)*(H**2) 
♦-(2»KHHA»KHA»MHA + KW*KHHA)*H - KW*KHA*KHHA 

DEN = H»»4 + {KHHA+MM0h)»(H»»3) + (KHhA*KHA + KHHA^MMOH - KW) 
»»(H*»2) + (KHriA»KHA»MMOH - KW*KHHA)»H - KW»KHHA*KHA 

VMOH = -VHA*NUM/DEN 

IF (VMOH) 1*2*2 

2 CALL FPLOT (-2*VM0H*PH) 

5 H - H - .1»H 137 

PH » -(ALOG(H) )/2.303 



H 



= 1. 
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hUt-'i = H»»/f + KHHA*(H»»3) + (KHHA*KHA - KHHA»MHA - KW)*(H**2) 
*-(2*i<;hha*kha*mha + KW»KHHA)*H - KW*KHA*KHHA 

DEN » H«*^ + (KHHA+MM0H)»(H*»3> + (KHHA*KHA + KHHA*MMOH - KW) 
»»(H*»2» + (KHHA»KHA*MMOH - KW*KHHA)*H - KW*KHHA*KHA 

VMOH « -VHA*NUM/DEN 

IF (35.-VMOH) 10»3»3 
3 CALL FPLOT (0»VMOH»PH) 

GO TO 5 

10 CALL YAXIS (0« *0« *1« ) 

WRITE (7»17) 
17 FORMAT ( •PH« ) 

CALL XAXIS (0«* 0«* 5«* 7* 7) 

WRITE (7»16) 
16 FORMAT ('MILLILITERS OF SODIUM HYDROXIDE*) 
CALL FINPL 
CALL EXIT 
END 



// XEQ 
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DO 20 I = 1, 2 



* a command to execute in sequence all the following statements 
to, and including, statement 20 — and then to return and repeat 
the entire sequence again. The integer I is a counter which is 
set initially at 1 and which increases by 1 each time through the 
DO loop. When the counter gets to 2, then the DO loop is executed 
for the final time, and the program continues past statement 20. 
A more general form of the DO statement is 



where N is a statement number (the number of the terminal state- 
ment in the DO loop) , J is an integer that gives the initial 
value of the counter, and K is another integer that specifies the 
test value of the counter, the value which I may not exceed. 



volume have been added to the sample. Ij{ all the protons from 
a , -articular functional group were removed before any of the 
protons of the other group in the acid molecule, then the point 



scopic half-equivalence point for titration of the first group. 
Such an extreme situation never is realized, but there is nothing 
to prevent us from defining these two places in a titration as 
the first and second partial -equivalence points. The relevant 
question is: "l4 tke^^e. a simple. n.Q.latlon6klv b^tWQ.Q.n thz pH valae^A 
at thz partial' ^qalvate^nae. points and thz valae.6 the. ti/oo 
e^qailib^iam aon6tant6 ^OK the. acid?" The answer: "Often, yes!" 



DO N I = J, K 



SOME TESTS OF THE 
VALIDITY OF THE METHOD 
FOR DIPROTIC ACIDS 



We shall define two partial -equivalence 



points as the points in the titration 
at which % and Vt> of the final titrant 



at which \ 



A of the titrant is added would represent the macro- 
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Inspection of the three titration curves, marked by the 
computer to show the two partial-equivalence point volumes and 
the two pH values equal to the pK values, gives us some assurance 
that this partial-equivalence method does indeed have some 
validity* On page 140 we see that for a 0.05 moiar acid with 
pK values of 5 and 7, excellent values of the two equilibrium 
constants can be obtained by direct inspection of the curve: 
that is, the two crosses have their points of intersection on 
the titration curve, at least within the limits of the accuracy 
of the plotter. With pK values of 3 and 7, the second partial- 
equivalence point cross lies on the curve, although there is 
some discrepancy with the first; perhaps pK = 3 is too small 
for the method to work precisely, although clearly this would 
be adequate for a first guess in a successive- approximations 
procedure. With pK values of 6 and 7 (page 141), the agreement 
is not so good, although again the values obtained would be 
adequate as first guesses for a procedure that could lead to 
refined values . 

Then what can we conclude? It may be dangerous to generalize 
from these three calculations, but it does appear that: 

HFor acids with certain values of the two equilibrium cons- 
tants, the two pK values are equal to the pH values obtained 
when one-fourth and three-fourths of the equivalence point volume 
of titrant have been added. 

Hlf the values of the two equilibrium constants are very 
nearly equal (differing only by a factor of 10 or 100) , then 
the method loses accuracy, but the values obtained may be ade- 
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quate for first guesses in a successive-approximations procedure 
lllf one of the equilibrium constants lies outside the range 
4 < pK < 10, then the method loses accuracy for that equilibrium 
constant. 

• Pe^jjo-^m catcutation6 dz6lQnzd to invQ.6tigatQ. the. eiJ^ect 
oi acid concentration on the validity o^ the partial- cquivalzncc 
method. 

The conctu6ion6 presented above are based on the results 
0^ just three calculations . Perform calculations designed to 
test these conclusions . What sort o^ values o^ equilibrium 
constants and concentrations uoould provide some critical tests? 

for the previously 'investigated case o^ a monoprotic 
acid, the algebra shovos that the macroscopic hal^- equivalence 
method is rigorously exact only ior the single value pK = ?• 
[G. M. Fleck, Equilibria in Solution > New York: Holt, Rinehart 
and Winston, Inc., 1966, pp. 68 - 70.] Is there any situation 
in (johich the partial- equivalence method is rigorously exact 
^or a diprotic acid? 

What sections o^ the titration curve are influenced most 
by the particular value of '^HaA^ ^HA*-^ What experimental 
data need to be most accurate for precise and accurate determina 
tion of these equilibrium constants? Design and perform some 
calculations to substantiate your statements . 
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CHANGES IN THE THREE DISTRIBUTION FRACTIONS 
DURING THE TITRATION OF A DIPROTIC ACID 

During the titration of a solution of the diprotic acid 

HjA, the relative concentrations of the three species H2A, HA', 

and A change. It is informative to examine in some detail how 

these concentrations change, because some of the differences 

between diprotic acids are revealed in the striking differences 

in the dependence of the species distribution fractions on pH 

and on V„ 

NaOH 

Equations 48, 49, and 51 can be combined in three different 
ways to yield three distribution-fraction equations: 
CHA') 1 



— (56) 

[A] (H ) K, 



+ 1 + HA 



(H2A) _ 1 

1 + ^ ^HzA^HA- ^^^^ 



(A ) _ 1 

[A] ~(hV 

+ + 1 



(58) 



FORTRAN Program 12 uses Equations 56, 57 and 58, together with 
Equation 55, to calculate the changes in the three distribution 
fractions that are predicted by our model. Three illustrative 
examples are presented on pages 149 — 151, and we shall examine 
each in turn. 

If both of the carboxylic groups on the acid H2A are 
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S71007.S51 



// FOR 

♦LIST SOURCE PROGRAM 

♦IOCS(CARD» 1132 PRINTER* TYPEWRITER* PLOTTER) 
♦ONE WORD INTEGERS 

REAL NUM* KHA* KHHA* MMOH * KW* MHA 

COMMON V* HH* KHA* KHHA* MMOH* KW* MHA* VHA 
C SUBROUTINE PL0T2 

C CALCULATION AND PLOTTING OF THE DISTRIBUTION FRACTIONS 

CALL PREP (.14* 4«* 1«* 1«} 

CALL YAXIS (0«* 0«* •!* 10* 2) 

WRITE 17*37) 
37 FORMAT ('DISTRIBUTION FRACTION') 

CALL XAXIS (0«* 0«* 1«* 50* 10) 

WRITE (7*56) 
56 FORMAT {'MILLILITERS OF SODIUM HYDROXIDE*) 

VMCH = V 

H « HH 

FHHA = H*H/{H*H KHHA*H * KHHA*KHA ) 
CALL FPLOT {-2*VM0H*FHHA) 
35 H = H - •1*H 

NUM « H**4 ■»• KHHA*(H**3) * (KHHA*KHA - KHHA*MHA - KW)*(H**2) 
«-(2»KHHA*KHA*MHA -f J<;W*KHHA)*H - i;W^KHA*KHHA 

DEN = H**4 {KHHA-»-MM0H)*(H**3) * {KHHA*KHA ■»• KHHA*MM0H - KW) 
««(H»»2) ■»• {KHHA*KHA*MMOH - KW«KHHA)*H - KW*KHHA*<HA 

VMOH = -VHA*NUM/DEN 

FHHA = H*H/(H*H * KHHA*H -f KHHA*KHA) 

IF(50«-VM0H) 30*33*33 
33 CALL FPLOT ( 0 * VMOH *FHHA ) 

GO TO 35 
30 VMOH = V 

H = HH 

FHA = KHHA*H/(H*H ■»• KHHA*H KHHA*KHA) 
CALL FPLOT (1*VM0H»FHA) 
CALL FPLOT { -2 * VMOH* FHA ) 
45 H = H - •1*H 

NUM = H»*4 KHHA*(H**3) (KHHA*KHA - KHHA»MHA - KW)*{H**2) 
«-(2*KHHA*KHA*MHA -f KW*KHHA)«H - KW*KHA«KHHA 

DEN = H**4 ■»• (KHHA-»-MM0H)*{H**3) ■»• (KHHA*KHA KHHA^MMOH - KW) 
««(H**2) ■»• {KHHA*KHA*MMOH - KW*KHHA;*H - KW*KHHA*KHA 

VMOH = -VHA*NUM/DEN 

FHA = KHHA*H/(H*H KHHA*H ■»• KHHA*KHA ) 

IF (50«-VMOH) 40*43*43 
43 CALL FPLOT (0*VMOH*FHA) 

GO TO 45 
40 VMOH ■ V 

H » HH 

FA = KHHA*KHA/{H*H ■»• KHHA*H ■•• KHHA*KHA) 
CALL FPLOT {1*VM0H*FA) 
CALL FPLOT (-2*VM0H*FA) 
55 H « H - .l^H 

NUM = H*»4 KHHA*{H**3) * {KHhA»KHA - KHHA»MHA - KW)*(H**2) 
♦-(2»KHHA*KHA*MHA KW*KHHA)»H - K^*Kttf<*KHHA -|46 
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DEN = H*»A + (KHHA+MM0H>»(H»»3) + (KHHA«KHA + KHHA«MMOH - KW ) 
♦»CH«»2) + (KHHA»KhA»MMOH - KW»KHHA)»H - KW»KHHA»KHA 
VMOH « -VHA»NUM/DEN 

FA « KHHA»KHA/(H»H ♦ KHHA»H + KHHA»KMA> 
IF (50. -VMOH) 50»53t53 
53 CALL FPLOT (OtVMOHtFA) 
60 TO 55 
• 50 CALL FINPL 

CALL EXIT — 
END 

// DUP 

♦STORE WS UA PL0T2 



// FOR 

♦IOCS(CARDt 1132 PRINTERt TYPEWRITERt PLOTTER) 
♦ LIST SOURCE PROGRAM 
•ONE WORD INTEGERS 

REAL KHAt KHHAt KWt NUMt MHAt MMOH 

COMMON Vt HHt KHAt KHHAt MMOHt KWt MHAt VHA 
C CHANGES IN THE THREE DISTRIBUTION FRACTIONS DURING A TITRATION 

VHA « 25.00 

MHA » 0.0500 

MMOH = 0.1000 

KHHA = l.E-5 

KHA » l.E-3 

KW = 1.008E-14 

CALL PREP <.l4t .25t l.t 5.5) 

WRITE C7tl5) 

15 FORMAT ('CHANGE OF THE DISTRIBUTION FRACTIONS DURING A TITRATION*) 

CALL YAXIS (0.t0.tl.tl5t5) 

WRITE (7tl7) 
17 FORMAT ( 'PH* ) 

CALL XAXIS (O.tO.tl. t50tl0) 

C CALCULATION OF THE FIRST POINT OF THE TITRATION CURVE 

H = 1. 
1 H = H - .1»H 

PH = -(ALOG(H) )/2.303 

NUM = H»»^ + KHHA«(H»»3) + (KHHA»KHA - KHHA»MHA - KW)»(H»»2) 
«-(2*KHHA»KHA»MHA + KW»KHHA)»H - KW»KHA»KHHA 

DEN = H»»<f + (KHHA+MM0H)»(H»»3) + (KHHA»KHA + KHHA»MMOH - KW) 
»»(H»»2) + (KHHA*KHA»MMOH - KW»KHHA)»H - KW»KHHA»KHA 

VMOH = -VHA»NUM/DEN 

IF (VMOH) It2t2 
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2 CALL FPLOT {-2»VM0H»Pri) 
V * VMOH 

C CALCULATION AND PLOTTING OF THE TITRATION CURVE 

HH = H 
5 H «= H - ,1»H 

Ph « -{ALOG(H) ) /2.303 

NUM ' H**U + KHHA»(H»»3) + tKHHA»KHA - KHHA»MHA - KW)*{H«»2) 
♦-(2»KHHA»KHA»MHA + KW»KHHA ) ♦H - KW»KHA»KHHA 

DEN = H»»4 + (KHHA+MM0H)»tH»»3) + (KHHA#KHA + KHHA#MMOH - KW) 
»*{H»»2) + (KHHA»KHA»MMOH - KW»KHHA)»H - KW»KHHA-»KHA 

VMOh = -VHA»NUM/DEN 

IFt50.-VM0H) I0t3t3 

3 CALL FPLOT tOtVMOHtPH) 
GO TO 5 

10 CALL FINPL 

CALL LINKtPL0T2) 
END 



// XEQ 
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identical and noninteracting, then the ratio of the two equi- 
librium constants is 1:4. All the binding sites for protons 
are equivalent, and for such a case the protons will distribute 
themselves in a random, statistically-determined arrangement on 
the binding sites* At most, only half of the molecules can 
ever exist as HA, and then a fourth are H2A, and the remaining 
fourth are A". This distribution results when there are just 
half as many protons as sites, so that each site has a probability 
of V2 of being occupied and a probability of V2 of being empty* 
The graph on page 149 shows this anticipated distribution at the 
macroscopic half -equivalence point, where the distribution frac- 
tion (HA'')/[A] reaches its maximum value of V2 , and the other 
two distribution fractions are equal to Vi^. 

• V^^^o^m 4Dme catcutatlon^ to Investigate, tko. 6hapz4^ oi 
the. dlstfilbiitloYi' {^Kactlon cuKv^s {^ok othzK hztM oi tu)o zqiiltlbfilam 
con6t(iYit6 that 6tand in the. Katlo U4. Thz pfizdlatlon, bah ad 
on 6tati6ticat conhldzKatlonh ^ that a Katlo o^ J:2:l maht 
K<L6att vohzn thzKZ afiz just hat^ e^noagh pKotonh to occupy alt the 
binding hltu. VoeJ^ thlh condition necessarily occur at the 
macroscopic hal^- equivalence point? Can you ^ind a situation 
in u)hich the fraction (HA')/[A] does not pass throuo- a maximum 



In order to achieve a solution in which more than half of 
the A-containing species are present as HA , we must chose an 
acid that has equilibrium constants that differ by more than a 
factor of four* Even when the two equilibrium constants differ 
by a factor of 100 (as is the case depicted on page 150), there 



oi V2? 
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remain significant amounts of species H2A and A" at the point 
in the titration where the fraction present as HA passes through 
its maximum value. 

• ln6pzctlon 0^ the. g/iaph6 on pagQ.6 149 — 151 6uggQ.6t6 that 
the. ^Kactlon6 (H2A)/[A] and (A'')/[A] are equal M;feeneve^ the. inac- 
tion (HA ) /[A] machz6 It^ maximum va£ae. A^e thzm any con- 
dltlon6 iindzn vohlch thl6 may not fae tnuo.? 

• Undtn u)hat condition^ dan 6olutlon6 fae pxapafizd In t^hlch 
two 6pdcld6 ha\)t dl6tfiibiitlon ^fiadtlonh zaah ^qaat to V2? Cfeecfe 
yoa^ an6u)^^ i/olth appfiopfiiatz hlmatatlonh . 

Two endpoints are seen in the titration curve presented on 
page 151 for a diprotic acid whose equilibrium constants differ 
by a factor of 10**. The condition for a sharp endpoint is 
generally that a few drops of titrant must change the qualitative 
nature of the chemical solution, and such is the case in this 
titration. In the vicinity of 12.5 m£ of added sodium hydroxide 
solution, the solution changes from a mixture of H2A and HA to 
a mixture of HA^ and A". There is essentially no place in the 
titration where all three species coexist. At 12.5 mi, the frac- 
tion (HA'')/[A] is essentially equal to 1.0, and we can say that 
at this point all the H2A has been titrated to HA" ; the corres- 
ponding statement is not valid for the acids whose titrations 
are simulated on pages 149 and 150. Likewise, we can speak of 
a solution of HA* for this acid, but not for the other cases, 

• ItifiltiL and ^x^cute. a compute.^ pn,ogn.am that glve.6 a nume.M'- 
Icat output 0(5 dl6t^lbutlon'' inaction vatuQ.6 during a titration* 
By hou) much do the. zqultlb^lum constant vatuQ.6 have to dl^c^ 
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In oKdo^fi that tho. inaction (HA')/ [A] can Kzach the. value. 0.9? 
the. value. 0.95? the. value. 0.99? the. value. 0. 999? l6 the. limit 
e.qual to zxattly 1, ok afie. thtfit KZ6tKlctlon6 that feeep the. 
maximum attainable value, o^ thl6 IntzKmzdlatt 6pz(ilz6 conczntJta- 
tlon le66 than 1? 

• Which oi the. vaKlou6 hypothetical ca6Z6 o^ dlpKotlc acld6 
pn,e.6tntzd fee^e coKfit^pond to fizal acld6? What, In iact, aKZ the. 
fiatlo6 OjJ zqulllb^lum con6tant6 that a^e zncounte.Ke.d In thz lab- 
o^ato^y (A)lth common dlpKotlc ca^boxyllc acld6? othe.^ dlpfiotlc 
acld6? Look up 6ome. zqulllbfilum- constant value.6 In thz llbKaKy 
to ge.t a izzllnq jjo^ thz ^angz o^ tqulllbfilum- constant ^jalue.6 
that can fae Kza6onably expected. 



-1 ' yO 
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CHANGES IN THE FOUR MICROSCOPIC-SPECIES DISTRIBUTION 
FRACTIONS DURING THE TITRATION OF A DIPROTIC ACID 

When talking about the equilibria involving a diprotic 
acid, it is often helpful to consider indepenuently the two 
different species that have the overall composition and charge 
indicated by the symbol HA'. Here we explicitly recognize the 
two species H^A' and 'A^H which differ because a different 
proton binding site is occupied. The; e is only one species with 
the formula H2A, and that is H^A^H. There is also only one 
species devoid of dissociable protons, and that is A"". 

The fully-protonated species can react to form the com- 
pletely-dissociated species by either of two distinct reaction 
pathways which differ in the order in which the two protons are 
removed. The four chemical steps and the associated mlc^06Coplc 
zqalllbKlim con6tant6 in this model for diprotic-acid equilibria 
are 

' + . (K*)<'A--H> 

H'^A^H «♦ H + A^H ""i = (59) 

(H^A^U) 

♦ (H*)<H-A") 
Vi^A^n ii H + n^A K2 = (60) 

(h*)<a"> 

H^A i* H + A" K3 = (61) 

nH-'A') 

(h*)<a"> 

A^W ♦» H + A" U = (62) 

("A'^H) 

where the individual st>ecies are written with indication of the 
position of a bound nroton: left and risht sides of the letter 
A indicate different Dositions in tne mf^lecule. The molar con- 
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centrations of these individual species (we shall call these 
concentrations micKo6i:,oplC'6pQ,clQ.6 concQ.ntKation6 whenever there 
is a need to make a careful distinction) are indicated by en- 
closing the formula by the symbols < >. In this section we shall 
indicate microscopic equilibrium constants by using subscript 
numbers, and macroscopic equilibrium constants by using sub- 
script formulas. 

The four microscopic equilibrium constants Ki , K2 , K3, and 
Kif are not independent. If the numerical values of any three of 
these K*s are known, the value of the lourth can be calculated 
without difficulty. 

• Comblm Equation^ 59 — 62 to obtain thz Kzlatlon^hlp 

YiiU = K2K3 (63) 
However, there are three independent microscopic equilibrium 
constants, even thoufyh there are only two mrcroscopic equili- 
brium constants whose values can be obtaine^l by analysis of 
titraj^ijjn ^d^tV. This fact implies that knowledge of>^the values 
of th^^ ^KpVrimentally-determined macroscopic equilibrium cons tan ts| 

is not sufficient to specify a unique set) of valued for the four 

*i . !*. 
microsoibpic equilibrium constants. We can, howpvei'^.. obtain some 

' W ' 'w % ' 

relatiiOU^Hips between the mafcro and micro const^f^ts and we shall 
proceed with that task at this time. ^^ 

The link between the macro and micro descriptions of these 
equilibria is the set of equations 

(H2A) = (H-^A-^H) , (64) 

(HA") = <'A-^H> + iH-^A') (65) 

(a') = <a"> (66) 



• Combine. Equations 48, 49, 59 — 62, and 64 — 66. to get 
the. following two Kzlationships batvoaan thz ma(in,o6copic zqai- 
t^bn^iam constants and the micn,o6copic e,qa<libn<um constants : 

^HaA = + K2 (67) 

To obtain algebraic equations relating distribution frac- 
tions for the four microscopic species to (H*) , we need the 
conservation equation 

[A] = (H^A^H) + CA^H) + <H-^A'> + <a"> (69) 
Equations 59 — 62 can be used to eliminate all but one of the 
microscopic species concentrations in Equation 69. Depending 
on which concentration remains, we can get any one of four 
distinct equations which can in turn be rearranged into a con- 
sistent form so that the four distribution fractions can be 
written as 

<H^A^H> (h'*')^ 



[A, '^•^ 


dhnominaioK 




, ^B*)K2 \ 


[A] 


dcnominaton. 

V 


(h^a') 


— i 


[A] 


denominator 


<A=> 


K2K3 


lA] 


de^nominatofi 



(70) 



\ 



1) 

(72) 
(73) 



where 



dznamlnatoK = (H*)^ + {Ki + KzlCH*) + K2K3 
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• Pe^jjo^m tho, atgzbn,a n,2,quln,2,d to obtain Equatlon6 70 — 
73* Vou may ^Ind It donv tniznt ofi nk.(LZ66aKy to U6Q. Equation 
63 to slmpZi^y dZfitaln 2.xpK2.66lon6 . 

• Shou) that Equations 56 — 58 can be vo^lttcn a6 thfiQ^a 
inactions volth a common denominator.. 

• Starting volth Equations 70 — 73, and U6lng relationships 
bctijoccn microscopic and macroscopic quantities, obtain Equations 
56 — 5S. 

• Verify that thz f^our microscopic distribution fractions 
together add to unity [that Is, that the sum ol the iour fractions 
Is equal to ^) . _ 

COMPUTER CALCULATIONS OF THE FORTRAN Program 13 provides a 
MICROSCOPIC DISTRIBUTION 

FRACTIONS DURING A SIMULATED convenient way to learn about 
DIPROTIC ACID TITRATION 

changes in the four microscopic 
distribution fractions during the course of a titration. The 
numerical values of three microscopic equilibrium constants are 
assumed, as well as a set of initial condi t ions \ (vol lime of p^d 
solution and concentrations of acid and base) . For ^tSie inf^raja- 
tion of the chemist who is using this program to simulate varxpus 
titrations, the^j is a printout of pK values for both the micro- 
scopic and macroscopic equilibrium constants. Then the program 
produces a huge array of numbers, a sampling of which is pre-^ 
sented on page 160. Actually, you probably cannot find any real 
acid that behaves like this simulation, since K3 was assumed to 
have a larger value than either Ki or K2 . Generally it is hardeV 
to remove the second proton, which means that the second equilib- 
rium constant is smaller. 




// JOB 571007^551 
/ / t .}i< 

*IwC.b(CArwf 113? PIUMLX) 

REAL .^UXf KHAf KHHAf KWt MMQhi t MHA t 1 1 K2 1 K3f 
C CALCULAflO^ OF FKACTIOMG OF TOTAL A-CO.\T A I \ L;G SPiCILS PHLSlihT 

C AS EACH Jr THE MICROSCOPIC SPECIES DUKIr>G THM 

C TITRATIO^i OF ThE ACID H2A 

.<1 = l»00t-4 

K2 = 3«0d-4 

K3 = 9.0i:-4 
= K2»<3/K1 

<YiHA = Kl + K2 

.<HA = <3*K4/ (K3 + K4) 

J'<i = -(ALOG(i^l) )/2«303 

P<2 = -{ALaG(K2 ) )/2#303 

P.<3 - -(AL0o(K.3 ) )/2«303 

P<,4 = -(AL0i.i(K3 ) ) /2«303 

PKHHA = -( ALOG(;<h!HA) ) /2«303 

^^.<hA = -( ALOO(KHA) )/2#303 

aHITL {3tll) PKlt PK.2t PK3f P^U 
11 HOKy.ATC'lPKl >='f K6#2f • PK2 K6#2f • f^K3 f F6#2f 

^•PK.4 , F6#2) 

//HUE {3f 14} PKHMA» PK.HA 

14 FCR*-1AT CIPKHHA =«» F6#2» • PK.HA =»f F6#2//) 
WRITE (3f 15) 

15 FOK.vaT (» CH2A) CHA) (AH) (A) VMOH PhM 
VHA = 25«00 

MHA = 0^0500 
'•'iv'.OH = C»10CC 
.<W = 1.008E-14 



il.. _ _ » Al .^'1/ lH 1 / T 0.-%0 i 



= -(ALOC^Un ) /2^303 
1 2^K<^r^;*i^*^/>r';nA + N.\»KHMA)*h - /;A*iO-tA*J<;HMA 

\ = ^^'*V4*^ UMHA + >ivl|)H)»(H»*3) + Vc;<HHA»rHHA + .^HHA^XMOH • i 

|; ♦*'tf^^*^2J ;4iMA*,<MA*i\l-iO'\ - K.»v».<.HhA)*M- K.vV*K.HA*KHriA 

\ IF (Vi'^OH)' if'^ZL'f 20 f\ 



i 

i2 IF {Ph - 7#i36) If It 30 • ' ^ 

2%} CONTIMJL t \ I ; 

^1 FDcis ^ H^^2 + (iCl + K2)^H + K3*J<;? 1' 
F^iAH = (h»*2 ) /FL;t.r, ^ \ 



FH/k = H»<2/FDEr; : 
FAN = H*K1/FDEN ' 
FA = <2»<3/FDEK 
i>M = -( ALOGC.H) )/2#303 

v%kI?c (3fl6) FHAHf FHAt FAHt FAt V^'Oh » PH 
16 FU^v.AT (6F10.6) 

^ GO TO i 
30 CALL tXIT 

END 

// XEO ,59 



ERIC 





(H2A) 




(HA) 




(AH) 


(A3 


V>'OH i-M 




Vll7^^ 






u • 


J 1 V 'i . - 


0« 




w.w:; .'9Vv 


aI . 3 ^ ^ /' 2 






Oi 


070699 


Oi 


C23t>33 


C« 




o.6^.e{-.9:; 


2.4r..315 


0« 




Ui 


V»V^ yO 


1 


0/d476 


Oi 




i.2o599!- 


V' 2 J 7 


C- 


-'•-ryi yj 


i 


lJ24v:6 






Oi 


J3dll« 




.1 . 1 i 


()• 




Oi 


121267 


0- 


040^22 


Oi 


036^24 


2.fi922^4l 


2.712993 


0- 




v.» 1 


I4l3/i? 




::47i27 


Oi 


0^22 30 




2 ..:.;v3;>t? 


0« 


6669 /6 


Ji 


101634 


Ui 


Oi>3b78 


Oi 


11 7b 11 


!;#2t^71i9 


2.906778 


0< 


59t)i)26 


Ji 


1B0399 


Ui 


060133 


Oi 


1639^1 


6.790356 


3.C0 367U 


0« 


t) 16792 


Ui 


19p6b6 


Oi 


06t?/2b 


'J i 




v:5ic9l5 


3.100563 


0- 




Ji 


20t?306 


Ui 


06B5J2 


w i 


29li;-j9 


10.495376 


3.197455 


Oi 


3b22'ib 


Ui 


206399 


')» 


U69a66 


J. 


36 7i-US 


12.530477 


3.29434fc 


0- 


27571i? 


Oi 


► 2039UO 


Oi 


067969 


0. 


432^07 


14.543225 


3.391241 


0« 


20b^8tJ 


Ji 


►192737 


Oi 


064245 


0. 


33^323 


16.441021 


3.488133 


Ot 


► 1527^*7 


v; i 


► 176i>09 


Ui 


058636 


0. 


611903 


iG. 127510 


3.5e5026 


Oi 


► loebf>9 


U 1 


»1572b7 


Oi 


►052419 


0. 


.681^53 


19.564746 


3.681919 


Oi 


► 075H2;; 


-Ji 


»136913 


Oi 


►045637 


0. 


.7^f 1620 


20.746391 


3.778811 


Oi 


»051b37 


Oi 


►116996 


Oi 


►038993 


0. 


.792167 


21.692054 


3.875704 


Oi 


►034917 


Ui 


►09HSUC 


Oi 


►032b36 


0. 


.S'5373S 


22.434622 


3.972596 


Oi 


»0232f>Q 


Oi 


•0B1991 


Ci 


►027330 


0. 


.c67'^2-j 


23.0113C6 


4.C69489 


Oi 


»U1534^ 


Oi 


i06763t5 


Oi 


i022546 


0. 


.fc9'»^72 


23.456131 


'♦.166331 


Oi 


i'JiaCf>7 


Ji 




Oi 


►018472 


0. 


.91605^ 


23.796408 


4.263275 


Oi 


1006^57 




.045163 


0< 


.01:>'j54 


0. 


.93322'* 


:'4.g6l9xil 


4.360167 


Oi 


*00U2bl 


Oi 


.'.)366t>9 


Ot 


.012219 


0. 


•9^6863 


24.265407 


4.457059 


Oi 


•002756 


c 


.029662 


Ot 


.009bd7 


0, 


.V37693 


24.422V43 


4.553952 


Oi 


.001779 


J 


.023943 


Ot 


.007981 


0, 


.966293 


24.5454C6 


4.650845 


Oi 


.O011'V7 


0 


.C19290 


Ot 


.006430 


0. 


.973132 


24.64096b 


4.747737 


Qi 


• 0(. J7'JS 




.Olbbie 


Ot 


.005172 


u, 


.978570 


24.71581? 


4.84463U 


Oi 


• UO(V*?'V 


\» 


.012469 


Ct 


.00'^ 156 


0, 


.932S99 


24.774627 


4.941522 


Oi 






. J 1 L J 1 J 


t 


.0^3336 


u. 


.9663^7 


24.820991 


5.u384l5 


Oi 


»U0ul9t) 


u 


. jJiiC3'J 


Li 


. wv?2676 


0 


.9c9wi96 


24.b^>7624 


5.13530O 


Oi 


»0 jC12t) 


u 


.006438 


Ot 


.002146 


0. 


.9912«9 


24.806C34 


5.232200 


Oi 






.003160 


Ot 


.001720 


0. 


.993039 


24.9CS656 


5.327092 


Oi 


»000051 




.3J4134 


Ct 


.001376 


0, 


.99^^36 


24.927948 


5.425966 


Oi 


iUOA)U33 


Ut 


.0033x0 


Oi 


.00lij3 


0, 


.993332 


24.942504 . 


.5*52287e 


Oi 


►00UC?i 


J i 


.002651 


Ot 


.00u^b3 


.0. 


.996^^'4 


24.954097 \ 


iJ. 619771 


Oi 


*O0U\jl3 
» JUUU'Jd 




.UJ2122 


Oi 


.Ouw7o7 


0, 


.997136 


24.963333 


5.716663 


Oi 


u 


.:'Mjl69b 


U t 


.s;Uu566 


0 


.997726 


24.970710 


5.813556 


Oi 






.Jul3b9 


Ot 


.000453 


0, 


.99clHl 


24.976589 


5.91044-9 


Oi 








0 


.u0u362 


u, 


.996345 


24.981288 


6.0^-7341 


s)i 


1 U'JUwv^2 


J 






. 0OU^9O 


u 


.99Hfi36 


24.9B5U42 


6.104233 


Oi 


1 UwU'J Jl 


'J 


.0'JU6*^6 


Ot 


.U0O232 


u< 


.99VJ69 


24.988044 


6.2Jil27 




i y'^-O ..^ 


- 


• . j'JS'v 7 




.wL.. Ir.O 


0 


.999253 


24y VVC.4-*'* 


-5^298019 








• J^i.^'*6 


c 




0 


.999/t04 


2^*.9^2 3i>6 


o. 394911 








• .;.y v}356 


0 


. C*; 0 11:; 


c 


.999523 


2^. Vv389.if 


6. ^91804 


C-i 




• I 'w^2Bt; 


J 


.CU*^. v5 


0 


.999619 


24.9Vi;l/:5 


i^r 5::b697 




• •)' ' ^ 


- 


• '^ ^C2/^e 


f 


» J -J . ' ^*^'^ 


0 


.999593 


24.V9;:1U 


:;.C*c:55dV 


0 






• ' Jwlb2 






0 


.99^736 


VVt: VUO 


u.7c24d2 


c 






.'Jw«Ul46 


i. 


»0U^^4:'; 


u 


.999804 


24.9v7535 


6.879374 


Oi 






.:;.;0il6 


- \ 




0 


.999a/+4 


2^. V ^0 w'* w- 


o .97t>2o7 


Oi 






• C JJt:93 




• COs.051 


c 


.999373 


24.99^462 


7.C73160 


J< 


» V J w ^ ^ 




• jj.i74 


w 


• w w w \/ ^* 


0 


.9999J'^ 


2H.9Vb7vr 


7.170052 


\' i 


• 00 w w </> sJ 






■.' 


..;JJ.;19 


0 


.999920 


24.v:?vC6'^ 


7.266045 


Gi 






,47 


c 




0 


.999936 


24.9v929r. 


7.36383c 


J 






• ^ w 3 B 


's/ 


• ^Js.^12 


0 


.99v94e 


24. Vvv4:;t: 


7.4o0730 


0 




.J 


• 0.iU*.3 J 


J 


• J'wwt^l J 


0 


.999939 


;'4. V996P6 


7.557622 


U. 






.JX.U..24 




• w u W ^ ^ O 


u 


.999967 


24. vvvb05 


7.654516 


0 






• J-j19 






u 


.9999 74 


2^.9v995j 


7. 7^1408 


r * 

V i 


• U W « *^ 




• ^ w w ^ 15 




» w u V 5 


0 


.9999 79 


2 5 . w w u r 1 


7.84830U 


^) 


• Jw..w 




• J^v.;12 








.999983 


25.0^w2H/ 


7.945193 


• Ov'V/U«« «^ 






■Jf 


» V u w ^ w 3 


y 


.9999,^6 


/5.Uv'J3b5 


5.042087 










c 


i u u J 


u 


.999959 


;;t>.Ov05t>6 


8.i3l'.976 



160 



This program calculates pairs {Vjj^^qjj ,pH} for a tit.ration, 
and then uses the pH value (actually, it uses the value of the 
concentration of H*) together with Equations 70 — 73 to calculate 
values of FHAH, FHA, FAH and FA, the four microscopic distribution 
fractions, at each pH value. Statement 20 has not been used in 
previous programs; it doesn't do anything, but it seemed to be 
a good place to branch for an IF statement. The CONTINUE state- 
ment could have been eliminated by numbering the following state- 
ment as Ms 20. 

One significant feature of the data from this -Jimulation 
should be noted. The ratio <H^A' >/<"A'^K> is constant throughout 
the titration. Since this ratio does not depend on pH, there is 
no information about the ratio that can be obtained from the 
shape of a titration curve. There is in fact information about 
just two equilibrium constants in the titration data. This fact 
means that analysis of a titration curve yields insufficient 
data for evaluation of the microscopic equilibrium constants. 

Chemists can and have evaluated the microscopic eq{^B^T±um 



constants for diprotic acids, often 'by very clever means. For 
a discussion of the evaluation of micro constants for amino 
acids, see J. T. Edsall, **Dipolar Ions and Acid-Base Equilibria," 
chap, h in E. J. Cohn and J. T. Edsall, Proteins ^ Amino Acids and 
Peptides (American Chemical Society Monograph th 90), Hew York: 



Reinhold Publishing corp., 19^3. 
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CONFRONTATION OF EXPERIMENTAL DATA WITH THE 

PREDICTIONS OF A DIPROTIC-ACID MODEL: 

A STRATEGY FOR EVALUATING TWO EQUILIBRIUM CONSTANTS 

The critical test of a chemical model is a confrontation of 
the predictions of that model with data from appropriately de- 
signed experiments* Such a confrontation simultaneously tests 
both the adequacy of the model to interpret the data, and the 
adequacy of the data to test the model ♦ With our model, we use 
the process of confronting data and model to find the set of two 
macroscopic equilibrium constants that bring data and model into 
the best quantitative agreements The result is that we can 
obtain numerical values for the two equilibrium constants* We 
are attempting to obtain a consistent description of the chem- 
istry of a diprotic acid in solution* 

The shape of the titration curve is most sensitive to the 
numerical value of K^^^ in the vicinity of the first partial- 
equivalence point (that is, ^near.the point where the volume of 
added titrant is % of the vplum^ i-e^uired to reach * the jequiva-* 
lence point of the titration!. T^iese are the data thatjare most 



useful in 



evaluating that equilibriui| constant. The nuirkrical 
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predictions of the model are most sensitive to the numerical 

4 i ■ 

< « 

value of Kjj^^ in the vicinity of the secbnd pjjrtial-equivljlence 
point (near the point where % of the equivalence-point volume 
of titrant has been added) , and these data are "^he most ustiwful 
in evaluating that equilibrium constant. Only a dafiziiilty- 
^zttdtzd £Kom youn. a.Kpa.^imQ.ntat titration data a^e 

needc.rf jJo/L the. 6ucct66lvQ.'app^oxlmation6 p/toceda/te that to.adh to 
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// FO.v 

♦ LIST SOUkCh PKOGRA.v 
*ICCb{CAKL>. 1132 PXIaTl:-.) 

C COuFRONTATlON OF ThE DATA xITH PREDICT I'JNS Fu-c A JlPi-cOTiC AGIO 

i<i£AL <»'A» <fJ'liA» KV. » NJ --.» MHAt ^■■i'.OU 

VHA » 25.00 

•iHA «" 0,05 
>!MOH = 0.10 

KHHA « •9E-6 

" i^HA i"l .Ttf-7 " 
K.W = l.i-:-14 

28 JieAD_ C2»2j>) V» PH 

25 FORMAT rF2b.2» F20.2) 
IF (V) 30t29»29 

29 H » 10**ij!!PH' 

" KUM «H»«4 + 'khMA«(H«»3) ♦ ^.<HHA«KHA^ - Ki)hA«--;HA - K.'a')« (M*»? ) 
*-( 2*KhmA»K»HA*MHA + A.W«.'<.HHA )«H - !(.W*KHA*K.HhA 
OEH^ H««4^ +_ ( KHHA + Mp^Om ) ♦ { H««3 ) ♦ (i^nHA«iCHA + KnHA*TAXOH - (^.v. ) 

■ ♦«TH**2T + < l<;:-^HA^K;MA*MyiOri - MV«K.HHA)>H - X,W*K.hA*K,riHA 

VMOH « -VHA*NUM /DEN ' - - . . . _ 

WRITt (3»10) Ph» VMOh. V 
1? FORMAT (3F7.2) 

^— CD"T0"'28"~' - 

30 CAUL EXIT 
END 

-/7~XEy— - • ■ ■ 

5»03 1. 14 1.2t 

^32 2^09' 2.2li ' ■ ■ ' 

5.71 4,36 4,65 
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namzfiJidat vaJLaz^ io^ the. ti^o zquA.libn.lum con6tant6 ioK the. acUd. 
Then the full set of data from throughout the titration can be 
compared with the predicted titration curve for the full titra- 
tion, using these "best values** for the macroscopic dissociation 
constants . 

FORTRAN Program 14 is a simple and convenient way to compare 
a few data valucs with the numbers that the model predicts, for 
an assumed pair of equilibrium-constant values- This program 
reads a card on which are punched paired values of pH and titrant 
volume. The pH value is converted into a value of (H*) , and 
Equation 55 is used to calculate a predicted value of V^^^q^- 
The experimental and predicted values of volume are then printed, 
together with the corresponding pH v^aiue. Several sets of cal- 
culations, each with a slightly different value of one or both 
equilibrium constant(s), can be performed until you feel that 
the best agreement has been achieved. 

• 1)z6ign a 6zt catcutatlon^ to vo^^iy {on m^ute.] the. 

athzfvitlon that ct^T.taln no^gloru o£ the. titnantlon ca^ve am tha 

• > * 1 'if 

^ r^.o^t 6Q.mitivQ. mdi(^a6 uT^h nl-iS^eci to vanXationA in the warned- 

icai vatuz6 oi ce^tta^w zqkitihfix^um ^on6tant6 . 

• an ^id t6^ <jkwd^ng|tfie ^et calcutation^ that givu the. 
bz6t ag^eemew^l voni^te. ^ p^oglam||tfeai compam6 the calculated and 
pncdictzd titnant voZumc6, pznhap6 by taking the di^^cncncc, on 
the natio, o^ the6e tvoo numben6. Ji you voant mone significant 
figunes in the output, what should you change in the pnognam? 

• Changing the value of K^^^ has some effect on the pnedic- 
tions of the model in the vicinity of the finst p an tial- equiva- 
lence point. Design a set of calculations {and ptnfonm the 
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// FOK 

♦ LIbT SOUKCb PkOC-.kA'V. 

♦I0C6(CAKDf 1132 PKlNTcKf TYPtwKITtKf PLOTTLK) 
C TITRATION CUKVh 

C CONPAKISOa CF LXPhKIML/HAL DATA WITH THL PKlJICTIOisS OF A 

C CHt:/'»ICAL .XODLL 

C DIHKCTIC AClJf ASSU.^'ED VALUES OF «CHA AM^ <^^HA 

C DATA CARDS hAVL Vr-Xrt PUNCHED v^ITHI^ THt FIkST 20 SPACES 

C DATA CARDS HAVE PH k^Ui\CmLD A'lThIN TmE btCO:<D 2C SPACES 

C IMCLUJE DECIMAL POIST Ii\ EACH VALUE OF Vi>;OH AM; PH 

C LAST l^ATA CAKU MUST HAVE A .NEGATIVE VALUE OF V.vOH 

KEAu KHA» KHHAf KWf liUt^^f '^*HAf f-MOH 
K.Hr'.A = l#lE-6 
KhA = 0#9E-7 
VHA = 25#C0 
.v:ha = 0#0&00 
MMOh = O.iOOO 
<%v = i#OOI5L-l4 
CALL PKEP (•2» .55f !•» 1.) 
WKiTE (7fl5} KHA» KHHA 
15 FOkMAT (•KHhA = ••E8«2»» <HA = S£e#2) 

C CALCULATION AND PLOTTI.'^iG OF THE TITRATION CURVE 

H = 1# 

1 H = •5»H 

PH = -(ALOGCH) )/2#303 

N\J>\ s M*«4 + K.HHA* ( H«<^") + {:<HHA*KhA - <nHA*MriA - ^V;)«{H««2) 
*-( 2»KHHA»KHA*MhA + Kw«KHHA)*M - KA«KhA«KnHA 

DEiN = H«»4 + (KHHA + A'MOH) ♦ {h«*3 ) + (KLHHA«<HA + f;HHA*yyCh' - <W ) 
♦ ♦(H**2) + (KLriHA«KHA«.MiV.OH - <W«KHMA)*h - JCW *^^iA*^C^tHA 

V.*'OH = -VHA«(^;UK/DE»\ 

IF (Vr-IOH) If2f2 

2 CALL FPLOT (-2fVM0HfPH) 
5 H = H - •1*H 

P.M s -(AL0G{H) )/2.303 

NUM = 11**4 + KHHA*(fi**3) + (KMHAKj^HA • <hHA*V.nA - n/;)*{H**2) 

2»KHHA*KHA*MfiA + iCW*?;HHA)*H - <W*<HA*ICHMA 
OEU - M**4 + (KnHA + rvyOrj)*(H**3) + (MiHA*i^HA + .<MHA*MyuH - KW ) 
nai-i^if^) + MA^Jrs'-./^'i'M.sOM ~ .<.A'*KH.iA)»lt - .^V.-^^s'iA^r.nfiA 

Ir i-jC.-V'-.jfn 10*3*3 
i CALL FPLOT ( 0 ♦VMOH»PH) 

(iO TC 5 
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C PLOTTIf.o DATA l'0I.,7S F\C.v THt TITUATI.. . 

10 COUTIriUL 

28 READ {2t25>) VtP 

2b FORMAT (F20.2.F20.2) 
IF (V) 30«29«;?9 

29 CALL FHLOT Jl.V.P) 
CALL FHLOT {-?»VtP) 
CALL POI.<T (1) 

00 TO 28 

C DRAWING AND LAHtLLIiNG THc AXES • 

30 CALL YAXIS (0.» 0.» l*t 15» 5) 
WRITE (7»i7) 

17 FORMAT { 'Ph* ) 

CALL XAXIS (0.*0.tl*>!jO*10) 

WRITt (7.16) 
16 FORMAT ( 'MILL ILI TENS OF SODIUM HYDROXIDE') 

CALL FINPL 

CALL EXIT 



// XEJ 



ERIC 



^9 



166 



tWA = 0-90E-07 mA-0-ll£-05 
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<iatcutation6) to a4>'tZ66 thz Q^Ho^tt oi int2.Kation6 bo^tijozzn the. 
tu)o zquitib^iam con6tant6 on the. /^hapz o^ the. titration cu^ve 

^eg^ow^ that a^e critical ^o^ the. zvatuation o^ tho/iZ equi- 
tib^ium coMtant/i . The. etJ^Jeci (jooatd be expeciecf to be di^^zKznt^ 
d^p^nding on thz ^ztotivz vatu^^ {and on thz ab^otutz valuer] o^ 
tha zqaitib^iam aon6tant6 . 

A GRAPHICAL CONFRONTATION FORTRAN Program 15 allows you to 
OF THE DATA AND THE MODEL 

compare data that you have punched 
on cards with predictions of the model in the form of data points 
placed on a calculated curve. Most people find the visual impact 
and effectiveness of this graphical presentation much greater 
than with a listing of numbers. 

All the data used in the calculations whose output appears 
on page 167 were simulated, using the same algebraic equations 
but slightly different equilibrium-constant values. The data 
used in FORTRAN Program 14 were also obtained from a similar 
simulation, but with a third set of equilibrium-cons tant values. 
The |i umbers used are as follows: 
♦ \ SIMULATED DATA PROGRAM 1^ PROGRAM 15 

1.0 X- 10"^^ 0.9 X 10"^ 1.1 X 10-^ 

1.0 X 10"^ 1.1 X 10-^ 0.9 X 10^^ 

The difference between data and predictions is just barely de- 
tectable on the graph, whereas the numerical output shows the 
difference much more clearly. The overall effect can be grasped 
more readily by looking at the graph. Probably both methods can 
be used to advantage. 
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ASSESSING THE PRECISION An experimental scientist must always 
OF EQUILIBRIUM-CONSTANT 

VALUES OBTAINED BY USING be concerned with the reliability and 
PROGRAMS 14 AND 15 

potential repeatability of his or her 
experiments and calculations. When a number is reported, he or she 
should indicate (by the number of significant figures in that 
number, or by a statement of the precision -^f that number) the 
probability that a repeat of the work by an independent investi- 
gator would yield the same number. 

Thus it is important to see if two different equilibrium- 
constant values give equally-good fit of the data, within the 
limits of experimental uncertainties of the data. After you 
have found your '*best value", you should repeat the calculation 
with ^ slightly altered value, checking to see if any significant 
variation occurs. If there is no discernible difference in the 
agreement obtained in using 1.0 x 10"*^ and 1.1 x 10"*^, then it 
is misleading (and therefore not good science) to report a value 
as 1.00 X 10"*^, since the two zeros imply a knowledge of the 
digits that belong in those two decimal positions. 

How can you estimate the '^limits of experimental uncertainties 
of the data*'? This is not always easy, but it is always important. 
For some guidance in assessing the limitations of burettes and 
pH meters, read over the information in Appendix II. Also keep 
asking yourself the question: **If I repeated the work, would I 
get the same answer a second (or third, or fourth) time?** If you 
don*t know, then repeat some of the work and find out! 
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CAN A CLEARLY- INCORRECT MODEL 
GIVE A GOOD DESCRIPTION 
OF AN. ACID-BASE TITRATION? 

Throughout this case study, we have used a simple and straight- 
forward criterion for testir^g a chemical or mathematical model: 
"Do thz predictions OjJ the. modzl agA.ee, uoithin the. limits OjJ 
e.xpz^me.ntat unce.n.tainty , taith zxpzkime.ntat data?*^ When there 
fails to be adequate agreement, this criterion signals the chemist 
that problems with either the model or the data (or both) exist. 
But what can be concluded when there is excellent agreement? 

We shall investigate a clear case in which data are inter- 
preted in terms of an incorrect model. In FORTRAN Program 16, 
we simulate data from the titration of a diprotic acid, marking 
these points with X marks on a piece of graph paper. Then, using 
the same axes, we plot as a smooth curve the predictions of a 
monopro tic- acid model. The results are shown on page 173. The 
data points lie precisely on the curve. 

There is no way to tell, by analysis of titration data, 

whether this acid is a monoprotic acid with K^^^ = 0.5 x 10"**, 

or a diprotic acid (present in half the molar concentration) 

with K„ = 1.0 X iQ-^ and K„.- = 0.25 x lO"'*. The data are 
H2A HA 

interpreted equally well with either model. 

This is a special case, in which the diprotic acid equilib- 
rium constants were chosen to stand in the ratio 1:4. But if 
we vere to introduce some experimental uncertainty (some random 
scatter) into the data points, we might be able to apply an in- 
appropriate model to a rather wide range of data. You can 
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// JOB" S71007#S51 



// FOR 

• LIST SOURCE PROGRAM 

•lOCSCCARO* 1132 PRINTER* TYPEWRITER* PLOTTER) 
•ONE WORD INTEGERS 

C OIPROTIC ACID DATA INTERPRETED VIA A MONOPROTiC MODEL 

REAL KHA* KHHA* KW* NUM* MHA* MMOH 
VHA > 25 .00 

MHA » O.O&OO 

MMOH > 0*1000 

KHHA • l«E-4 

KHA ■ •25E-4 

KW > 1«008E-14 

CALL PREP {•2* .SS* 1.) 

WRITE (7*1$) 

15 FORMAT (•DIPROTIC ACID DATA INTERPRETED VIA A MONOPROTIC ACID MODE 
♦L« ) 

CALL YAXIS (0«*0«»1«»15*5) 
WRITE (7*17) 
17 FORMAT ( 'PH* ) 

CALL XAXIS (0«* 0«* 5«* 7* 7) 
WRITE (7*16) 

16 FORMAT ('MILLILITERS OF SODIUM HYDROXIDE*) 

C CALCULATION AND PLOTTING OF THE DIPROTIC TITRATION DATA POINTS 

H ■ 1« 

1 H ■ H - •1*H 

PH > -(ALOG(H) )/2«303 

NUM ■ H*»4 + KHHA»(H»«3) + (KHHA»KHA - KHHA»MHA - KW)»(H»»2) 
♦-(2»KHHA»KHA»MHA + KW»KHHA)»H - KW»KHA»KHHA 

DEN - H»»4 + (KHHA+MM0H)»(H»»3) + (KHHA»KHA + KHHA»MM0H - KW) 
♦♦(H**2) + (KHHA*KHA»MM0H - KW»KHHA)»H - KW»KHHA*KHA 

VMOH ■ -VHA»NUM/DEN 

IF (VMOH) 1*2*2 

2 CALL FPLOT (-2»VM0H»PH) 
5 H ■ H - •2»H 

PH > -(ALOG(h; )/2*303 

NUM « H»»4 + KHHA»(H»»3) + (KHHA»KHA - KHHA»MHA - KW)»(H»»2) 
♦-(2»KHHA»KHA»MHA + KW»KHHA)*H - KW»KHA»KHHA 

DEN ■ H»»4 + (KHHA+MM0H)»(H»»3) + (KHHA»KHA + KHHA»MM0H - KW) 
«#(H«*2) + (KHHA«KHA«MM0H - KW»KHHA)»H - KW«KHHA»KHA 

VMOH ■ -VHA»NUM/DEN 

IF (24«5-VM0H) 10*3*3 

3 CALL FPLOT (1»VM0H»PH) 
CALL FPLOT (-2»VM0H»PH) 
CALL POINTdl 

60 TO 5 

C CALCULATION AND PLOTTING OF A MONOPROTIC TITRATION CURVE ON SAME AXES 

10 MHA « 2t*MHA 
KHA «0»5E-4 



PAGE 



2 



5X H • .l^H 

PH • -(AL06(H) }/2«303 

NUM • H»»3 ♦ ICHA«(H»»2) - (KW ♦ (KHA»MHA))«H - KW»KMA 

OEN • H»»3 ♦ (KHA ♦ MM0H)»(H»»2) ♦ (KHA»MMOH - ICW)»H - KW»KHA 

VMOH • -VHA*NUM/OEN 

IF (VMOH) 91t52*52 

52 CALL FPLOT (1*VM0H*PH) 
CALL FPLOT (-2*VM0H*PH) 

55 M • H - •1»H 

PH - -(AL06(Hn/2«303 

OEM • H»»3 ♦ (KHA ♦ MM0H)»(H«»2) ♦ (KHA»MMOH - t;W)»H - KW*KHA 
NUH • H»»3 ♦ ICHA»(H»»2) - (KW ♦ (KHA«MHA})»H - KW*KHA 
VMOH ■ -VHA*NUM/OEN 
IF(39«-VM0H) 50*53*53 

53 CALL FPLOT (0*VMOH*PH) 
GO TO 95 

50 CALL FINPL 
CALL EXIT 
END 

// XEO 
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DIPROTIC ACID DATA INTERPRETED VIA A WNPROTIC ACID MDDQ 




check on this possibility by incorporating some random scatter 
in the simulated data, and running FORTRAN Program 16 for various 
values of equilibrium constants. Check the IBM instruction 
manual for ways to generate a random variable to use in adding 
or subtracting random "errors" from the simulated data. How 
would you decide what value to use for the monopro tic- acid 
equilibrium constant? 

This is a special case, as presented, but the problem is a 
general one throughout experimental science, and perhaps in all 
areas of knowledge. A careful investigator can often design a 
confrontation that shows a^ model for reality to be invalid, or 
inadequate. Such a model probably cannot be asserted to be 
"true." But the fact that a model is faithful (in all investigated 
respects) to reality does not provide proof or assurance of its 
truth. 
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WHAT IF Kh^- were LARGER THAN K^^ ^? 



The removal of a proton from a molecuxe results in an over- 
all increase by one unit charge in the negative charge on that 
molecule This increase in negative charge (which is equivalent 
to a decrease in positive charge) results in a greater attraction 
(which is equivalent to a reduced repulsion) for remaining protons 
still bound to the molecule. These electrostatic considerations 
would be expected to make the successive microscopic equilibrium 
constants for removal of protons decrease in numerical value. 
The limiting situation would be the case of non- inter xting, 

dentical binding sites in which the sites are so far apart in 
space that electrostatic interactions are neglible; then the two 
successive microscopic dissociation constants are equal » the 
successive macroscopic dissociation constants differ by a factor 
of four, and K^^^ > K^j^. . 

Removal of a proton does more than just increase negative 
charge. Changes in electronic configuration of the molecule are 
expected, and hanges in the three-dimensional shape of the 
molecule are possible. And it may be reasonable to ask the 
question: "What voould happen 1^ Aome change.6 dttZK^d the. motccutz 
60 much that the. Aacce^A-cve mlc/i06coplc e^qulllbKlum cOfUtant6 
IncKza^z In namzKlcal ^^aluo.?" A search through tabulated values 
of macroscopic e^^uilibrium constants has not revealed reports of 
any such instances. Maybe there are no such acids. Or maybe 
there is an alternative explanation. 



CHEMISTS MAY HAVE BEEN We shall present some results of a 

MISLED BY THEIR OWN DATA . . . 

computer simulation that indicate 
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PAGE 1 PARTIAL LISTING OF PROGRAM USED TO PREPARE GRAPH ON PAGE 177 
// JOB T S71007#S51 

SMITH COLLEGE IBM 1130 MONITOR VERSION 2 LEVEL 11 CORE SIZE 8K 
// FOR 

• LIST SOURCE PROGRAM 

VHA - 25.00 
MHA a OaOSOO 
MMOH « 0*1000 % 
KHHA • l«E-$-<-l 

KHA « i,E-3 -^f Diprotic K's 

KW - 1.008E-14 / 

CALL PREP (•2ff ,55* X«» 1.) 

WRITE (7*15) 

IS^FORMAT CDIPRCTIC ACID DATA INTERPRETED VIA A MONOPROTIC ACID MODE 

CALL YAXIS (0«*0«*l«*15f5) 

WRITE (7.17) 
17 FORMAT CPH*) 

CALL XAXIS (0«* 0«» 5«* 7* 7) 

WRITE (7*16) 
16 FORMAT ('MILLILITERS OF SODIUM HYDROXIDE*) 

C CALCULATION AND PLOTTING OF THE DIPROTIC TITRATION DATA POINTS 

H ■ !• 

1 H ■ H - .l^H 

PH ■ -(AL0G(H))/2«303 

NUM - H»»4 ♦ KHHA*(H**3) ♦ (KHHA»KHA - KHHA*MHA - KW)»(H»»2) 
»-(2«KHHA»KHA»MHA + KW»KHHA)«H - KW»KHA»KHHA 
DEN • H»»4 ♦ (KHHA4MM0H)*(H**3) ♦ (KHHA»KHA ♦ KHHA*MMOH - KW> 

* (KHHA»KHA«MMOH - KW»KHHA)»H - XW»KHHA»KHA 
VMOH ■ -VHA«NUM/DEN 
IF (VMOH) 1*2*2 

2 CALL FPLOT (-2»VM0H»PH) 
5 H • H - •2»H 

PH ■ -(AL06(H) )/2«303 

J^^^^ZJi*.*"* * »^HHA*{H»#3) * (KHHA»KHA - KHHA»MHA - KW)»(H«»2) 
♦-(2»ICHHA»KHA»MHA * KW*KHHA)»H - KW»KHA»KHHA 

•2?l!l.r,?**^,J (KHHA+MM0H)»(H»«3) * (KHHA»KHA + KHHA»MM0H - KW) 
Vunl ^^ * (KHHA»KHA»MMOH - KWKHHA)»H - KWKHHA»KHA 
VMOH « -VHA»NUM/DEN 
IF (24.5-VMOH) 10»3»3 

3 CALL FPLOT (l.VMOH.PH; 

CALL FPLOT (-2»VM0H»PH) 
CALL POINT(I) 
GO TO 5 



It^A^'i^Viul.^^^ PLOTTING OF A MONOPROTIC TITRATION CURVE ON SAME AXES 
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the possibility that a chemist could be led astray in inter- 
preting experimental data from titration experiments if the acid 
were indeed diprotic with the second dissociation constant larger 
than the first. The computer-drawn graph on page 177 shows data 
calculated for an extreme case in which the equilibrium constants 
are in the atypical order, and differ by a factor of 100. The 
simulated data are indicated with X marks, and are compared with 
with a continuous curve, drawn on the same axes, calculated for 
a monoprotic acid, with twice the molar concentration, and with 
an equilibrium constant with pK the average of pKfe of the diprotic 
equilibrium constants. (Some other value for the monoprotic K 
might have been more appropriate.) 

If, as has been common in many research laboratories, only 
the half -equivalence point pH value (s) is used in evaluating the 
equilibrium constant(s), then this acid might be incorrectly 
interpretied as being a monoprotic acid with pK = 4, or as a 
diprotic acid with pK values in the normal order, equal to about 
3.7 and 4.2. So maybe some such anomolous acids do exist, but 
have been missed by earlier investigators. 

PREDICTIONS OF CHEMICAL One of the most important features of 
PROPERTIES OF AN ACID 

WITH ATYPICAL ORDERING simulation methods is the ability to 
OF DISSOCIATION CONSTANTS 

make quantitative predictions about 
systems that have never been studied experimentally and which may 
in fact not exist. But the^/te. vot/tt an acid with pK^^^ - 5 and 
pKjj^^ = 3, then we can make some predictions on the basis of the 
information contained in the graphs on pages 177 and 179. The 
following features are worthy of note: 
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f A solution of this acid and its srdium salt would be an 
excellent buffer near pH 4. The slope of the titration curve 
in the interval 5 m£ < V^^^jj < 20 ml is closer' t:> zero than for 
any other of the titration curves in this case study. 

IT Since HA accounts for at most only five percent of the 
A-containing speciss in solution, the titration can be represented 
by the reaction 

H2A i» A° + 2H* 
without serious error in the concentrations of H2A and A". For 
discussion of the equilibrium properties of these solutions, 
considering only the reac;ants and products, there would seem to 
be no need t-j even mention the intermediate species. Can you 
write down the appropriate chemical model? What equilibrium 
Constant would you use? Is there really any difference in the 
predicted shape of the titration curve? 

f A Guess: Probably the shape of the titration curve would 
not be very much different if pK^^- were 2, or 1, or any smaller 
number. For larger values of this equilibrium constant, the 
distribution fraction (HA') /[A] becomes even smaller. A numer- 
ical evaluation of this equilibrium constant from titration data 
becomes increasingly more difficult and less reliable as its 
value become larger. Can you check these guesses? How? 
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Chapter 3 

Polyprotic Acid 
Equilibria 



THE TITRATION OF CITRIC ACID IN AQUEOUS SOLUTION 



Citric acid is the triprotic acid with the structure 




H 0 H 



H 



We shall investigate the chemistry of citrxc acid from the point 
of view of the competition in water for protons. The competition 
for the limited number of available protons is among OH and each 
of the three carboxylate groups on citrate. This detailed exam- 
ination involves a quantitative description of the distribution 
of chemical species throughout the course of a titration of a 
citric acid solution with a solution of sodium hydroxide. Our 
strategy in formulating a plausible, reasonable and verifiable 
description of the chemical system takes the form of the now- 
familiar sequence in which we 

1. Propose a chemical model. This chemical model consists 
of a statement of the chemical species assumed to be 
present, and a collection of balanced chemical equations 
sufficient to permit interconversion of the species. 

2. Formulate a mathematical model which embodies the relevant 
features of the chemical model. This mathematical model 
consists of a set of simultaneous algebraic equations. 

3. Combine the algebraic equations to get an equation which 
relates experimentally-observable quantities, thus 
permitting a comparison of prediction with experiment. 

4. Perform quantitative experiments to obtain numerical data 
for comparison with predictions. 

If the predictions are consistent with experimental data, then 
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we shall consider that the model is an acceptable description 
of the chemical system. With citric acid, we shall find that 
the chemical model that we propose is not unique, and that many 
unreasonable models will give rise to the same prediction of the 
experimentally-observed titration curve. Some chemical judgment 
is required, and some data from non-titration experiments, as 
well as comparative titration data from certain related com- 
pounds, becomes very helpful in making the test of reasonable- 
ness* 

The proton dissociation equilibria for citric acid can be 
schematically diagrammed as in Figure 1. It is possible to 
define even more microscopic equilibria, but the collection of 
explicitly-written reactions and associated equilibrium constants 
given in Figure 2 includes more than enough K*s to describe 
completely the distribution of species in a solution of citric 
acid at any pH value. It includes many more K*s than are re- 
quired to calculate a titration curve, and in fact we shall soon 
find that only three equilibrium constants are required in order 
to dz6cAlbz thz maaKoscopZc tltKatlon phznomQ.na. The many ad- 
ditional equilibrium constants are introduced in order to 
ZnttKpKZt the phenomena. 

PREDICTION OF THE SHAPE For the prediction of the titration 
OF THE TITRATION CURVE 

curve, we need a chemical model written 
in terms of only the macroscopic equilibria. The following 
chemical equations suffice. Written with each chemical equation 
is the associated equilibrium-constant equation that forms a 
part of the mathematical model for the titration. 
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CH2C )(H ) 

HjC i! H2C + H K. = (74) 



(HC'MCH ) 

H2C 4* HC"2 + H Kp » (75) 

^ (H2C') 



(C'')(H ) 

HC"* j» C"' + H Kp » (76) 

^ (HC-2) 



H2O i» H* + OH" = (H'')(0H') (77) 

The two conservation equations are 

IC] = (H3C) +"'(H2C") + (HC^) + (C-') (78) 
[Ml = (M*) (79) 

The electroneutrality equation is 

(H*) + (M*) = (OH') + (H2C') + 2(HC-^) + 3(C-') (80) 

An equation for the titration curve can be obtained by using 
Equations 74 — 79 to eliminate the individual species concentra- 
tions from Equation 80, producing an equation in which only the 
concentrations [C] , [M] , and (H*) appear; and then introducing 
volume variables. A great deal of algebra is simplified if we 
ja&kQ use of the results of the next section on distribution 
fractions. So we interrupt this derivation for a digression. 

DiSTRIBUTION FRACTIONS IN TERMS Proceeding exactly as we have 
OF THE MACROSCOPIC EQUILIBRIA 

before in the cases of monc - 
protic and diprotic acids, we find the distribution fractions by 
starting with Equation 78, and then using the equilibrium-con- 
stant equations to eliminate all but one of the species concen- 
trations. Thus we can write Equation 78 as 
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Figure 1 

A SCHEMATIC PROTON DISSOCIATION SCHEME FOR CITRIC ACID 
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Figure 2 

MICROSCOPIC EQUILIBRIA IN THE CHEMICAL MODEL FOR CITRIC ACID 
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IC] = (H3C) + (H3C)K^/(H*) + (H3C)K^Kg/(H*)2 + (H3C) K^Kj^K^,/ (H*) 
- (H3C){1 + K^/(H*) + K^Kp/(H*)2 + K^KpK^7(H*)n 
(H*)'(C1 = CH3n{(H*)' + K^(H*)2 + '^'a'^IjC'*^ * ^aVc^ 181) 
In like manner, we get Equation 78 in terms of the concentration 
(H2C') to give 

Ka(H*)MC] = (H2C"){(H*)' + J^A^"*^' * Vb^"*) * hVc^ ^^^^ 
It is significant that the quantity in { } is the same in both 
Equations 81 and 82. In fact, this same factor appears— again 
when we write Equation 78 in terms of (HC"^) or in terms of 
(C'). Let us therefore define this factor as F^jg^^* 

^den = ' h^^'y ' Vb("') * Vb^ ^^^^ 

and then write each of the distribution fractions as a fraction 

with the denominator • The four resulting equations are 

den 

(H,C) (H*)» (HjO (hVKa 

(HC-') (H*)K^Kg (C) _ K^KgKg 

ICl " ' "der, 

RETURN TO THE We now use Equations 77, 78, 79, and 84 to trans- 
TITRATION CURVE 

form Equation 80 into 

(H^) * [M] • K^/CH*) -{(H*)'K^ * 2(H*)K^K3 ♦ SK^K^Kj^l/F^^^ 

Multiplication of both sides of the equation by ^^^^^ followed by 
rearrangement of terms , yields 

(H*)» * (H*)M[M) * K^} ♦ (H*)MKa[M) - K^[C) ♦ K^Kg - K^) 

* (H*)MVb[M) - 2K^Kg[C] * K^KpKc - K^K^} (gS) 

* (H*){K^K3Kj,(M) -3K^KgK^,(C) - V^K^l ' VgKj,K^ - 0- 
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Recognizing that the two concentrations [M] and [C] are not 
independent during a titration, we decide to change variables 
so as to have a single independent variable. The natural 
choice for this independent variable is the volume of titrant 
which has been dispensed from the burette. The change of vari- 
ables is accomplished by utilizing the relationships 

[M] . ^MOH^OH . VhA (86) 

^MOH * ^HA ^MOK * ^HA 

Introduction of Equations 86 into Equation 85 results in 



*mom""*ha 



* (h*)MVb«hoh ♦ hhh ■ W 

Equation 87 is in the usual form which we have been using for 

computer-assisted calculation of titration curves. 

Without any consideration of the microscopic equilibria, we 

have been able to make a quantitative prediction of the data 

which would be obtained in an experimental titration. The 

description of the observable phenomena has been made in terms 

of Equation 87 which is a phcnomznologlcaJt equation written in 

terms of four phe^nomznotoglcat equilibrium con6tant6. We consider 

as synonyms the terms 

phznormno logical tqulllbrlum con6tant 
macro6coplc dqulllbKlum con6tant 
txptKlmtYitcilly-ob6zn.v(iblt tqulllbfiium con6tant 
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Each macroscopic equilibrium constant may describe a 
process that is a superposition of several of the microscopic 
chemical reactions of the original model. The microscopic reac- 
tions ^ and the individual chemical species used in formulating 
the model for the reaction of protons with the variously-proton- 
ated citrate ions, are not directly, explicitly or unambiguously 
observed by casual inspection of the titration curve of citric 
acid. To analyze a titration curve, we first decide to treat 
citric acid as a triprotic acid, and obtain Equation 87. We 
then extract numerical values of three macroscopic equilibrium 
constants from an experimental titration curve. Then, supplying 
additional information not contained in the titration data, we 
try to interpret these three macroscopic equilibrium constants 
in terms of the many reactions in Figure 2. 

Equations 74, 75, and 76 are not 
written with the same symbols as the 
dozen equations in Figure 2. We 
need to make a translation between the macroscopic language and 
the microscopic language, and we now do so via the definitions: 

(H3C) 



RELATIONSHIPS BETWEEN THE 
MACROSCOPIC AND THE 
MICROSCOPIC EQUILIBRIUM 
CONSTANTS FOR CITRIC ACID 



sum of the concentrations 
of all triply-protonated 
citrate species 



(H2C') = sum of the concentrations 
of all doubly-protonated 
citrate species 



(HC^) = sum of the concentrations 
of all singly-protonated 
citrate species 













H 






+ 


H 










IhJ 
















H 


+ 















(88) 



(89) 



(90) 



(C) H 



sum o£ the concentrations 
of all unprotonated 
citrate species 



(91) 
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Direct substitution of definitions 88 — 91 into Equations 74, 
75 and 76 results in 



fHl 










il 
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H 










[hJ 



H 
H 



(HC"2)(H*)2 
(H3C) 



(92) 



r 




















+ 


H 


+ 





















(93) 



(C ')(H^) 
(HjC) 



(H ) 



(94) 



Inspection of Figure 2 reveals various ways in which the micro- 
scopic equilibrium constant equations can be combined with Equa- 
tions 92, 93, and 94. For instance, we find that 



= K, + K2 + K3 

K^Kg = K,K^ + K2K6 + K3K, 
= KiKw + K1K5 + K2K8 
= K2K6 + K2K8 + K3K7 

» other combinations of microscopic equilibrium constants 
K^KgK^ — KiKi^Kio 
= K2KgKio 

« many other combinations of microscopic constants 
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There are more independent microscopic equilibrium constants 
than there are macroscopic equilibrium constants, and therefore 
experimental titration data alone contain insufficient informa- 
tion for numerical evaluation of the microscopic constants. There 
are many sets of values of the microscopic constants that are 
fully consistent with the experimental values for the three macro- 
scopic equilibrium constants. 

DISTRIBUTION FRACTIONS The fraction of all citrate-containing 
IN TERMS OF THE 

MICROSCOPIC EQUILIBRIA species present as the microscopic 

species 

H 
H 

is given by the first .of Equations 84, since by definition 88 

= (HaC) 



We thus have 



HsC 



H 
H 

H 

[C] 



den 



(95) 



Combination of Equation 95 with the defining equations for Ki , 
K2, and Ks from Figure 2 gives . * 



(96) 



fH 



H 



%C2 '= — = FHaC'^^/f"*) 



(97) 
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• Vlnd tquatlon6 ioK computing tht dlhtKlbutlon iKattlon^ 
ioK tht Ktmainlng ml<iK06<Loplc 6ptclt6. CompaKt youK Kt6utt6 
with tht tquatlon6 u6td In FORTRAM PKogKam 17 ^OK catcalatlng 
tht6t dl6tKlbutlon ^Kactlon6 . 

A FORTRAN PROGRAM FOR CALCULATING FORTRAN Program 17 gives a 
THE MICROSCOPIC DISTRIBUTION 

FRACTIONS DURING THE TITRATION prediction of each ot the 

OF CITRIC ACID WITH SODIUM HYDROXIDE 

microscopic distribution 
fractions during a simulated titration of citric acid with a 
solution o£ sodium hydroxide. Numerical values of seVen of the 
microscopic equilibrium constants are introduced. The molecular 
symmetry of citric acid is invoked, allowing us to equate each 
of the remaining five micro constants with one of the initial 
seven. Equations such as those developed on pagp 190 are used to 
calculate values for the three macroscopic equilibrium constants, 
and these constants are used with Equation 87 for the calculation 
of the titration curve. At each pH value on the titration curve. 
Equations such as 96, 97, or 98 are used to calculate each of the 
microscopic equilibrium constants. 

Format statement 11 has two programming features that make 
the output a little neater and more legible. In our previous 
format statements for the printer beginning with information 
within * * marks, the first character with the ' ' has been a 
blank. This time the numeral 1 appears, and this 1 is an instruc- 
tion for the printer to advance to a new page of paper. Appearing 
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// FOK 

* Llbf SOUKCl HkOvsKAi'' 
*I0CS(CAKD» 1132 PKlNTtK) 

HEAL .MUM» K.A» KB* KC» K.W» MMOH» .VHA » Kit K2 » K3» K.4 » f.5 1 K,6» K.7» 
«Kfl» K9» K10» Kilt K12 
C CALCULATION OF FRACTIONS OF TOTAL C I TKAT h-COM FA I NIi\G SPECIhS 

C PREStNT AS hACH OF THt£ MICROSCOPIC SP£CIbS DUKINu THt TITRATION 

C OF CITRIC ACID 

Kl = 10**(-3,e5) 

K2 = 10**(-3.35) 

K4 3 10**(-4,60) 

K5 ■ 10**(-4.40) 

K6 « 10**(-5,10) 

KIO * 10**(-5.85) 

Kll » 10**(-6.05) 

K3 = Kl 

!«;7 = K5 

K8 = K6 

K9 = K4 

K12 = KIO 

KW = 1.008E-14 

KA « 2.*K1 + K2 

Kb = (K1*K5 + 2.*K2*K6)/KA 

KC « (Kl*iC4*K10)/(KA*Ka) 

PKA « -(ALOGCKA) )/2.303 

PK8 = -{ALOG(KB) )/2.303 

PKC « -(ALOG(KC) )/2.303 

WRITE (3»11) PKA» PKB» PKC 
11 FORMAT CIPKA F6,2» • PK8 =•, F6.2. ' PKC »•» F6.2//) 

VHA = 10.00 

iMHA = ,100 

MMOH = .100 

H = 1, 

1 H = H - •1*H 

NUr': = ti**b + (H**4)*f;A + (M**3)*{NA»Kb - K;v - KA*MHA ) + 
*(H**2)*(KA*Kd*KC - MV*.<A - 2.*K.A*Ke*r/nA) - H*(KA*Kb*KW + 
*3.*KA*i<.B*KC».vlhA ) - KA *:Cb>KC*KW 

OCn » (H**5) + «H»*4)«(MK0H + KA ) + ( H** 3 ) * ( -^A«.>- :-.UH - KW + K.A«k.B) 
»+ (H**2) *(KAa!<.B*,v,.v,GH - KK*KA + KA«Ku*nC) + H* ( K.A*i<.8*KC*f-iM0H - 
*KA*KB«KW) - r'xA«K.b*KC*Kw 

V.VOH = -VHA*i\Uf-i/Dtt\ 

IF (VMOH) 1»2»2 

2 fOLH = (H*«3) + r;A*H*H + KA*KB*H + KA*K,lj«l<;C 
FH3C = (H**3)/F0EN 

f-H2C = KA*(H«*2 )/FDEi-4 
FH2C = KA*{H*#? )/FDEN 
FHC = .<A*KB*H/FDEN 
FC = KA*KB*KC/FOEN 
FH2C1 = K1*FH3C/H 
FH2C2 = K2*FH3C/H 
FH2C3 = K3*FH3C/H 

FHCl = FC*H/K10 193 
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FhC? = PC»H/K11 
FHC3 = FC»H/K.12 
PH c -(ALOGCfi) ) /2.30^ 

wRITfc(3tl2) fH3C» HH2C1* hH2C2» HI2*^3» i=Hv.i» h-iCZ* PhC3» »-C» V^Om» 

1? FOKf^AT (F9,5t3x.3H9,t)»3X,3F9.i)»3X,h9.5»3X.Ffi.2.f-3,2) 
H = H -•!>«H 

.<Ui'. = H««5 + (!-:««4)*KA + (H**3)«{KA*Kb - K a - K.A«MHA ) + 
♦ {H««2)«(KA«K.b»^C - ^■a'*K.A - 2««K.A«i^B*v,HA) - H» t i;A«K.B*S.W + 
»3«*K,A*Klj*K.C*i«'HA ) - K.A«K,b*KC*K-W 

DEN = (H**5) + (H««4)«(MX0H ~ K.A ) + ( H««3 ) * ( KA«vj..OH - KW KA«»iB > 
»♦ (H«*2 )*lKA*K.ti*MMOH - KW*kA + |iiA*l<^b*rNC) + H* ( <v.A»:xd*KC*f-;.' Oh - 
»KA«K.B«K.W) - K.A«K.b-'*KC«!<;W 

VMOH = -VHA«NUy./DEN 

IF (VMOM) 30»2»2 
30 CALL EXIT 

END 



// XEO 
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at the end of the statement is //, telling the printer to advance 
two lines . 

The computing capacity of the IBM 1130 system is exceeded 
during the calculation of NUM and DEN for very small values of 
(H*). As an example, at pH 11, (H*^) = 10^^^ and (H"*") ^ = lO"". 
This computer cannot handle such a small number, but happily the 
failure is fail-safe failure for us, since the computer treats 
the very small number as zero. This would pose serious problems 
for some calculations, but not for us at this time. Each of these 
very small numbers enters as an insignificantly-small term in an 
indicated sum, the sum containing other terms that are properly 
calculated, and no significant error results in the sum. 

The output would be easier to read if a caption had been 
printed over the columns of numbers. Such a caption could have 
been printed by using a WRITE statement just after format statement 
11. 

ABOUT COMMONALITY OF The mathematical model for thir^ system is 
CHEMICAL SPECIES 

formulated by describing the equilibrium 
system as involving many independent chemical reactions, each 
reaction characterized by a microscopic equilibrium constant. 
Coupling of these independent chemical reactions is accomplished 
through the requirements of 

HConservation of mass 

IfElectroneutrality 

IFCommonality of chemical species 
The phrase commonality 6pe,cie.6 means that a chemical symbol 
has the same meaning in each of its appearances in chemical equa- 
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tions, and that the corresponding concentration has the same 
numerical value in each of the algebraic equations in which it 
appears. This means, among other things, that the pool of that 
species is available at the same concentration to all reactants. 
When that species is produced by a reaction, it becomes distributed 
throughout the solution. The requirement of coramonality of 
species might fail to be satisfied in an irtra-molecular reaction 
during which a group was transferred between two regions of a 
molecule without ever being released into the main bulk of the 
solution as a free and independent species. 

A valuable research paper that should be read in conjunction 
with this section is R. B. Martin, "A Complete Ionization Scheme 
for Citric Acid," The Journal of Physical Chemistry ^ vol. 65 » 
pp. 2053 - 2055 (l96l). Details are given of the experiments and 
the assumptions required for estimating values of the microscopic 
equilibrium constants for citric acid. 

Equilibria involved in the titration of the amino acid 
tyrosine involves twelve microscopic equilibrium constants, and 
all twelve have been evaluated from experimental data and a set 
of plausible assumptions. Details are given in R. B. Martin, 
J. T. Edsall, D. B. Wetlaufer, and B. R. Hollingworth , J ournal 
of Biological Chemistry , vol. 233, pp. ll*29ff (1958) 
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DESIGNING YOUR OWN RESEARCH PROJECT INVOLVING 
CARBOXYLIC ACIDS AND MULTIPLE EQUILIBRIA 

One way to plan a research project is to begin by reading 
some relevant research papers relating to the system that you 
intend to study. Some excellent sources of references to such 
papers are given on page 61. Many common diprotic and polyprotic 
acids have been studied only in passing, as parts of research 
projects that other orientations and goals. Only a small number 
of polyprotic acids have been studied in sufficient detail so 
that all microscopic equilibrium constants can reliably be 
assigned numbers. Even fewer have been studied with a view to 
obtaining enough independent data so that assumptions made can 
be checked for consistency. 

Another way to plan a research project is to pick a system 
that has never been studied carefully, and perform enough pre- 
liminary experiments to get some crude data. Then, with these 
data before you, you can design appropriate experiments that can 
be interpreted by some chemical model. 

It is important that you subject your data and your inter- 
pretative model to stringent critical examination. Many of the 
experimental studies reported in chemical research journals are 
flawed in this respect. Uses to which computers have been put 
on preceding pages of this case study are typical of the sorts 
of quantitative confrontations of model and experiment that can 
provide convincing objective tests. And you have the tools 
needed to make such tests 4 

It is often helpful to think about the data you would expect, 
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if the tentative model you propose were indeed operative. Com- 
puter simulations can be helpful in this regard, allowing you to 
determine if certain experiments could be helpful, even in prin- 
ciple. 

Many sorts of chemical systems can be informative, challenging, 
or just a little bit different, and can be fun to investigate. 
You may want to find out what is known about the carbon dioxide — 
bicarbonate — carbonate equilibria, both in distilled water, and 
in systems such as blood, or soda pop. You may wish to unravel 
the multiple equilibria in a polyprotic acid such as EDTA, EGTA, 
or another of the common complexing agents. Mixtures of acids 
are important in many ways: can you design a mixture of acids 
that v/ould make an especially-good buffer over a wide range of 
pH, or devise a way of titrating two different carboxylic acids 
in the same solution, obtaining values for the concentration of 
each? What is the effect on the observed titration behavior of 
various caids when various metal ions are added to the solutions? 
How can these effects be explained? Is there ever any difference 
in the observed equilibrium constants between various optical 
isomers (the left-hand and right-hand varieties of isomers) , or 
mixtures of both? Or of other kinds of isomers? 

Use your imagination, be creative, and show yourself (and 
maybe your instructor, and perhaps the wider world of other 
scientists) that you can make some critical, informed, and sub- 
stantiated judgments about the nature of carboxylic acids in 
solution . 
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APPENDIX I 

PROGRAMING FOR THE PLOTTER 

STANDARD FORTRAN STATEMENTS 

DRAWING A STRAIGHT LINE: FPLOTCl >X,Y) 

The smooth solid curves drawn by the plotter are in fact 

constructed by instructing the plotter to draw a series of very 

short straight-line segments. The procedure is to bring the pen 

to one end of the desired line, with the pen up so that the pen 

is not writing. The pen is then lowered, and moved while in the 

down position to the other end of the straight line; a straight 

line is thereby drawn. An instruction to move to still another 

position yields a continuation of the line if the pen is down.. 

The necessary information is given to the computer in the form 

of FPLOT statements of the form 

CALL FPLOT (I^X^Y) 

I is an integer controlling the up-down pen 
position as follows': 

= 0 => no change 

= positive => control pen before movement 

= negative => control pen after movement 

= odd =^ raise pen 

= even =^ lower pen 

X and Y are the coordinates of the graph position 
to which the pen is moved. 

PRINTING LABELS ON THE GRAPH 

If you want to place labels on the graph, use the FPLOT 

statement to bring the pen to the position where you wish the 



This appendix is reprinted from TEACHING rUIDE TO ACCOMPANY 
CHEMICAL REACTION MECHANISMS by George M. Fleck. Copyright © 1971 
by Holt, Rinehart and Winston, Inc* weprinted by permission of 
Holt, Rinehart and Winston, Inc. 
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lower left corner of the first character to be* The pen should 
be in the up position. Then use the statement 
WRITE(7,J) 

where J is the number of a format statement of the form 

d FORMAT C *) 
The desired label goes between the ' ' marks • 
PLOTTING A POINT ON THE GRAPH 

If you wish to pl^t a point on the graph, you may use the 
statement 

CALL POINTCl) 

This statenTent cause the plotter to draw an X mark at the 
current pen position. It expects to find the pen down and it 
leaves the pen down when finished. 

SPECIAL PLOTTING SUBROUTINE 

[These subroutines were written by Professor Bruce 
Hawkins, Department of Physics, Smith College.] 

PREPARATION SUBROUTINE! PREPCXS . YS >X0 . YO) 

This subroutine determines the scale of your graph, sets 

the position of the origin, and leaves the pen in a convenient 

position to write a title. 

XS is the scale you have chosen for the x-axis, 
in inches per unit; 

YS is the scale you have chosen for the y-axis, 
in inches per unit; 

XO is the distance in inches from the left-hand 
edge of the graph to the origin; 

YO is the distance in inches from the bottom of 
the graph to the origin. 

PREP types a reminder to set the pen position and waits for you 

to do so. PREP should be called before any other plotting 
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subroutine. After calling it, a title may be written using 
an ordinary FORTRAN WRITE statement and its associated FORMAT 
statement. Since this subroutine uses the typewriter, your IOCS 
card must mention the typewriter. 

A listing of this and each of the other special plotting 
subroutines appears as the conclusion of this appendix. 

SUBROUTINE YAXI SCX . Y >U , N ,NLAB ) 

This subroutine draws the y-axis the length you have chosen, 
draws evenly-spaced tick marks, places a numerical scale beside 
the line, and leaves the pen in an appropriate position to 
write a label for the axis. 



X and Y are the coordinates in your own units of the 
point where the y-axis begins; 

U is the distance in your own units between tick marks 
on the axis ; 

N is an integer such that NxU is the length of the axis; 

NLAB'is the number of tick marks which are to be 
l^^.belled with the value of the scale at the 
point. NLAB should be equal either to N or 
to N divided by an integer. 



This subroutine should normally be called before the XAXIS 
subroutine. It leaves the character writing process set to 
write characters on their side. If you follow this statement 
with a FORTRAN WRITE statement, you can label the axis. 

SUBROUTINE XAXI SCX , Y ,U . N>NLAB) 

This subroutine behaves exactly like the YAXIS subroutine and 
like that one can be followed by a WRITE statement to place a 
label on the axis. 
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FINISHING THE GRAPH! FINPL 

At the end of a graph, or after plotting several graphs 
on the same axes, it is convenient to use the statement 
CALL FINPL 

This subroutine plots a cross at the lower left hand cotner of 
the graph pn top of a plus sign put there by PREP* This is a 
check on machine errors during the plotting* If the two symbols 
do not land exactly on top of each other, the graph should not 
be trusted* FINPL then raises the pen and moves it to an 
appropriate position for beginning a new graph* 
NORMS^ SAVES. GETS 

These three subroutines are listed because they were used 
within the plotting subroutines* 
ON GETTING THE DESIRED GRAPH SIZE 

These programs were used with an IBM 1627 x-y"strip chart 
plotter* The system used produced a half-size graph, and so we 
introduced some factors of 2* into the PREP program* For the 
standard system, one replaces the statement 

CALL SCALF (2 * *XS^ 2 * *YS^ -XO/XS^ - YO/YS ) 
with the statement 

CALL SCALF (XS^ YS^ -XO/XS^ -YO/YS) 
in the PREP program listed on page 86* The student user is never 
aware of these considerations* 
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// FCR 

«nMr riTEGC-^s 

♦LIST SOURCE PROGRA.v 

SUBROUTINE PREP ( XS » YS » XO » YO ) 
C SETS SCALE. ORIGI.%. AND LOCATES P£.\ TO ft A TIILE 

C IT PLOTS A REGISTlR POINT AT THE I.NITIAL POSITION OF THE PLi\ AuO 

C IT ADJUSTS THE SCALE TO THE HALF SUE STEPS OF TmIS HLOTTEU 

C XS» YS ARE X AND Y SCALES 1 1\ liNiCHh.S PER USER'S UiNlT AiMD .'.UST lit 

C VARIABLES SINCE PREP CHANGES THEM 

C P'iEP ASSUMES THAT THE PEi\ IS AT THE LOwER LEFT HAND COKwER OF 

C THE GRAPH 

C XO IS THE DISTANCE IN INCHES FROM THE LEFT HANw SIDE OF THE 

C GRAPH TO THE ORIGIN 

C YO IS THE DISTANCE IN INCHES FRO|V THE BOTTOM OF THE GKAPh TO THE 

C ORIGIN 

C 

IF (r\SW-15324) 6»10»6 
6 vV'^ITE (1»1) 

1 FORMAT ( 'PLACE PEN AT RIGHT EDGE OF PLOTTER AND PUSH START') 
PAUSE 

CALL SAVES ( XS » YS » XO » YO ) 
10 CALL SCALF ( 1 . » 1 . » 0 . » 0 . ) 

CALL FPLOT ( -2 15 » .05 ) 
00 5 1=1.20 

CALL FPLOT(0»-.15».15) 
CALL FPLOT (0.-.C5i.l5) 
CALL FPLOT (0»-.05».05) 
5 CALL FPLOT (0»-.15».05) 

CALL FPLOT ( 1 .0. »0. ) 
CALL FPLOT (-2»0.».2) 
CALL POINT (0) 

CALL SCALF ( 2 . *XS » 2 . *YS »-XO/XS »-Y0/ YS ) 

CALL FCHAR ( ( 1 .-XO ) /XS » ( 10 .-YO ) / YS » . 2 » . 4 » 0 . ) 

NSV; = 16324 

RETURN 

END 

// DUP 

♦ STORE .WS„ UA PREP 



Written by Prof. Bruce Hawkins, Department of Physics, 
Smith College. Used by permission. 
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// FO"? 

*ONE W0'.0 INTEGERS 
♦LIST SOURCE PROGRAy 

suer^ouTiME yaxis (x.y.u..\.nlab) 

C PfTS A STANDARD YAXIS ON PLOT STA.xTiNG AT X.Y. rtlTh in SEGi>'itiMrS 

C OF LENGTH U SEPARATED BY TICK MARKS A.\D LABELS \LAB +1 OF THE 

C TICKS WITH THE SCALE VALUE AT THAT POINT. X AND Y ARE l\ USERS 

C UNITS AS DEFINED BY X5 AND YS a^HICH ARE INCHES/ USER'S UaIT 

C ^'UST BE PRECEEDED 3Y A CALL TO ?'-?EP OR SCALF 

C LEAVES PEN POSITIONED FOR LEGEND 

C 

C'aLL gets (XS.YS.XO.YC) 
CALL FGRID (l.X.Y.U.N) 
CALL FChAR (X. Y. .2. .3. 1.570795) 
IE = NLAri + 1 

TENLG = .A342945*AL0G( Y+u*FLOAT(.N) ) 
DO 5 1=1.1- 

YP = Y+U*FLOAT( ( I-1)*(N/MLAB) ) 
CALL FPLOT ( 0 » X- . 1/XS . YP- . 5/ YS 5 
IF (AfiS( YP*YS)-.05) ^.4.6 

4 aRITE (7.10) 

10 FORv.AT (5X.'0') 

GO TO 5 
6 CALL NORM'S ( TENLG. YP) 

5 CONTINUE 

CALL FPLOT ( 0 .X-.3 /XS . Y+ ( U»FLOAT ( N ) ) / 4 * ) 

RETURN 

END 

// DUP 

♦STORE \'!S UA YAX.IS. 



Written by Prof. Bruce Hawkins, Department of Physics, 
Smith College. Used by permission. 
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// FOP. 

*Q'AE 1%0'?D I.>{TEGEK3 
»L1ST SOURCE PROGRA.-i 

SUaKOUTINE XAXIS ( X . Y »U..\ ..\LAd ) 

C PUTS A STA.\DARD XAXlS ON PLOT STARllNo' AT X.Y» aITh :,Lo.^.L.\Ti) 

C OF LENGTH U SLPAHATED BY TICK. MARKS AM; LAIJLLS .NLAli +1 UH InL 

C TICKS WITH THE SCALE VALUE AT THAT P0I.->1T. X AND Y ARE I A Uiitiii) 

C UMIT5 AS DEFINED BY XS ANO YS WHICH ARE INCHES/ USER'S UNIT 

C .'-UST BE PRECEEDED BY A CALL TO PKEP OR SCALF 

C LEAVES PEN POSITIONED FOR LEGE><D 

/" 
V. 

CALL GETS(XS»YS.XO.YO) 
CALL FGRID(0»X.YtU.N) 

Call fchar (x. y. .2.. 3.0.) 

IE = .\LAB+1 

TEiMLG = .4342945*AL0G(X+U*FLCAT(i\) ) 
DO 5 1= l.IE 

XP = X+U*FL0AT( ( I-l )*(N/NLAB) ) 
CALL FPLOT (O.XP-,5/XS.Y-.25/YS) 

IF (ABS(XP*XS)-.05) 4.4.6 

4 WRITE (7tl0) 

10 FOR«AT (5X. •0' ) 

GO TO 5 

6 CALL NORMS (TENLG.XP) 

5 CONTINUE 

CALL FPLOT ( 0 . X+ ( U*FLOAT ( N J ) /4. . Y-. 5 / YS ) 

RETURN 
END 

// DUP 

♦STORE WS. UA XAXIS 



Written by Prof. Bruce Hawkins, Department of Physics, 
Smith College. Used by permission. 
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// FOR 

»Oi\e W0:^0 INTEGERS 
»LIST SOURCE PROGRAM 

SUBf^OUTINE FIiNPL 

C FINPL PLOTS A REGISTER POINT IN TOP OF THE O.NE .MADE bY PREH ANi; 

C -MOVES THE PEN TO BLANK PAPER (ASSUN.ING AN d.5 3Y 11 OKAPnJ* 

C 

CALL GETS ( XS»YS»XO»YO) 

x=-xo/xs 

Y=-YO/YS 

CALL FPLOT (1»X»Y) 
CALL FPLOT (-2»X»Y) 
CALL POINT (1 ) 

CALL FPLOT (1 »9, /XS+X » Y-. 2/ YS ) 

RETURN 

END 

// DUP 

♦5T0SE WS UA FINPL 



Written by Prof. Bruce Hawkins, Department of Physics, 
Smith College. Used by permission. 
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// FOR 

*ONE WORD INTEGERS 

SUBROUTINE NORKS (TtNLGiP) 
C WRITES SCALES ON AXES ADJUSTING THE FORMAT TO THE SIZE OF Tut 

C SCALE 
C 

NPOW « IFIX (TENLG+.G!) 

P « P+SIGN { •00005*{ 10«**NPOW) f P) 

IF (TENLG-3.^) 5»5»30 

5 IF (TENLG+2. ) 30*6»6 

6 NPOW « NPOW + 3 
NPsP 

GO TO (11»12»13»14»15»16) »NPOW 

11 WRITE {7»21) P 
21 FORMAT {F7.4) 

GO TO 50 

12 WRITE (7»21) P 
GO TO 50 

13 WRITE {7»23) P 

23 FORMAT (F6.3) 
GO TO 50 

14 WRITE (7*24) P 

24 FORMAT {F6. 2) 
GO TO 50 

15 WRITE {7»25) P 

25 FORMAT {F6,l) 
GO TO 50 

16 WRITE (7»26) NP 

26 FORMAT (15) 
GO TO 50 

30 WRITE {7»35) P 
35 FORMAT (E10.3) 
50 RETURN 
END 

// DUP 

*STORE WS UA NORMS 



Written by Prof. Bruce Hawkins, Department of Physics, 
Smith College. Used by permission. 
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// ASf-' SAVES AND GETS ROUTI.^lE TO STOKE SCALE Ai\.D lt><0 FOK PLcT .iwuTI^ts* 
»LIST 





LiNT 




SAVtS 






L\T 




GETS 






DFC 




1. 


DEFAULT VALUES 


YS 


DEC 




1. 


DEf-AUL! VALUtb 


XO 


DEC 




0. 




YO 


DEC 




0. 




SAVES 


DC 




* 






SIX 


1 


SVXl+1 






STX 


2 


SVXl+3 






LDX 


11 


SAVES 


GET PA.^AM ADDR 




MDX 


1 




FOR Ii'iDEXlD PIC<-UP 




STX 


1 


SVLP+i 






STX 


1 


SVOUT+1 


IS ALSO EXIT ADJRESS 




LDX 


1 




FOR PICK UP LOOP 




LDX 


2 


0 




SVLP 


LDD 


11 




FILLED BY END OF PAKAM AUUK 




STD 


L2 


xs 






yDX 


2 


2 






MDX 


1 


1 






MDX 




SVLP 




SVXl 


LDX 


LI 








LDX 


L2 






SVOUT 


BSC 


L 


*-* 


EXIT ADDRESS SUPPLIED 


GETS 


DC 










STX 


1 


GTXl+1 






STX 


2 


GTXl+3 






LDX 


11 


GETS 






MDX 


1 




LIKE SAVE PROCEEDUKt 




STX 


1 


GTPR+1 






STX 


1 


GTCUT+1 






LDX 


1 








LDX 


2 


0 




GTLP 


LDD 


L2 


XS 




GTPR 


STD 


U 




DESTIfMATlOiM ADDR FILLED 




.VDX 


2 


2 






,MDX 


1 


1 






VDX 




GTLP 




GTXl 


LDX 


LI 


*"* 






LDX 


L2 


*»* 




OTOUT 


bSC 

END 


L 


»-# 


EXIT ADDRESS SUPPLltO 



Written by Prof. Bruce Hawkins, Department of Physics, 
Smith college. Used by permission. 




APPENDIX II 

NATIONAL BUREAU OF STANDARDS BUFFERS 
FOR CALIBRATING A pH METER 

Before using a pH meter for measurement of pH, the meter 
must be cal3,brated using a solution of accurately-known pH. 
The information in this appendix, giving instructions for pre- 
paring such standard solutions and the pH values of those solu- 
tions, is taken from the following publications of the United 
States National Bureau of Standards: 

flPotassium Tetroxalate (pH = 1,679 at ZS'^C) , N,B,S, Certif- 
icate for Standard Sample 189, January 10, 1964- 

UPotassium Hydrogen Tartrate (pH = 3,557 Csaturated solution], 
pH » 3,639 CO. 01 molal solution: at ZS^'C) , N,B,S. Cextificate for 
Standard Sample 188, January 10, 1964 • 

fPotassium Hydrogen Tartrate (pH = 4,008 at ZS^'C), N,B,S. 
Certificate for Standard Reference Material 185d, July 21, 1966, 

f Potassium Dihydrogen Phosphate and Disodium Hydrogen Phos- 
phate, equi-molal mixture (pH = 6.862 at 25*'C), N,B,S, Certifi- 
cate for Standard Reference Materials 186-I-c and 186-II-b, May 
1, 1969. 

^Borax (pH = 9,180 at 25'*C) , N,B.S, Certificate for Standard 
Sample 187a, January 10, 1964, 

For highest accuracy, N,B,S, Standard Reference Materials should 
be used and the directions followed exactly. Comparable reagents 
and comparable procedures will give satisfactory results for pH 
equipment which can be read to ±0,01 pH units, since the greater 
precision of the N.B,S, values is not needed, 

POTASSIUM A 0,05-inolal solution of potassium tetroxalate di- 
TETROXALATE 

hydrate is recommended by use as a pH standard. To 
prepare the solution, weigh 1,261 grams of the salt into a 100-mil 
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volumetric flask. Dissolve the salt and fill to the mark with 
distilled water. The salt should not be dried before weighing. 
Values of the pH of this solution at various temperatures are 
given in the following table: 

°C pH °C pH pll 



0 1.666 30 1.683 55 1.715 

5 1.668 35 1.688 60 1.723 

10 1.670 38 1.691 70 1.743 

15 1.672 40 1.694 80 1.766 

20 1.675 45 1.700 90 1.792 

25 1.679 50 1.707 95 1.806 



The uncertainty of these values is estimated not to exceed 
±0.005 pH units for temperatures from 0 to 60°C and ±0.01 units 
from 60 to 95°C* The liquid- junction potential of the common 
H cell displays a considerably greater cariability in solutions 
of less than 2.5 than in solutions of pH between 2.5 and 11.5. 
For this reason, experimental pH values may differ by as much 
as 0.02 to 0.05 pH units from the values given above. 

POTASSIUM A solution saturated with potassium tartrate near 25^C 
HYDROGEN 

TARTRATE is recommended as a standard for the calibration of 

pH equipment at temperatures between 25 and 95°C. A 0.01-molal 

solution is also recommended as a pH standard for the temperature 

range from 0 to 60^C. 

The saturated solution is prepared by adding an excess of 

the salt to distilled water contained in a glass-stoppered bottle 

or flask. Then shake vigorously. With a 100 percent excess of 

the salt, a few minutes of shaking is sufficient for saturation. 
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One hundred ml of water will dissolve about 0.7 g of the salt 
at 25**C. Allow the solid to settle and decant the clear solution, 
or filter if necessary. Store the solution in a glass-stoppered 
Pyrex or Kimax bottle. For an accuracy of ±0.001 pH units, the 
temperature of saturation must lie between 24 and 26°C. 

To prepare the C.Ol-molal solution, weigh 0.1878 grams of the 
salt to a lOO-mA volumetric flask. Fill to the mark with dis- 
tilled water having a conductivity not exceeding 2 x 10'^ ohm'^cm'^ 
at 25**C. The salt need not be dried before use. Shake well. 

These tartrate solutions are very susceptible to mold growth 
which is usually accompanied by an increase of a few hundredths 
in the pH value. For the most accurate results, tartrate stand- 
ards should be prepared fresh each day. If they must be kept 
longer, they should be stored in a refrigerator. 

Values of the pH of the saturated solution are given in the 
following table. It is assumed that the temperature of saturation 
is in each case 25^C, and that the temperature of measurement of 
pH is the temperature listed in the table. 

""C pH ^'C pH pH 



25 3. 557 45 3.547 70 3.580 

30 3.552 50 3.549 80 3.609 

35 3.549 55 3.554 90 3.650 

38 3.548 60 3.560 95 3.674 

40 3.547 

The uncertainty of these pH values is estimated not to exceed 
±0.005 pH units from 25 to 60''C, and ±0.01 units from 70 to 95''C. 
Values of the pH of the 0.01-molal solution at various 
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temperatures are given in the following table: 



Op 




0 p 


pH 


Op 


pH 






20 


3.647 


40 


3.632 


0 


3.711 


25 


•J 9 \J %J iJ 




• V «/ ^ 


5 


3.689 


30 


3.635 


50 


3.639 


10 


3.671 


35 


3.632 


55 


3.644 


15 


3.657 


38 


3.631 


60 


3.651 



The uncertainty in these pH values is estimated not to exceed 
±0.005 pH units. 

POTASSIUM A 0.05-molal solution o£ potassium hydrogen phthalate 
HYDROGEN 

PHTHALATE recommended for the standardization of pH 

equipment. The salt should meet the specifications of the 
American Chemical Society for reagent grade chemical , and should 
be dried for 2 hours at llO'C before use. To prepare the solu- 
tion, weigh 1.012 grams of the dried salt into a 100-m£ volu- 
metric flask. Add sufficient distilled water to dissolve the 
salt, and then fill to the mark with distilled water. Mix the 
solution thoroughly by shaking. The distilled water should have 
a conductivity not exceeding 2 x 10'^ ohm'^cm"^. The solution 
should be protected against evaporation and contamination by 
molds. It should be discarded if mold growth occurs. Values of 
the pH of this standard solution at various temperatures are: 
•C pH "C pH ' pH 



0 4.012 30 4.014 60 4.089 

S 4.005 35 4.023 70 4.12 

10 4.002 40 4.033 80 4.16 

IS 4.001 45 4. 045 90 4.20 

20 4.003 50 4.058 95 4.22 

25 4.008 55 4.073 
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The uncertainty of these pH values is estimated not to exceed 
±0.005 pH units for temperatures from 0 to dO^'C, and ±0.01 unit 
from 70 to gS^'C. 

PHOSPHATE A 0.05-molal phosphate solution (0.025 molal with 
MIXTURE 

respect to both KHaPOn and NaaHPOn) is recommended for 
the calibration of pH equipment. Both salts should meet the spec- 
ifications of the American Chemical Society for reagent grade 
chemicals, and should be dried for 2 hours at 110 to IZO^'C before 
use. To prepare the solution, weigh 0.3388 grans of the dried 
potassium dihydrogen phosphate and 0.3533 grams of the dried 
disodium hydrogen phosphate into a 100-m£ volumetric flask. Dis- 
solve the salt and then fill to the mark with carbon-dioxide-free 
distilled water. The water can be freed from carbon dioxide by 
boiling distilled water for 10 minutes and guarding it with a 
soda-lime tube while cooling. The distilled water should have a 
conductivity no greater than 2 x 10"^ ohm"^cm"\ Although elab- 
orate precautions to prevent contamination of the buffer solution 
with atmospheric carbon dioxide are usually unnecessary, the con- 
tainer should be kept tightly stoppered at all times when a sample 
is not actually being removed. The solution should be replaced 
after a few weeks, or sooner if molds or sediment appear, or if 
it has been exposed repeatedly to air containing carbon dioxide. 
^'C pH ' pH pH 

0 6.981 30 6.850 55 6.832 

5 6.948 35 6.841 60 6.836 

10 6.920 38 6.837 70 6.845 

15 6.897 40 6.835 80 6.859 

20 6.878 45 6.831 90 6.877 

25 6.862 50 6.830 95 6.886 
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The uncertainties in the pH values listed in the preceding table 
are estimated not to exceed ±0.005 pH units from 0 to dO^'C and 
±0.01 units from 70 to 95°C. Minor variations of the order of 
a few thousandths of a pH unit may be expected to occur between 
different lots of phosphate salts. 

BORAX A 0.01-molal solution prepared from NazBifOy-lOHaO is 
recommended for the calibration of pH equipment. The salt should 
meet the specifications of the American Chemical Society for 
reagent-grade chemical. The water content of this salt, stored 
under ordinary conditions, is less than indicated by the formula. 
This does not affect the use of the salt as a pH standard. How- 
ever, the salt must not be dried in an oven before use. To pre- 
pare the solution, crush gently any large lumps of the salt, and 
weigr* 0.380 grams into a 100-m£ volumetric flask. Dissolve and 
then dilute to the mark with carbon-dioxide-free distilled water. 
The water can be freed from carbon dioxide by boiling distilled 
water for 10 minutes and guarding it with a soda-lime tube while 
cooling. The distilled water should have a conductivity no 
greater than 2 x lo"^ ohm'^cm'^. To avoid contamination of the 
buffer solution with atmospheric carbon dioxide, keep the stopper 
in place except when removing a portion of the solution. If 
desired, the solution may be protected with a soda-lime tube. 

pH °C pH pH 



0 9.464 30 9.139 55 8.985 

5 9.395 35 9.102 60 8.962 

10 9.332 38 9.081 70 8.921 

15 9.276 40 9.068 80 8.885 

20 9.225 45 9,038 90 8.850 

25 9.180 50 9.011 95 8.833 
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The uncertainties in the pH values listed in the preceding table 
for standard borax buffer solutions have uncertainties estimated 
not to exceed +0.005 pH units from 0 to 60°C, and +0.01 units 
from 70 to 95°C. 
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